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Abstract 


This work discusses an extension to conventional low-coherence interferometry 
by the introduction of dispersion-encoding. The extension facilitates the mea- 
surement of surface height profiles with sub-nm resolution. The selection of 
a dispersive element for encoding allows for tuning of the axial measurement 
range and resolution of the setup. The approach is theoretically designed and 
implemented for applications such as surface profilometry, the characterization of 
polymeric cross-linking and as a tool for the determination of layer thicknesses in 
thin-film processing. During the characterization of the implemented setup, it was 
shown that an axial measurement range of 79.91m with a resolution of 0.1 nm 
was achievable in the evaluation of surface profiles. Simultaneously, profiles of 
up to 1.5 mm length could be obtained without the need for mechanical scanning. 
This marked a significant improvement in relation to state-of-the-art technologies 
in terms of dynamic range. It was also shown that axial and lateral measurement 
range can be decoupled partially. Additionally, functional parameters such as sur- 
face roughness were characterized with the same tool. The characterization of the 
degree of polymeric cross-linking was performed as a function of the refractive 
index. Here, the refractive index could be acquired in a spatially-resolved manner 
with an index resolution of down to 3.36x 1075. This was achieved by the devel- 
opment of a novel mathematical analysis approach. For the acquisition of layer 
thicknesses of thin-films an advanced setup was developed which could be used to 
characterize the thickness of thin-films and its (flexible) substrate simultaneously. 
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Introduction and Motivation 


The electronics industry with all its branches such as semiconductors, organic- 
electronics and the photovoltaics industry, is continuously growing in terms of its 
economic as well as its technological influence, [1]. This trend is fostered by the 
ongoing integration of various electronic functionalities in fields such as energy gen- 
eration & distribution, transport & mobility as well as in consumer goods. Over three 
decades, electronic and semiconductor products as well as processes were driven by 
MOORESs law, [2]. This paradigm focused nearly exclusively on the miniaturization 
of structures and therefore on the increase of the number of transistors per chip area. 
While this paradigm was directly attributed to an increase in computing power, the 
industry’s scope broadened in recent years. With the introduction of the so-called 
More-than-Moore paradigm, the enhancement of functionality and cross-modality 
integration has become more important, [3]. The implementation of this approach 
has enabled the development of novel products where functionalities of the analog, 
digital as well as of the power domain were combined. Furthermore, sensors and 
other MEMS have been used to enhance these products, [4]. Opposed to traditional 
geometrical scaling of chips this is also known as performance scaling, [5]. The 
implementation of these ideas at the product level has enabled the transformation of 
traditional analog components from board level to package (SiP—system in pack- 
age) or chip level (SoC—system on chip). This includes a significant increase of 
complexity in semiconductor manufacturing, [6, 7]. New and advanced processes 
have been developed to enable three-dimensional production and packaging. Addi- 
tionally, novel concepts for the management of disturbances such as thermal and 
electro-magnetical influences have been developed. All these product and process 
developments drove research and applications of novel materials and material com- 
binations, [8-11]. On the one hand, this enabled bold innovations in power semi- 
conductors regarding performance and energy efficiency which are main drivers 
behind electro-mobility and sustainable energy, [12, 13]. On the other hand, the 
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novel and advanced manufacturing technologies imposed great challenges towards 
quality assurance, reliability and therefore metrology, [14]. According to Leach et 
al., [15], the currently available metrology still lags behind advanced production 
regarding precision, measurement speed and cost scaling. Where processes such as 
lithography could significantly scale on a cost per unit basis, metrology hardly did. 
This was partly due to technological reasons but also partly due to high special- 
ization of metrology approaches with a strong focus on single applications. On the 
technological side, the main challenges lie in the need for systems that deliver large 
measurement ranges with high precision (high-dynamic range) while maintaining 
fast speeds to ensure process-integrated operation. Appropriate metrology has to 
enable the determination of geometrical parameters like critical dimensions or sur- 
face profiles alongside functional parameters such as surface roughness or defect 
detection. A typical case for the requirements of a process-integrated metrology tool 
is given by the stack of a modern solar cell module, [15]. The cell itself consists of a 
number of active and barrier layers, deposited by thin-film processes, which require 
the layer thicknesses to be monitored in order to ensure performance. Interconnects 
or divides such as bus-bars, laser scribes and vias need to be observed regarding 
positions and dimensions. Optics with high aspect ratios on top of the cell increase 
the efficiency and must be measured with respect to their surface profile. Addition- 
ally, polymeric adhesives as well as a polymeric top laminate protect the compound 
against environmental hazards. Monitoring of the degree of cross-linking of these 
materials enables optimal performance over 20 and more years of operation, [16]. 
The main scope of this work is the development of a metrology system which 
is capable of meeting the main requirements posted by this exemplary application 
in terms of resolution, dynamic range, measurement speed and flexibility. More 
precisely, the metrology approach aims to measure surface profiles, thin-film layer 
thickness as well as the degree of cross-linking in polymers. A review of relevant 
research literature paved the way for the development of a novel approach based on 
the principles of low-coherence interferometry. The enhancement of the principle 
by the introduction of dispersion encoding is performed in order to facilitate higher 
resolutions while maintaining large measurement ranges, known as high-dynamic 
range metrology. This approach will be called DE-LCI. With regard to surface 
profilometry, the aim is to cover an axial measurement range of nearly 100um 
while having sub-nm resolution. Additionally, the capabilities of the novel approach 
regarding the characterization of thin-film thickness as well as polymeric cross- 
linking will be assessed. This work evaluates the necessary developments of the 
setup and the analysis algorithms facilitating all three possible measurement modes. 
The design and implementation of the novel approach will furthermore focus 
on the robustness for thermal and mechanical influences to maintain the resolution 
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during long measurements. Another design criterion is the preparation for process 
integration which can be translated to acquisition and data processing speed. This 
work will demonstrate efforts to capture large sections of a sample comparatively 
fast to conventional methods. Theoretical limits of measurement range and reso- 
lution will be investigated on a more detailed level. Also, the occurrence of these 
limits will be assessed under practical conditions. In particular, it will be investi- 
gated which transients influence the axial resolution and measurement range and 
which technological measures can be taken to partially decouple axial and lateral 
measurement ranges. Furthermore, this work will evaluate techniques to gather two- 
and three-dimensional information of surfaces and bulk materials without the need 
for mechanical scanning to reduce influences from movement and increase mea- 
surement speed. 

In terms of the characterization of polymeric cross-linking, it is investigated 
how the wavelength-dependent refractive index can be used as a measure for cross- 
linking. Consequently, the resolution of the DE-LCI approach will be examined with 
regard to the refractive index. An exploration of the capabilities to resolve nm-sized 
layers of thin-film materials will conclude the work. 
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2.1 Profilometry 


The characterization and measurement of surface profiles is one of the most basic 
metrology tasks in industrial manufacturing. What started with mechanical stylus 
profilometers has developed with ultrasonic transducers towards optical instruments. 
These are capable of appropriate resolutions to enable nanometrology, [14]. As sur- 
face metrology is well established, basic terms and parameters are defined in cor- 
responding norms such as ISO 4288 (assessment of surface texture) and ISO 4287 
(terms, definitions and surface texture parameters), [17, 18]. Established optical 
instruments are classified as areal integrating, line profiling and areal topography 
instruments in ISO 25187, [19]. According to Leach et al. [14], sensors with the 
capability to record areal information about surface shape, waviness and roughness 
simultaneously are the most desirable in nanometrology. These sensors are usually 
defined as 3D-sensors where the topography information is gathered as the local 
distance z at a specific coordinate (x, y). Apart from the classification mentioned 
in [19], numerous variants exist, [20]. On the most basic level, the common denom- 
inator for optical approaches is that all of them follow a simple interaction model 
where light from a light source interacts with an object of interest. The result of this 
interaction can be observed separately on an observation plane. The properties of 
the light source such as intensity, polarization state, temporal behavior, coherence 
and spectral behavior can be utilized to encode information about the object sur- 
face topography. Depending on the surface structure of the object of interest and 
the design of the observation method, different limitations or artifacts occur. One 
example is the occurrence of so-called speckles on optical rough surfaces under 
coherent illumination. While most metrology approaches suffer from this effect in 
terms of a decreased signal-to-noise ratio (SNR), some techniques such as coher- 
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ence scanning interferometry or speckle interferometry make use of this effect. 
Production integrated metrology becomes increasingly important in areas where 
the critical dimensions decrease from the um- to the nm-range. Industrial sectors 
such as MEMS, power semiconductors or photovoltaics require a continuous, pre- 
cise monitoring of production to maintain quality, [14]. A variety of measurement 
technologies to characterize e.g. surface roughness, topography or film thickness 
have been reported, [21]. Apart from specific measurement machines with pm res- 
olution for research purposes which require high efforts in building, calibration and 
maintenance, a broad range of technologies is known which is more applicable as 
production-accompanying tool, [22]. The following section is intended to give an 
overview of these approaches. 

Haussler and Ettl classified optical 3D sensors in [21] mainly according to their 
limitations and their dominant noise sources. Following that, four main classes of 
sensors can be associated: 


Type I sensors (e.g. laser triangulation) 

Type II sensors (e.g. coherence scanning interferometry) 
Type III sensors (e.g. phase-shifting interferometry) 
Type IV sensors (e.g. deflectometry) 


As Type I sensors mostly incorporate technologies which measure a lateral perspec- 
tive shift on local details, the uncertainties are determined by the uncertainty of this 
shift measurement. This means that the uncertainty scales with the inverse square 
of the distance to the detail, [21]. Techniques include laser triangulation, fringe 
projection, confocal microscopy and others. For Type II sensors such as coherence 
scanning interferometry, the measurement relies not on the detection of the phase of 
the signal but on the correlation of single speckles reflected from a rough surface. 
Therefore the statistical noise, mainly composed of the standard deviation of the 
object surface, is the primary source of noise, [23, 24]. In contrast, Type III sensors 
involve classical interferometric technologies such as phase-shifting interferometry 
[25] and digital holographic microscopy [26]. These are phase measuring technolo- 
gies which enable sub-nm axial resolution and are only limited to photon noise. 
Furthermore, [21] separates deflectometry as a Type IV sensor. With limitations 
of photon noise and developments towards nanoscale accuracy, this technology is 
promising for industrial applications, [14]. 

This section analyses important technologies from every sensor class as well 
as atomic force microscopy (AFM) as one important, high-resolution, stylus-based 
technology in order to evaluate open questions in nano-profilometry. All of them 
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have been widely qualified for the measurement of nm-scale surface parameters in 
a lab environment, [27, 28]. 


2.1.1 Atomic Force Microscopy 


AFM, although not an optical metrology approach, has become an important and 
highly precise tool in nano-engineering and science, [29-31]. It is a technological 
development which goes beyond the capabilities of classical stylus instruments in 
terms of its resolution. It provides an axial resolution down to 0.01 nm while the 
lateral resolution can be in the range of 0.1-10nm, [32]. These features made it 
exceptionally important in areas such as cell biology and molecular sciences where it 
enables the precise study of molecular processes and interactions, [33, 34]. Industrial 
and technical applications make use of it for the characterization of nano-structured 
materials and features, [35-37]. AFM can deliver interesting additional features 
like force sensing and surface modification, [38]. Several developments such as 
multi-cantilever arrangements and resolution-enhancing methods also have enabled 
the technology to investigate cm?-large structures with high aspect ratios. In this 
context, aspect ratios of a axial measurement range of 15 um with a deviation of + 60 
nm could be achieved while the lateral measurement range was 15x 100 um?, [39]. 
Using an array of cantilevers in a 32x1 chip design, Minne et al. [40] were able to 
scan large areas of 2000 x 2000 um? within 30 minutes while the axial measurement 
range was 2.5um at a resolution of 3.5nm. Most of these developments have the 
need for complex setups and measurements last several minutes. Additionally, AFM 
is sensitive to systematic errors and contamination of samples. These circumstances 
prevent using AFM-based approaches for production accompanying tasks. 


2.1.2 Confocal Laser Scanning Microscopy 


Another established method for the characterization of surface profiles is confocal 
laser scanning microscopy (CLSM), [41, 42]. By aligning the focal plane of the 
illumination of a single point on a sample with the imaging plane of that same point 
on a pinhole, information is solely captured from this point, [43]. By scanning a 
depth range, the height profile of a sample can be acquired with a high SNR even 
if the surface is strongly scattering. As scanning has to be performed in the lateral 
dimension as well, the acquisition of a full areal height map requires a long time, 
meaning several minutes. In order to achieve a high axial resolution, the so-called 
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auto-focus method is applied. In this method, the focal plane is adjusted for the 
highest intensity for a given point. The method is repeated in a given axial range 
of the sample while the position of the focal points is recorded correspondingly, 
[44]. A faster, yet less precise, option is the so-called intensity method. This method 
correlates a measured intensity to a calibrated height-intensity curve. According 
to [44], this results in measurement times in the order of 1/50 compared to the 
auto-focus method. Also, due to the limited linearity of the calibration curve, the 
measurement range and precision is also significantly reduced. Furthermore, the 
calibration curve has to be recaptured for every new sample material in combination 
with the microscope objective used for the investigation. 

Depending on the translation stages and magnification, a lateral resolution of 
(100-500) nm is typically achieved to capture surface parameters such as roughness, 
waviness and form error in an areal fashion, [45]. Larger lateral measurement ranges 
such as 1.3x 1.3mm are usually captured with a lower resolution of e.g. 10 um in 
order to speed up the measurement process. Buajarern et al., [42], have shown that 
the axial resolution is not only limited by the resolution of the translation stages 
used to perform the auto-focus determination, but even more by the depth of focus 
/ resolution of the microscope objective used for imaging. Using high numerical 
aperture (NA), large magnification objectives, the axial resolution is optimal but 
restricted to about 150nm. 

Different works have shown the applicability of CLSM to characterizations in 
biology, for dental as well as engineering materials [46-48]. Besides the need for 
mechanical scanning and the interdependence of measurement range and resolu- 
tion, a major disadvantage of CLSM is the high effort necessary for adjustment of 
parameters in order to achieve optimal results. Tomovich et al., [49], have demon- 
strated that the adjustment of parameters is highly dependent on the reflective and 
structural properties of the sample. From this perspective, CLSM is a versatile tool 
in research labs but is limited in the usage for production accompanying tasks. 

In addition to classical confocal microscopy, several developments have been 
made in order to circumvent some of the drawbacks such as the need for scanning. 
Chromatic confocal microscopy is one approach. It particularly targets scanning 
along the z-dimension. While dispersion in the optical system leads to different 
foci for each wavelength, height information can be decoded in spectral variations, 
[41, 50]. In terms of measurement range and resolution, both — the spectral range 
and spectral intensity stability—are important. Typical light sources include halo- 
gen and xenon lamps as well as LED sources. Other works have demonstrated the 
usability of supercontinuum sources as well, [51]. Using the linearized wavelength 
calibration, the calculation of height can be performed in different ways. The fastest 
and computational most efficient technique, yet the one with the lowest resolution, 
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is to determine the maximum intensity of every spectral line. The highest computa- 
tional effort is necessary when fitting the intensity distribution with e.g. a Gaussian 
approach, though these approaches enable the highest possible resolutions. A possi- 
ble compromise between speed and accuracy is the so-called bary-center calculation 
to estimate the height, [21]. The limitations regarding the lateral resolution and the 
resolution of slopes are very similar to other microscopy techniques. Depending on 
the light sources, the NA of the objective and the size of the pinhole, typical spot 
sizes are 5—10 um for vertical ranges of < 1mm and 10-30 um for vertical ranges 
> 1mm, [21] can be achieved. The resolution of slopes is usually defined by the 
half aperture angle where the range is between + 18° to + 44° for NAs of 0.3 to 0.7 
respectively. 

The extension of this method in order to capture areal information normally relies 
on scanning approaches. Some problems occur due to dynamic stitching errors and 
spherical aberrations from scanning, [52, 53]. 

The method can also be used to measure the profile and thickness of multiple, 
semi-transparent layers. In this case, the analysis has to take into account the refrac- 
tive index of the transparent material at the different wavelengths. Faster approaches 
utilize multi-probe setups in order to parallelize data acquisition. Nonetheless, they 
are strongly limited regarding the lateral resolution, [54, 55]. 


2.1.3 Digital Holographic Microscopy 


Digital holographic microscopy (DHM) is a fast, robust and full-frame-capturing 
technology to evaluate the surface topography of samples. While surface data from 
an interferometric approach is captured by a camera as full-field information (spa- 
tially, intensity, phase) within the integration time of the device, all aberration cor- 
rection, exact focusing and surface reconstruction can be done digitally within post- 
production. Stroboscopic illumination can be used additionally to speed up acqui- 
sition from conventional frame rates such as 20-60 fps to 25 MHz, [21]. The axial 
accuracy is dependent on the calibration of the light sources wavelength/frequency 
which leads to achievable axial resolutions of 0.1 nm with stabilized laser sources, 
[21]. One of the limiting factors is the axial measurement range which can be 
extended by the use of multiple wavelengths to several ums while keeping a sub- 
nm resolution, [56]. Furthermore, the combination with other techniques, such as 
reflectometry, enables the analysis of parameters such as layer thickness, refractive 
index profiles and topography of multilayered structures, [57]. The signal formation 
and analysis is generally described in two steps where one is the data acquisition and 
the other is the so-called reconstruction, [26]. A typical holographic setup contains 
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e.g. aMach-Zehnder interferometer as its main component. In this setup, the object 
wave, which is either reflected or transmitted from the object, is imaged close to the 
camera plane. The reference wave interferes with the object wave and its intensity 
can be detected and analyzed. The subsequent reconstruction is performed in two 
steps. First, the complex wave is reconstructed in order to separate the real and vir- 
tual hologram from the zeroth-order signal. This can be performed in hardware by 
introducing defined phase jumps (in inline-holography) or by introducing an angle 
between the object and the reference wave (in off-axis holography). Furthermore, 
the separation of the zeroth-order term [58] or the twin-image term [59] can be 
done by FFT approaches in software. Afterwards, the filtered hologram is digitally 
illuminated by a reference wave. During a second process step, numerical propaga- 
tion is performed by approaches like the single Fourier approach [60], the angular 
spectrum approach [61] and the convolution approach [62]. In typical implementa- 
tions, two approaches to define numerical lenses for the correction of higher order 
aberrations: 


e recording of a physical reference hologram on a plane sample and by transmission 
through air, [60, 62] 
e fitting of Zernike polynomials in an assumed flat area of a measured sample, 


[63-65]. 


It is quite common to use combinations of both solutions in order to achieve optimal 
results. The reference wave is used to correct all phase and tilt aberrations as well 
as for curvature of waves and de-focus correction. In particular, this technique has 
the advantage that only a single hologram needs to be acquired while the correct 
focal point can be tuned by numerical propagation. This can also be used to increase 
the depth of field, [66]. Furthermore, it indicates that DHM is a full-field metrology 
where all relevant information is captured in one image acquisition. 

A common problem in DHM is 2r/phase ambiguity as it limits the usable mea- 
surement range. An efficient method to overcome this problem is the so-called 
multiple-wavelength holography, [67, 68]. In this approach, data is captured with 
slightly different wavelengths which allows the numerical construction of a new, 
artificial wavelength which can be used to ‘construct’ unambiguous measurement 
ranges. DHM is also capable of recording dynamic events which occur in a repeti- 
tive fashion. For this purpose, the events of interest are observed using stroboscopic 
techniques. It could be shown that an acquisition of events with a repetition fre- 
quency of up to 25 MHz is possible, [69]. Extensions to DHM in order to evaluate 
functional parameters such as reflectometry have been reported, [21]. However, 
complex wavefront reconstructions or fitting procedures are necessary to separate 
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the information of different layers. Currently, this prevents the use of this approach 
in industry. 


2.1.4 Phase-shifting Interferometry 


Phase-shifting interferometry (PSD, is an more elaborate variant of traditional 
(high coherence) interferometry for the determination of surface information, [25]. 
Depending on the light source, the detector and the optical setup, lateral resolution 
of 2 um and axial resolutions up to Inm are reachable, [21]. PSI enhances classical 
interferometry by introducing multiple controlled phase shifts in order to determine 
the phase-dependent surface height of an object, [70] . Furthermore, it overcomes 
the directional ambiguity of interferometric signals by generating quadrature sig- 
nals, since back-reflected light from a surface in combination with light from the 
reference arm generates an areal interferometric fringe pattern, [21]. Due to some 
simplifications data analysis of the generated quadrature signals can be accelerated 
by fitting sine and cosine signal components to determine the phase with high accu- 
racy. A common approach for the data analysis of these signals is the so-called 
auto-correlation or synchronous detection, [21]. Both, sine and cosine components 
are described as integrals of a full cycle of the phase shift. Depending on the number 
of introduced phase shifts, various algorithms can be used, [25]. Furthermore, the 
use of non-linear e.g. sinusoidal manipulation of the phase shifts has been reported 
to be accurate and computational more efficient than linear manipulation, [71]. In 
addition to temporal PSI, other works have described alternative methods to encode 
and detect phase-shifts such as polarization encoding, [72, 73]. 

A common problem in interferometry occurring in PSI is the determination of 
phase using an arctan function which is only possible in the range of 0-27. Sur- 
face heights corresponding to multiples of 27 will therefore be calculated within 
the given range resulting in phase jumps, known as a wrapped phase. Several algo- 
rithms are currently being researched, [74-77]. More recent approaches on phase- 
shift interferometry extend the standard setup with multiple sensors, [73, 74, 78]. 
By capturing interference signals at different polarizations, phase angles and wave- 
lengths, an axial resolution of 2nm on a 3 um height step in an area of 0.5x0.5 um 
could be achieved. Due to the need for multiple channels, high efforts have to be 
made to synchronize and handle the data of up to 6 cameras. The evaluation of the 
phase shift of the total interference contrast (TIC) is an alternative phase-shifting 
technique which is capable of providing sub-nm axial resolution in a full-field man- 
ner, [79]. As it is not based on coherence properties of the light source it is reliant 
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on the knowledge of an material model of the tested sample. Therefore, the method 
is limited to mainly laboratory usage, [80]. 


2.1.5 Coherence Scanning Interferometry 


Belonging to so-called type II sensors, coherence scanning interferometry (CSI) 
distinctively deviates from classical interferometry. This approach utilizes light of 
low temporal coherence not to measure the phase of the signal, but to estimate 
the correlogram of individual speckles on a surface, [21]. The technique tradition- 
ally scans one interferometer arm in order to evaluate the signal contrast which is 
equivalent to the envelope of the correlogram, [70]. The coherence length of the 
light source is the main limitation as it restricts the possible measurement range. 
Thus, only height variations within the coherence length can be detected continu- 
ously. This can be an issue on very rough surfaces or steep slopes, [81]. Usually 
setups are designed to capture the full lateral information using a camera which 
detects the correlogram envelope at every point P(x, y), while a scanner ensures 
the axial movement in the height range, [82, 83]. While the typical scanning range 
is (10-200) um with piezo-based scanners and > 1 mm for mechanical scanners, 
the interferometer part of the system is most commonly designed as an microscope 
objective of the MICHELSON, MIRAU or LINNIK type. The approach is also known 
by the names of coherence radar, coherence scanning, white-light interferometry, 
vertical scanning interferometry and others, [70]. 

According to [84], the simplest model for the description of the signal in CSI is an 
incoherent superposition. The signal shape in this case is determined by the spectral 
and spatial distribution of light in the pupil plane. Using this model, several boundary 
cases can be found. For example the properties of a system with low NA (e.g. 0.2), 
narrow bandwidth (e.g. AA=100nm) will basically be dominated by the spectral 
distribution of the light source with high fringe visibility at the zero scan position. 
If the system has a high NA (e.g. 0.6) and narrow bandwidth (e.g. AA =20nm), in 
contrast, the coherence properties are dominated by the spatial distribution in the 
pupil plane, hence the focusing. In real systems both effects will be present, whereas 
it can be stated that in low NA systems the spectral contribution and in high NA 
systems the spatial contribution of the light source primarily influences the signal. 

One of the distinguishing features of CSIis the lack of a 27 ambiguity. This makes 
it well suited for measurements of rough surfaces and steep slopes. In general, it can 
be said that the slowly modulated envelope of the interference fringes is an indicator 
of the signal strength as well as a measure for the surface profile. But it has to be 
clear that the peak of the envelope is not equal to the maximum height of the sample. 
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It represents only the point where the optical path difference (OPD) equals zero. 
More precisely, it is the point where the dispersion affected path length difference 
in the setup equals zero. Therefore, the approximation of the envelope signal is only 
an idealization as in reality dispersion effects and scan dependent distortions have 
to be taken into account. Signal processing in CSI is usually performed in multiple 
steps, where the start of the scan position is determined by the maximum signal 
strength or fringe visibility. The scan range is determined by sectioning the images 
according to the expected height of the sample. Several methods to perform the 
envelope detection including peak finding of the maximum have evolved over the 
years. While [85] have shown the application of digital filters and an demodulation 
technique, other groups such as [86] and [87] have shown how to remove the carrier 
signal to separate the envelope. According to approaches of [88] and [89], who 
performed a centroid determination of the square of the signal derivative, a higher 
resistance to noise on the signal is achievable. 

Noise, e.g. from optical aberrations, the scanning process, diffraction or vibration 
can significantly influence the sensitivity ofthe technique. Therefore, more advanced 
signal processing algorithms use the result of the envelope detection only as arough 
estimate of the fringe order and call the result a first topography map. Based on this, 
a fine estimation using a phase analysis of the fringes results in a phase map. This 
method increases the accuracy significantly although other problems such as 27 
ambiguities on thin film structures, surface roughness or sharp edges occur during 
phase analysis. Some algorithms take this into account, [90, 91]. 

More sophisticated approaches developed methods for the simultaneous enve- 
lope detection and phase analysis by the correlation of the intensity data with a 
complex kernel. These kernels have been derived fromPSI algorithms and use dif- 
ferent approaches such as wavelet techniques or least-squares fitting, [88, 92-94]. 
Furthermore, recent works have also demonstrated that sub-nm axial resolution is 
possible through advanced signal modeling techniques, although there remain lim- 
itations due to scanning and the relatively small axial measurement range, [95, 96]. 
Coherence Scanning Interferometry is well suited for the determination of the sur- 
face topography on rough surfaces, [97]. One major disadvantage is the necessity 
for scanning along one physical dimension of a sample as this causes errors and 
slows down the measurement procedure. 


2.1.6 Low-coherence Interferometry 


One of the most advanced technologies for high resolution surface profilometry 
and tomography is low-coherence interferometry (LCI). Optical-coherence tomog- 
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raphy (OCT) is a commonly known implementation, [98, 99]. Although OCT was 
initially used in medicine, developing it as a tool for industrial purposes has become 
increasingly important, [100-103] 

The range of different OCT approaches reaches from time-domain optical- 
coherence tomography (TD-OCT) and spectral-domain optical-coherence tomog- 
raphy (SD-OCT) to frequency-domain optical-coherence tomography (FD-OCT), 
[104]. While most approaches are implemented as point sensors with the possibility 
to mechanically scan a surface, they vary in speed, measurement range and accuracy. 
More fundamentally, the accuracy and measurement range in the axial dimension 
are defined by the coherence length /, of the light source used, [105], 
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where àc denotes the center wavelength of the light source and AA its spectral range, 
assuming that the light source has a Gaussian-shaped spectrum. The typical axial 
measurement range, using broadband light sources, is a few hundred um with a 
resolution of about 1 um, [106]. These values are usually fixed within the design of 
the specific setup and can not be adjusted during operation. In order to achieve higher 
axial measurement ranges as well as a higher resolution some hardware problems 
must be solved. The most notable influences on the resolution are environmental 
disturbances on the interferometer arms as well as the repeatability of the scanning 
system used to gather areal information, [107]. An extension of the measurement 
range in TD-OCT is achieved by the introduction of multiple reflecting surfaces in 
the reference arm, [108]. This modification enables the measurement of a large axial 
range while keeping actual mechanical scanning to a minimum. Another approach 
merges both reference and sample arm in a FD-OCT configuration, [107]. In this 
approach, a reflective surface in the vicinity of the sample is used as a reference. 
Environmental changes have the same influence on both optical paths. This hardware 
adaption in combination with algorithmic frequency and phase evaluation leads to 
an accuracy of 0.1 nm in a measurement range of about 3 mm in the axial dimension. 

Other approaches aiming to extend the measurement range and to increase the 
resolution rely on dispersion compensation and are known by the term dispersion- 
encoded full-range OCT (DEFR-OCT), [109, 110]. The authors applied a numerical 
dispersion compensation to remove complex conjugates and therefore extend the 
measurement range. Another advantage is that the point of highest sensitivity of the 
setup is shifted to the center of the measurement range. As the initial algorithm relied 
on multiple iterations of Fourier transforms, new approaches have been developed 
to speed up processing times. Some of them introduce artificial dispersion into the 
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setup which can be compensated by the algorithm with only one Fourier transform 
and aconvolution, [111]. Current works applying these ideas show high accuracy in 
the tomographic analysis of nano-structured domain walls in ferro-electric media, 
[103]. 

In contrast to medical applications, industrial measurement tasks often demand 
the characterization of relatively large areas (mm? instead of um?). Scanning of 
samples is therefore necessary, but introduces issues regarding the accuracy and 
repeatability of results, [112]. Approaches to overcome these limitations are based 
on a setup incorporating a camera and a high resolution translation stage; summa- 
rized by the term full-field OCT (FF-OCT), [113-116]. In this particular approach 
scanning is still required in the reference arm in order to observe changes in the 
interference data due to different path lengths. 

More recently, full-field approaches that avoid scanning altogether have been 
developed. Based on a hyperspectral imager and an etalon, Zhu et al., [117] have 
shown high precision surface measurements. While one dimension of the sample 
is encoded in the bandwidth of the etalon-generated spectral slices, the second 
dimension is directly imaged on the camera. The surface profile is calculated from 
the phase information of the interference fringes. The resolution and measurement 
range of this setup is predominantly defined by the wavelength spacing of the etalon 
as well as by the resolution of the spectrometer. Both tuning options are opposing 
each other. 

A more advanced method substitutes the etalon for a microlens array in order 
to decode one areal dimension on the hyperspectral imager which increases the 
measurement range and light efficiency, [118]. The setup was capable of measuring 
an axial range of up to 880 um with an resolution of 0.49 um while having a lateral 
measurement area of 3.5x 3.5mm. Further developments by the same group have 
demonstrated the possibility to acquire 2500 independent probing points which 
increased the measurement range, light efficiency as well as the tilt angle acceptance, 
[119]. One drawback is that the method is only able to make use of about 50% of 
the detector size to image uum?-sized samples which decreases its lateral resolution. 
The axial measurement range was about 825 um where a resolution of 6nm was 
achieved. According to these works, the dynamic range (DR) is defined as the 
inverse ratio of the resolution to the measurement range. The analyzed works of 
recent, high-dynamic range approaches to LCI and OCT demonstrate a progression 
in technology and dynamic range, Tab. 2.1. 

In direct comparison, it can be deduced that one aim is to develop a full-field, areal 
approach to surface profilometry which is capable of measuring large distances with 
sub-nm resolution as Koch et al. have shown in a point sensor. Recent works such 
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Table 2.1 Comparison of current LCI approaches regarding measurement range, resolution 
and dynamic range 


Authors Axial Axial resolution | DR Remarks 
measurement 
range 
Pavlicek / 900 um 0.05 um 18000 point sensor, 
Häußler ([120], scanning 
2005) necessary 
E. Koch et al. 3mm 0.1 nm 3 x 107 point sensor, 
([107], 2005) scanning 
necessary 
Zhu et al. 250 um 0.1 um 2500 etalon-based 
([117], 2012) hyper-spectral 
imaging 
Ruiz et al. 880 um 0.49 um 1871 microlens-based 
([118], 2017) hyper-spectral 
imaging 
Reichold etal. |825 um 6nm 137 500 microlens-based 
([119], 2019) hyper-spectral 
imaging 


as of Reichold et al. have already demonstrated the potential of spectrally-encoded 
interferometry but still lack sub-nm resolution. 

The current work aims to develop a system which is able to detect surface profiles 
with sub-nm resolution and a high-dynamic range in the axial as well as lateral 
dimension incorporating sub-nm resolution. 


2.2 Polymer Cross-linking Characterization 


The control of mechanical, electrical and optical properties of polymers during fab- 
rication is necessary to ensure their performance, [121]. As cross-linking is a crucial 
process step in order to optimize properties and fabrication parameters, it is nec- 
essary to monitor its degree, [122]. Long-term mechanical resistance, temperature 
stability as well as functional parameters such as refractive index are adjusted with 
the cross-linking process, [123-125]. 

Interesting industrial applications of polymers are centered around lithographic 
processing such as resist coatings and optical waveguides, [126, 127]. Polymer- 
based optical waveguides are usually processed by patterning (photoresist-based or 
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direct lithography), soft lithography or printing techniques [128, 129] in order to 
achieve defined cross-linking and refractive index differences. From the analysis 
of reaction kinetics of polymers it is known that the density of a material as well 
as its refractive index changes during cross-linking, [123]. This relationship can be 
described by the Lorentz-Lorenz equation, [130, 131], 


u (n? = 1)My 


Mh m an. 22 
(n? + 2)p =) 


where Rm represents the molar refractivity, n the refractive index, Mọ the mass- 
averaged molar mass and p the density. Although it is known that the Lorentz-Lorenz 
relation is only an approximation, it has been proven applicable to a variety of poly- 
mers and was correlated with other methods such as hardness measurements, [131, 
132]. For their fabrication, optically-cured polymers have to fulfill several require- 
ments such as optical transparency and chemical as well as thermal stability [133, 
134]. Advancing from conventional thermoplastics such as polymethyl methacry- 
late, polystyrene, polycarbonate and polyurethane, research has been geared towards 
the development of new polymers which exhibit lower absorption losses and higher 
stability, [128]. Promising classes of polymers are halogenated polyacrylates [135], 
fluorinated polyimides [136] or polysiloxanes [137]. In particular, applications such 
as polymeric waveguides or direct laser writing on wafers make use of this effect to 
generate functional properties with refractive index changes of about 10~7, [138]. 
Zukauskas et al. applied this effect to generate gradient-index lens elements with a 
size of 50x 50x 10 um?, [139]. 

In order to characterize these functional properties alongside with the degree 
of cross-linking and their spatial distribution, different metrology approaches are 
known from literature. 


2.2.1 Soxhlet-type Extraction 


A very common method to determine the degree of cross-linking is the Soxhlet-type 
extraction, [140]. For that purpose a sample is exposed to a solvent, typically a xylene 
isomer, in which non-cross-linked material will dissolve after a few hours. After 
drying, the degree of cross-linking can be calculated from the respective weights 
of cross-linked and non-cross-linked material. The result is very precise whereas 
the process is time-consuming, destructive and not spatially resolving. Oreski et al. 
found that the time for the extraction is at least 18 hrs while drying takes another 
24 hrs, [141]. Furthermore, Hirschl et al. found that the extraction time and other 
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process parameters can have a huge influence on the repeatability of the measured 
degree of cross-linking, especially in weakly cross-linked samples. They determined 
that the repeatability ranges from 2-4 %, [142]. 


2.2.2 Differential Scanning Caliometry 


A method to determine cross-linking in polymers as a measure of thermal fea- 
tures such as heat-flow, melting point and reaction enthalpy is differential scanning 
caliometry (DSC), [121, 143, 144]. Similar to chemical-based methods, this method 
does not allow spatially resolved measurements, works destructively and is time- 
consuming. A typical measurement cycle in the so called dual-run mode takes 
2x45 min during which a defined heating profile is applied, [141]. A comparative 
study has shown that different approaches for referencing the measurements to other 
methods might apply and also that errors in the repeatability can be 10% and larger, 
[145]. Especially weakly cross-linked samples require slower heating profiles, hence 
longer measurement times, and produce inherently larger errors, [143, 146]. 


2.2.3 Dynamic Mechanical Analysis 


Dynamic mechanical analysis (DMA) is a well established laboratory approach to 
characterize thermal and mechanical properties of materials. It is especially well 
suited for visco-elastic materials such as polymers, [121]. During analysis a small 
dynamic force (e.g. in tension, compression, bending or shear mode) is applied to 
a well defined sample, typically in a sinusoidal fashion, [147]. The measurement 
time, the temperature as well as the oscillation frequency are typically variable 
parameters dependent on the sample types to be measured, [148]. In response to 
the introduced stress ø, the samples strain £ is measured dynamically. Here, these 
values are measured as the amplitudes of the sinusoidal signals. Furthermore, the 
phase difference between both signals is measured as ®. From these, the complex 
modulus of the sample E* = © can be calculated. It holds information on the real 
and imaginary components which are used to characterize the elastic properties of 
a material. The real component E’ is known as the storage modulus which is pro- 
portional to the elastic deformation work stored by the sample during deformation. 
The imaginary component, known as the loss modulus E”, is a measure for the 
thermo-mechanical loss of deformation work due to internal friction among other 
effects. Both parameters can be referred to the phase difference ® by 
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E' = |E*|- cos ® (2.3) 
E" = |E*|-sin® (2.4) 


From this relation, it becomes clear that the phase difference between excitation 
and response is crucial for the understanding of certain material properties. While 
fully elastic materials, such as steel, typically show a phase difference of ®=0, 
fully viscous materials show a phase difference close to ® = 90°. In the case of 
visco-elastic polymers any value in between may occur. Another important measure 
in DMA is the determination of the glass-transition temperature Tg. It describes the 
transition from a state where molecular networks are stiff and allow only elastic 
deformations to a state where non-elastic deformations are also possible. The Tg 
can be calculated with a variety of approaches, some of which are comparable to 
DSC measurements, [121]. 

In order to calculate the applied stress to a specific sample, the knowledge of its 
geometry is important. Usually standard sample sizes are fabricated for the purpose 
of DMA measurements, [148]. The investigation of fabricated products or compo- 
nents is not possible. As the dimensional parameters depend on the kind of stress 
applied, their input can be significant. Any error in the determination of the samples 
dimension will have a strong influence on the resulting stress on the sample as well 
as on comparability between samples. 

The time for a measurement is very significantly dependent on the chosen exper- 
imental parameters. Usually the temperature range as well as the heating rate are 
subject to changes while the excitation frequency is kept fixed. Typical investiga- 
tions of polymeric materials work in ranges from -75 - 150°C with a heating rate of 
2 K/min, [149] to ranges of -150- 200°C at a heating rate of 1 K/min, [143]. 

One of the main advantages of DMA is the ability to directly measure the degree 
of cross-linking, [145]. Typically the degree of cross-linking is determined as a 
weighted ratio of the storage modulus, the sample density and the monomer molec- 
ular weight, [150]. Other works propose the calculation as a quotient of the loga- 
rithmic storage modulus of cross-linked and non-cross-linked material, [151]. 

DMA is, amongst others, a standard testing procedure in material development 
for novel polymers, fabrication technologies or the understanding of cross-linking 
mechanisms, [152-154]. With its high accuracy, the capability to determine the 
degree of cross-linking directly and the long measurement times, DMA is a pre- 
dominantly laboratory-based method. Applications in production accompanying 
tasks are not established. 
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2.2.4 Spectroscopy-based Methods 


Non-destructive measurements of cross-linking can be obtained by using optical 
metrology such as Raman spectroscopy, [155] or luminescence spectroscopy, [156]. 
The analysis of spectral features of reactional groups and bonds or the lumines- 
cence intensity of characteristic peaks can be utilized to calculate the degree of 
cross-linking. Recent works have shown that these technologies are able to charac- 
terize cross-linking of coatings on solar cells. In a comparative study, Hirschl and 
co-workers, [157], have demonstrated that Raman spectroscopy gives comparable 
results to classical methods like Soxhlet-extraction. Although it has to be noted that 
the measured errors of the degree of cross-linking were up to 15%, especially for 
samples with weak cross-linking. Furthermore, acquisition times for Raman spec- 
tra depend very much on the SNR of relevant spectral intensity peaks and hence 
require a large amount of averaging. Recent studies report acquisition times for 
single-point measurements between 50-100 s, [139, 157]. Peike et al. [155] point 
out that Raman analysis is very material-specific and can be complex with differ- 
ent peaks overlaying each other. Additionally, they found that the SNR decreases 
with peaks at higher wavelengths as Ipeqx ~ 1/4 . This can be critical for weakly 
cross-linked material or materials with a low number of reactional groups. 

A recent work by Schlothauer et al. has qualified luminescence spectroscopy as a 
tool for cross-linking characterization with an accuracy of 4-6%, [156]. However, 
the method requires a large amount of averaged spectra in a point-by-point scanning 
fashion. Acquisition times for a 16x 16cm? were about 80 minutes. 


2.2.5 Low-coherence Interferometry and Other Optical 
Methods 


As refractometry is a well established method to measure refractive indices, it is 
also suited to evaluate cross-linking in polymers, [158]. It has been used to directly 
determine cross-linking progress during curing within thermoset polymers used as a 
matrix material for composites, [159]. For the purpose of these examinations, optical 
fibers are declad over some probing area and integrated in a to-be cured polymeric 
matrix in order to measure the change in refractive index in terms of a change 
in transmission intensity, [158]. Another common approach is the integration of 
cleaved optical fibers in a polymeric matrix where the Fresnel reflection is measured 
in relation to the change in refractive index, [160, 161]. Both approaches are well 
suited for the characterization of thermoset polymers as the change in refractive 
index is rather high during cross-linking with An=2-6 x 1077. 
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Classical optical-coherence tomography has been used to examine structural 
defects such as bubbles or phase separation during cross-linking by scanning a sam- 
ple within a few seconds, [162]. Other interferometric techniques such as spectrally- 
resolved white-light interferometry, frequency domain interferometry or digital 
holographic interferometry have been utilized to measure refractive indices with 
accuracies in the range of 1075-1076, [163, 164] as well as mechanical deforma- 
tions on the nm-scale in material and biomedical engineering, [165-167]. Methods 
based on phase-sensitive OCT have been used to characterize photo-elasticity on 
polymeric composite materials and might also be suitable for cross-linking char- 
acterization, [168]. More recent works made use of the combination of LCI and 
confocal approaches in order to measure refractive index of transparent media from 
a distance, [169]. 


2.2.6 Spatially-resolved Approaches 


None of the technologies for cross-linking characterization presented so far are 
inherently spatially resolving. Spatial resolution for chemical methods such as DSC 
or DMA is usually realized by cutting samples into defined sub-samples, [170, 171]. 
Especially spectroscopic technologies such as Raman spectroscopy gain spatial res- 
olution by scanning over a sample. Due to the need for integration over multiple 
spectra at each point, the measurement times of a sample having 330 x 150 prob- 
ing points can be as long as 82 minutes, [170]. Although this technology can be 
used to cover an area of a few cm”, the typical lateral resolution limit is about 2 
um, [172]. Consequently, other groups have shown a significant increase in lateral 
resolution down to 35 nm by combining Fourier-transform infrared spectroscopy 
(FTIR) with AFM. This approach increases the measurement time significantly so 
that an area of 1 x 1.5 um? was measured in 7.4 hrs, [173, 174]. None of these tech- 
nologies have been used for cross-linking characterization so far. Among the LCI 
techniques, Guerrero et al., [175] published an approach which was able to reach a 
spatial resolution of about 17 um. It was based on a phase estimation of intensity 
extremes and shows a theoretical refractive index resolution of 1074 which did not 
take any thermal or noise influences into account. Shortcomings of the method are 
the restriction to measurements of the differential refractive index as well as its 
dependence on intensity measurements which are influenced by noise. 

Beside these technologies, some other imaging approaches might be interest- 
ing for cross-linking analysis. Singh et al. have shown that corneal cross-linking 
can be indirectly imaged by using an OCT-based approach, [176]. An evaluation 
of the mechanical stiffness was performed by analyzing the damping vibrational 


22 2 Related Works and Basic Considerations 


response of the cornea by an OCT system. Other works have shown the applica- 
bility of nuclear magnetic resonance (NMR) imaging for polymer characterization, 
[177]. The resolution of NMR is strongly dependent on the natural line width of 
the molecules that are measured, and the magnetic field gradient. The method is 
primarily suited to characterize samples in a laboratory environment. A different 
approach would be time-of-flight secondary ion mass spectroscopy (ToF-SIMS) 
where spatial resolution in the axial domain is achieved. By structuring the sample 
under test with an ion beam, the cross-linking depth profile is acquired in depths up 
to a few nanometers, [178]. 

None of these technologies is suited for the in-line characterization of polymers 
in production or for products. 

The main ambition in the characterization of polymeric cross-linking is to estab- 
lish a measurement technique which is capable to measure non-destructively on 
product level, is capable of a high refractive index resolution (better than 107°), to 
measure fast enough to be process-integrated and offers spatial resolution at the same 
time. None of the technologies known from literature combine these characteristics 
so far. 


2.3 Film Thickness Measurement 


Appropriate metrological tools for in-line characterization have to fulfill different 
requirements in distinction to classical lab-based technologies. In particular, tools 
are often required to measure a certain area (~ some um?) with high axial resolution 
in accordance with the speed of the respective processing step and perform measure- 
ments autonomously for a large variety of materials. As the amount of data gathered 
is usually large, appropriate algorithms have to monitor key values constantly and 
only report on deviations. This applies especially to the production of high volume, 
complex multi-material systems in roll-to-roll processes as in the production of 
thin-film systems, [179]. 

State-of-the-art technologies for thin-film characterization most often incorpo- 
rate variants of ellipsometry, [180] and reflectometry,[181]. While both methods are 
capable of being production accompanying tools [182], they usually cannot provide 
high lateral resolution, [183, 184]. As experimental works showed, higher resolu- 
tions and the ability to measure multi-layer systems, one of the most challenging 
problems in both technologies, is the knowledge of material constants and models 
in order to find or fit correct start parameters as well as converging criteria for fits, 
[185-187]. 
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Alternatively, optical and mechanical properties of thin-films have been char- 
acterized by interferometric approaches, [188]. While reaching high lateral resolu- 
tions, these technologies tend to be slow, as they rely on z-scanning of the OPD 
between the sample and reference arm, [189]. With the development of k-scanning 
approaches this disadvantage was overcome as mechanical scanning in the axial 
dimension was substituted by a tuneable light source to scan through k-space, [190, 
191]. However, these approaches are limited at both ends of the thickness measure- 
ment range. On the one end, the typical spectral scanning range of 20nm limits 
the maximum measurable thickness to about 15 um. On the other end, the mini- 
mal resolvable thickness is limited by the distinguishability of Fourier peaks which 
are unique for every reflection of a material interface. The thickness as well as the 
amount of dispersion of every material contributes to the width and position of 
each Fourier peak. Ghim et al. demonstrated that 500 nm is a typical minimal film 
thickness which can be resolved with this technique, [192]. 


2.3.1 Spectral Reflectometry 


A common method for the determination of film thickness and material parame- 
ters of thin-films in the semiconductor industry is spectrally-resolved reflectometry, 
[181]. This technique captures the spectrum of a broadband light source back- 
reflected from a sample’s surface and tries to fit it to a model spectrum that was 
calculated beforehand. During measurement, a reference spectrum, from e.g. a pure 
silicon surface, is captured first in order to reference absolute intensity values later. 
Afterwards, the reflection spectrum of a substrate with a thin-film is captured by 
illuminating the sample with a broad spectrum under a defined angle 9 (typical 
0 = 0°). The back-reflected spectrum is typically collected with a fiber and spec- 
trally decomposed using a grating spectrometer. For the determination of the film 
thickness, the real and imaginary part of the substrate’s refractive index as well as 
that of the film material one needs to know. By applying e.g. a brute-force fitting 
routine, the error between measured and calculated spectrum is minimized to obtain 
the film thickness. In order to calculate all possible reflection spectra for different 
thicknesses of a material system, the transfer-matrix approach is typically applied, 
[193]. Using this approach, the E-field is described as a matrix for every layer of 
material as well as for every transition between different materials. Each matrix M 
accounts for the reflectivity and transmission characteristics of every material and 
transition, [181]. 

In practice, a variety of theoretical spectra can be calculated with an assumed set 
of different film thicknesses in order to evaluate the measured film thickness. Usu- 
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ally, the root-mean square error (RMS) is determined between measured and calcu- 
lated values. Additionally, it is possible to perform calculations for other parameters 
like the refractive index (real and imaginary parts) in relation to the wavelength on 
the basis of knowledge of the film thickness. The technique is advantageous because 
of the low instrumental effort, the lack of sample preparation and the fast data acqui- 
sition times. However, the method is tied to the knowledge of the material model of 
a sample and caution is necessary when handling measurements from samples with 
changing material parameters which can occur in a production environment. Some 
research has shown, that the parallel acquisition of data under a range of polarization 
angles in combination with a camera might enable the usage of spectral reflectom- 
etry in process-accompanying for complex multi-layer samples, [194]. Currently, 
research has proven the possibility to combine technologies such as spectral reflec- 
tometry with other techniques such as low-coherence interferometry, spectroscopic 
ellipsometry and Raman spectroscopy, [195-199] 


2.3.2 Spectroscopic Ellipsometry 


A widely used method to determine the dispersion properties of thin-films is spec- 
troscopic ellipsometry (SE). This technique evaluates the change of polarization of 
light reflected from one or more layers of a thin-film of material. In particular, the 
technique is based on the determination of the complex quotient of the reflection 
coefficient of both polarization components. Typically, the amplitude quotient Y 
and the phase difference A of s- and p-polarized light are captured as a measure, 
[200]. The data of these measures can be gathered either at a discrete wavelength or 
at a broad spectral range using a spectrometer as detector. In order to determine the 
dispersion of a thin-film, the geometrical parameters such as film’s thickness and 
angle of incidence as well as the refractive index of the substrate material must be 
known. By using an appropriate mathematical model of the material system, the real 
and complex part of the films’ refractive index can be calculated. Under the assump- 
tion of a measurement in a broad spectral range, this approach can be used to derive 
the dispersion of the film material. As [180] proved, appropriate models can be used 
to determine the dispersion of various materials in multi-layer material systems as 
well. Park, [201] has shown that the determination of the real and imaginary part of 
the dielectric function, hence the refractive index, the extinction coefficient and the 
dispersion in the range of (190—1000) nm is possible. The measurements were per- 
formed in an angular range from 45° —75° in 5° intervals. In the model, the authors 
took into account the material properties of the surroundings, a rough surface, a thin- 
film of cadmium sulfide (CdS) as well as for a substrate material. The film thickness 
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of the rough surface and the CdS layer were determined by minimizing the error 
between calculated and measured values for the amplitude quotient W and the phase 
difference A. Additionally, the gathered data in relation to the measured spectral 
range could be used to determine other material characteristics such as characteris- 
tic peaks which are dependent on the crystal structure of the material. Furthermore, 
the determination of the dispersion (according to the Wemple-DiDomenico model 
[202]) could be used to evaluate the bandgap of the material. Other works, such as 
[203], demonstrated the determination of optical properties on a variety of thin-film 
materials under the usage of different dispersion models. Obviously, the knowledge 
of material properties and the application of the correct model are important in 
order to obtain precise results using SE. Some developments proved that an appli- 
cation of the technology as a process monitoring tool could be possible. Fried et 
al. [182] demonstrated that by a modification of the setup and the simplification of 
the fit function, significantly faster data acquisition and analysis are possible. For 
example, in order to evaluate copper indium gallium selenide (CIGS) solar cells, 
a correlation of the fit values with the amount of gallium (Ga) in the material was 
performed. This reduces the number of adjustable fit parameters to one and speeds 
up the calculation. The authors show that a modification of the fit model can be done 
for other materials which are relevant for solar cells as well. It was also shown that 
a modification of the mechanical setup from using a collimated beam to the usage 
of a so-called expanded beam can decrease the time needed for data acquisition. 
This setup allows for the illumination and data acquisition on a large area and over 
various angles of incidence at the same time. Also, the simultaneous acquisition of 
data under several wavelengths and angles was demonstrated. When the sample is 
translated under the measurement spot, the fast measurement of areas with several 
cm? becomes possible. The authors of [183, 184] showed maps of film thickness for 
various material systems on a wafer. Additionally it was shown that a line projec- 
tion of the expanded beam can be utilized to perform film thickness measurement 
as process monitoring in a roll-to-roll production environment. Although the results 
are promising, it has to be clarified that this approach only works in settings were 
the material is well known and the production focuses on a small number of material 
systems (in [182] the manufacturing of thin-film solar cells). 

According to [187], the typical usable spectral range of SE is (130-2000) nm 
whereas most practical setups use a reduced range. While the lateral resolution of 
production accompanying is rather low, experimental works have shown that a lateral 
resolution of 4um can be reached using a microscope objective, [186]. However, 
this makes scanning of the sample a necessity in order to measure reasonably large 
samples. 
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Modifications of spectroscopic ellipsometry have been combined with scatterom- 
etry in order to perform the control of optically critical dimensions in the semicon- 
ductor industry, [204]. The measurement of the smallest lithographic structures like 
the 7-nm-node, introduces new challenges to metrology as e.g. dielectric functions 
for layers thinner < 10nm. According to [187], future challenges for SE lie primar- 
ily in the increase of accuracy of the optical models as well as in the generation of 
new models. Furthermore, approaches for the combination with other techniques, 
[205], as well as for parallel acquisition of different parameters over a large spectral 
and angular range exist, [206]. In this area, it will be necessary to develop new and 
adapted components such as light sources with shorter emitting wavelengths,[207]. 

Some current research has tried to minimize the disadvantages in terms of the 
need for scanning by developing one-shot techniques using broadband light sources, 
a modulated carrier-frequency from an interferometer and a common spectrometer, 
[208]. Other works have tried to extend known analysis models by new mathemati- 
cal approaches or by the combination of data from SE and other technologies such 
as reflectometry, [209-211]. Some works have shown imaging SE techniques which 
suffer from poor lateral resolution (~ 60 um) and the need for temporal effort to 
perform the measurement and analyze the data (~ 8 s per wavelength in AA = (400 - 
700) nm), [212, 213]. Additionally, some work has been done to combine spectral 
ellipsometry with low-coherence or phase-shifting interferometry in order to get 
information on the surface profile of the sample alongside with the thin-film thick- 
ness, [214, 215]. While the resolution for film thickness and surface profile were 
in the nm-range, the combination of multiple technologies imposed new obstacles 
in terms of data fusion. Spectral ellipsometry is especially powerful in laboratory 
situations where the material model of a sample is well understood. In this case, 
thickness resolutions in the sub-nm can be achieved. In situations like an produc- 
tion environment, SE suffers from its low-lateral resolution and time-consuming 
data acquisition/analysis. 


2.4 Material dispersion 


The refractive index, Nabs, of a material is a measure of the refraction of light 
traveling through a material. Since it is defined as the relation of the vacuum light 
velocity co to the light velocity in the material c and can be described as the optical 
resistance of the material, [216] pp. 94, 


co 


Nabs = —- (2.5) 
È 
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As most practical applications are operated in air rather than in vacuum, the term 
relative refractive index, Nye], is more common. This relative index is defined as the 
refractive index in relation to the refractive index of air, ngj;, 


Nabs 


Nrel = G (2.6) 


Nair 


The dispersion of a material describes the relation of the refractive index, and there- 
fore the phase velocity of the light, with respect to the wavelength. This implies that 
the refractive index is wavelength dependent n = n(A). Normal and anomalous dis- 
persion can occur. In the case of glass or most polymers, normal dispersion is present 
in the visible wavelength range. In this case the material shows comparatively high 
refractive indices ranging from (300-500)nm and significantly lower refractive 
index values in higher spectral regions (e.g. (500-900) nm). Anomalous dispersion 
consequently shows an opposing behavior. In materials with relatively low optical 
density such as glasses or polymers, the dispersion characteristic is closely related 
to the absorption behavior. In general, the material behavior can be described using 
a quantum mechanical model. The application of a simplified description utilizing 
an electro-magnetic model is feasible for materials with low optical density, [193]. 
This model describes the influence of an electro-magnetic wave on the moleculary 
structure of a material. Each bond charge of a material’s molecular structure has a 
specific resonant frequency wj; which defines its absorption behavior. Dependent on 
the material properties, an incoming electro-magnetic wave excites the molecular 
structure which than leads to a specific refraction behavior. The resulting behavior 
of the real and imaginary part of the refractive index, n and k’, can be described 
mathematically using a resonator model of the following form dependent on the 
angular frequency w, [217], 


(in — ie)? = Nae” y i ; 2.7) 
Eom (o — o?) +iyjo 
j j Yj 


where Nn is the number of molecules per unit volume, e and m are the charge 
and mass of the electron, £ọ is the vacuum permittivity, F; is the strength of the 
absorption and y; is a measure for the frictional force at the resonance frequency. 
The KRAMERS- KRONIG model takes this relation into account, [218]. It describes 
the relation of absorption and dispersion of light in a material by combining the real 
and the imaginary components into one complex model, Fig. 2.1. 

This relation has a general validity for all materials. For spectral ranges with 
negligible absorption of the specific material (k’ — 0), the resonance frequency 
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Figure 2.1 Depiction of the slope of the real part n(w) and the imaginary part k’(@) of the 
refractive index for a simplified resonator model with a resonant frequency wj 


tends to differ significantly from the frequency of the incident light. For this case 
the relation can be simplified and expressed in terms of the wavelength A/A j 


N 2 F; N 2 ay 
n-1= > 2 — = En >, : 2? (2.8) 
Eom 7 w — w An*c-eom 7 à2— À 


where one can abbreviate the material specific terms with A ; 


F Nme? F ji T 
J An2c2egm’ = 
which leads to the simplified equation 
Aji? 
2 = J 
n> — Ir (2.10) 
j j 


When analyzing real materials like glasses or polymers, the absorption characteris- 
tics differ very much from the simplified model. Most notably, multiple absorption 
peaks exist in contrast to only one single absorption peak at one resonant frequency. 
Furthermore, these peaks inherit a distinct structure which leads to a complex dis- 
persion characteristic. For this reason, a number of approximation relations exist. 
These relations are always valid in a well defined spectral range which is most 
likely far away from a strong absorption band. An extrapolation outside the valid 
spectral range is not considered valid. Over the years a broad range of approxima- 
tion relations have emerged for different materials, spectral ranges and achievable 
accuracies, Tab. 2.2. 
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The measurement of the refractive index of a medium is usually defined in 
relation to the refractive index of air, Eq. (2.6). For this reason, the conditions 
of the surrounding air should be well defined in terms of temperature, pressure 
and humidity. Calculations are based on typical values for laboratory conditions, 
Tab. 2.3. 

In order to perform these measurements for refractive indices at defined single 
wavelengths with high accuracy, methods like refractometry (accuracy + 1076) 
and interferometry (accuracy + 10-7) are preferred. Measurements are typically 
performed at one or few defined wavelengths of a gas discharge lamp such as the 
sodium line at Ap =589.592nm. The influence of the temperature T, pressure p 
and humidity w on the refractive index of air and therefore on the relative refractive 
index measurements of materials can be described as follows 


dnair Nair — 1 


= = 40.268 x 107% hPa! (2.11) 
dp p 
dnai pai 
Pair La ia (2.12) 
dT l+aT 
dnair 41 x 10-°h Pa! 
nar a AR BES = 0010 K 10 hPa. (2.13) 
dw lt+aT 


where the thermal expansion coefficient is introduced as a. From these equations, 
it can be calculated that minor changes in the environmental conditions such as a 
temperature change of dT =+5K or a pressure change of dp=+ 20hPa can have 
a significant influence on nair. For comparability, the joint commission for spec- 
troscopy has developed an engineering equation for the determination of the refrac- 
tive index of dry air (with 0.03 % C O2 volume content) nair at T=15°C and the 
reference pressure po =760 Torr in the spectral range between (200-1350) nm, as 
well as for other temperatures of air under the same pressure in the visible spectral 
range, [217]. 


2.4.1 Thermo-optic coefficient 


Analogous to the behavior of air, the refractive index of a material is dependent on 
its temperature where the discrete form can be noted 


n (Ak, Tiga) —n (Ax, Ti) _ An (Ax, Ti,i+1) 


(2.14) 
Ti41 — Ti ATi i+1 
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Table 2.2 Models for the description of material dispersion according to [217] 


Properties Equations Remarks 
CAUCHY precision in (VIS) n(A)=at+ 3 + 7 first model from 

1074 1830 
HARTMANN precision in (VIS) n(à) = no + ao constants have to be 

10-4 — 1074 found empirically, 

05<B<2 

SELLMEIER precision nay = based on physical 

10-5 —5- 10° 1 N ajx model, very good fit 

+ oe > 
Yj (P25) characteristics 
HELMHOLTZ- based on SELLMEIER, | 2(A)? = see also [193] 
KETTELER- DRUDE | precision similar ao + XLi ER 
: i 


SCHOTT 


precision 
1075 — 1076 


n? (A) = Ao+A14?+ 
Ana? + Azı* + 
Agr + Ası 3 


derived from 
SELLMEIER equation; 
valid over large 
spectral range and 
for nearly all glasses 


HERZ- BERGER 


precision 1074 (VIS) 


n(A) = Ao + Ata? + 
A 


based on analytical 


A2 3 z aha 
aot + or investigation; 
oo OJ ag = 168nm 
GEFFCKEN precision 1076 n(à)=1+[1— combination of two 
D(A)|[n(1) — 1] + | functional relations 
B()[nQ2) — (spectral behavior of 
n(A1)]+ D(A)d(A) |a normal glass and 
abnormal dispersion 
from that reference) 
BUCHDAHL precision 1—2- 10-4 | n(w) = reference 
no + viw + vw? + | 49=587.6nm (VIS); 
v3@ +--+ vno” with| @=2.5=const. for 
a(A) = a) glasses; fast 


conversion for w of 
j=2orj=3 


For practical reasons this dependency is often evaluated by measuring the refractive 
index at different FRAUNHOFER lines A, at different temperatures of interest T;, 
[219-221]. An appropriate fit of this equation has to determine six parameters when 
calculating a single temperature. When incorporating the correct thermal depen- 
dency, that number of necessary fit parameters increases to eighteen. In order to 
simplify the calculation, some boundary conditions have to be implemented. As the 
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Table 2.3 Typical environmental conditions for the calculation of influences regarding the 
refractive index according to [217] 


parameter value 
temperature T 20°C 
thermal expansion coefficient a 0.00367 K~! 
wavelength 589.592 nm 
air pressure p 760 Torr 
partial pressure of vapor w 10 Torr 


influence of temperatures on the refractive index is relatively low, it can be described 
by using a simplified SELLMEIER model. Hoffmann et al., [222], have shown, that 
a reduction to a one-term model (i = 1) is sufficient where measurements are 
performed at only one wavelength X in reference to a base wavelength A 


dn(a,T) _ n?(A, To)— 1 


= 2.15 
dT 609° A 
Eo ROE -T 
(+ 2010 -m +3007 paee 2), 
Rei 
which results in 
2(A, To) — 1 
AnA, T m= 0) 
n(A, To) 
Eo(T — To) + E1 (T — To)? 
(Dur 104 Dur -m+ Dor ree > = 2), 
0 
(2.16) 


after integration. Usually, calculations are performed with a reference Temperature 
To and a single measurement temperature T in order to fit the values for the thermo- 
optic coefficients Do, D1, D2, Eo, E1. For commonly used glasses these coefficients 
can be found in tables, e.g. in [220]. 
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2.4.2 Photo-elastic influences 


A second important influence on the refractive index and therefore on the dispersion 
is (mechanical) stress acting on a material. While materials like glasses or polymers 
are usually showing isotropic behavior, the application of mechanical stress o can 
lead to anisotropic behavior of the refractive index. The reason for this is that the 
refractive index is dependent on the electric field vector in relation to the stress plane 
in a sample. For a calculation, the relation of the refractive index to the stress plane 
can be described in parallel, ny and perpendicular, nı orientation, [217], 


dn 
nj =n+—o=n+k, (2.17) 
do 
and 
dn, 
ny =n+——-o=n+K|. (2.18) 
do 


If the deformation as a result of stress is elastic, HOOKEs law can be applied. In 
this case dn/do can be written as the stress optical coefficients Ky and K1. These 
coefficients are known for most common materials and are to be determined for 
example in a 4-point-bending test. In materials like glasses, the refractive index 
typically changes equally for both directions in relation to the stress plane. In the 
case of hydrostatic pressure p, not only the refractive index of the material itself but 
also of its surrounding medium will change. In this case the change is determined 


as 
dn 


— =Kı+2Kı. (2.19) 

dp 
Acousto-optical modulators are applications that make use of this effect. The stress- 
optical constant K can be found for a variety of materials in tables such as in [223]. 
In some cases, it might be necessary to determine it experimentally for a specific 
material. An experimental setup would consist of asample under uni-axial load such 
as pressure or tension where a polarized beam of light illuminates the sample. The 
angle of the polarization will be 45° towards the main stress axis. The polarized 
beam of light can be denoted as a superposition of one component perpendicular 
to the main stress axis and one component parallel to this axis. By applying stress 
on the sample, the two components will experience an optical path difference As 
dependent on the sample length / which can be described with 


d d 
As = (n+ Zlo -n - ei = Ken, (2.20) 
do do 
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This can also be written as a phase difference A® between the two components 
with 


27 27 
A®= x (Ki Kı)ol= a (2.21) 
with K = K| — K]. (2.22) 


Clearly, a material’s stress optical constant is also dependent on the dispersion of the 
material, [224]. Furthermore, the stress optical coefficient K shows a temperature 
sensitivity. Hoffmann et al. [225] showed that usually only the stress-thermo-optical 
coefficients Ag, Aj, A2, B of the equation: 


Kà, AT) = Ap + A1- AT + Ap AT? + ER (2.23) 
AD 
with AT =T- Tp (2.24) 


are relevant for glasses. For most materials, the temperature tends to have a minor 
influence until the glass transition temperature is reached, [226, 227]. 


2.4.3 Characterization of dispersion 


The knowledge of the dispersion of materials and optical components plays a sig- 
nificant role in several areas of photonics. For example, the slope of dispersion in 
optical fibers in communications, due to different mechanisms like waveguide and 
material dispersion, determines the bandwidth and range of a transmission system, 
[228]. During the construction of laser sources with ultra short pulses the knowl- 
edge and control of the dispersion of components is crucial in order to ensure the 
generation of short pulses with high energy in a determined spectral range, [229]. 
Components such as fibers and mirrors have high demands regarding the determi- 
nation of their dispersion behavior, [230]. According to the specific requirements 
of each application, different metrology approaches are common. 


Time-of-flight measurement and phase-shift method 

In the field of optical communications, methods which utilize the direct measure- 
ment of the propagation time of an optical pulse as a measure for the dispersion 
characteristics are common and known as time-of-flight measurements. As com- 
munication lengths are typically very long (tens to hundreds of kilometers), dis- 
persion has a heavy impact on the propagation characteristics. Dispersion can lead 
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to time delays of pulses and to spectral broadening of the optical pulses, [231]. In 
order to increase the transmission speed and data rates, modern systems implement 
approaches such as wavelength division multiplexing (WDM) where information is 
sent as spectrally fine separated pulses. The occurrence of dispersion influences the 
ability to clearly separate the pulses and therefore the data sent. Additional compen- 
sation mechanism have to be implemented, [231]. As a method of characterization, 
the direct measurement of time delays between sent pulses is typically used in opti- 
cal communications to characterize the dispersion behavior. A major advantage of 
such methods is that it can be used also in already existing fiber transmission install- 
ments, [232]. In this way, not only the dispersion of the fiber, but also of the system 
as a whole can be determined. In the case that fiber sections are 4km or longer, a 
significant amount of dispersion is present, so conventional metrology can be used 
to determine the time delay between pulses. The measured delay is a superposition 
of different dispersion mechanisms such as material and waveguide dispersion. 

Typical light sources used in this context are RAMAN lasers, superluminescent 
diode (SLDs), erbium-doped fiber lasers and semiconductor lasers. According to 
[232], the main requirements for these light sources are a broad spectral range, a 
good tunability of a certain center wavelength and its full width half maximum as 
well as a high spectral power density. 

One major disadvantage of this method is the requirement of appropriate hard- 
ware in order to achieve the necessary temporal resolution. In the case of samples 
with relatively low dispersion (circa 100 fs/nm-km) it is necessary to have samples 
with a substantial length in order to characterize them with standard equipment. 
Fast photo detectors are able to work at rise times of a few picoseconds (circa 15 
ps according to [233]). Although [233] could implement an autocorrelator in order 
to realize the measurement of very small phase differences, the shortest measurable 
fiber length was not shorter than 100m. An application of the time-of-flight method 
to characterize samples with very low dispersion or very short lengths in the range 
of a few mm or um is not possible. 

A further development of the time-of-flight method is a phase-sensitive detection. 
In this approach, light sent through a sample is modulated with a sine signal which 
can be detected by an appropriately designed phase-selected amplifier circuit, [234]. 
This technique, in contrast, allows to delay measurements, the detection of much 
smaller differences in the dispersion related delays (with a resolution of about 10 ps). 
Problems arise from low light intensities and a bad SNR, [235]. Despite these 
disadvantages, this method is widely used as a resolution enhancement of time- 
of-flight measurements with rather long fiber lengths in the telecommunications 
industry. 
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Time-domain interferometric measurements 

A method to cope with the requirements regarding the temporal resolution in time- 
of-flight measurements is interferometry in the time domain. For this approach, an 
interferometer (e.g. of the MICHELSON type) with one fixed and one movable arm 
is set up, [236]. Usually broadband light sources with low coherence lengths, such 
as SLDs and gas-discharge lamps filtered using a monochromator, are used in these 
setups. In a typical experiment, the movable reference arm is translated in a con- 
tinuous fashion while the corresponding combined intensity signal of both arms is 
recorded with a photo diode. The shape of the recorded interferogram depends on the 
coherence properties of the light source. This method is similar to coherence scan- 
ning interferometry used for profilometry, see subsection 2.1.5. The execution of 
similar experiments using different center wavelengths is the basis for the dispersion 
measurements. After these experiments in a dispersion-free interferometer, a sample 
of interest can be introduced in one of the interferometer arms. After the recording 
of several new interferograms at interesting wavelengths, a temporal delay of the 
maxima of the interferograms in relation to the dispersion-free measurements can 
be done, [237-239]. A fit, using an appropriate dispersion model such as the Sell- 
meier model, can be performed using the maxima of the recorded interferograms. 
The temporal resolution in this approach is strongly dependent on the mechanical 
resolution of the translation stage. The step width and its repeatability define the 
possible measurable temporal differences between different interferograms. 

The main disadvantage of this approach is the effort necessary to perform the 
measurements. Although broadband light sources can be used, the dispersion char- 
acterization over a large spectral range with a high spectral resolution takes time. 
Also, the appropriate filtering for each measurement has to be carried out. Fur- 
thermore, the accuracy of common translation stages limits the temporal resolution 
of the approach. A characterization of small dispersion, i.e. of thin structures (um 
to mm) is not possible due to instrumental constraints. Therefore, the approach is 
mostly used for the characterization of optical fibers with at least a few centimeters 
of length, [240-242]. 


Frequency-domain interferometric measurements 

Analogous to time-domain measurements, interferometry in the frequency domain 
uses broadband light sources as well. In contrast, these sources are not spectrally 
constricted. Furthermore, the interferometer is usually built using two fixed arms. 
The detection of the signal is performed on the complete spectral bandwidth of 
the light sources using a spectrometer, [243]. A typical signal of the recombined 
intensity of both arms leads to a modulated signal over a certain spectral bandwidth 
with a characteristic frequency. The frequency is dependent on the difference in 
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the arm lengths of the interferometer as well as on the dispersion of the system, 
[244]. The insertion of an dispersive element leads to a characteristic change in 
the modulation and its frequency. The determination of dispersion can be done by 
different methods in this approach. On the one hand, the modulation can be analyzed 
using a Fast Fourier-transform (FFT). On the other hand, the analysis can be based 
on the determination of the so-called stationary phase point. 

The usage of an FFT-based approach enables the fast determination of very 
small dispersion which can occur due to material characteristics or small amounts 
of media. The analysis is similar to FD-OCT, although the broadening of the Fourier- 
peaks due to dispersion can be neglected if dispersion is very low. Otherwise, 
complex dispersion compensation methods have to be implemented. Liebermann 
et al. [245] showed an experiment to determine the dispersion characteristics of 
distilled water using a MICHELSON interferometer in a free-space configuration, 
utilizing a supercontinuum light source as well as two combined spectrometers 
—(350-1100)nm and (900-1780) nm—to cover a large measurement range. For 
the analysis, two FFT-based approaches were compared. The so-called indirect 
approach transforms the signal initially from the frequency into the time-domain 
using zero-padding in order to achieve an even distribution of the signal after back 
transformation. The transformed signal is filtered afterwards and Fourier trans- 
formed for a second time in order to extract the spectral phase. A Taylor expansion 
enables the calculation of the dispersion coefficients. The second analysis method 
discussed in [245] is the so-called direct approach. In this approach, the spectral 
data is Fourier transformed into the non-even distributed time-domain and filtered 
using a Heaviside-Filter. After the back transformation and a Taylor expansion the 
dispersion can be determined. According to the authors, the indirect approach is 
computationally more time-consuming but also more suitable to characterize sam- 
ples with large dispersion. Opposing to this, the direct approach is relatively fast 
and suitable for samples with small dispersion in process-accompanying problems. 

Another approach to characterize refractive indices, is to fit dispersive interferom- 
eter data from experiments, [246]. The authors assumed that the frequency chirp of 
the spectral intensity was a measure for the dispersion. By fitting this data the fringe 
periodicity was analyzed for different sample configurations. Under the utilization 
of a mathematical model (e.g. Sellmeier) the wavelength-dependent refractive index 
was estimated. In order to increase the fitting quality, more complex fit methods or 
correction algorithms could be used, [247, 248]. 
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Surface Profilometry 3 


3.1 Experimental Setup 


As outlined in the previous chapter, existing technologies for surface profilometry 
show certain drawbacks in terms of resolution, dynamic measurement range, three- 
dimensional measurement capabilities and speed. The following chapter introduces 
a novel approach to surface profilometry which aims to provide solutions to the 
problems named. The basic setup for all experiments is centered around a two- 
beam interferometer of the MICHELSON type, Fig. 3.1. 

A WLS—white-light source emits a collimated beam which is typically split in 
a 50:50 ratio by a BS—cube beamsplitter. The reference arm typically consists of a 
DE—dispersive element with the thickness tp g. Following the transmission through 
the DE, the light is reflected off of a REF reference mirror which causes the light 
to transmit through the DE a second time before it is guided by the beamsplitter. In 
the sample arm, the light is reflected on a sample surface before it is recombined 
with the reference arm light by the beamsplitter. The optical path difference of both 
arms 6 is fixed as any change in ö is usually measured as a change in the samples 
surface profile! . The recombined signal is detected by a grating spectrometer. As the 
dispersive element plays a significant role in gathering relevant information on the 
topography of a sample, the term dispersion-encoded low coherence interferometry 
(DE-LCD is used throughout the work for this approach. 


! In practical implementations, one of the arms can be equipped with a translation stage 
for adjustment purposes. 
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Figure 3.1 Principle of a profilometry setup with WLS—white light source, BS— 
beamsplitter, DE—dispersive element (with thickness tpg), REF—reference mirror (trans- 
latable in one dimension for adjustment purposes), SMP—sample with a surface profile and 
SPEC—spectrometer 


3.2 Measurement Range and Resolution 


In general, the mathematical model of a two-beam interferometer applies in the case 
of the setup described in Fig. 3.1 according to [249] 


Epu = Ep- ei(otteo) — Ej- EER, (8.1) 


where the resulting electric field Eour is composed of the initial electric field Eo in 
combination with the oscillating portion defined by the angular optical frequency 
œ or the optical frequency f, the time f and the phase pọ which consists of the 
optical path difference ô, the refractive index of the surrounding medium n and the 
wavelength X. In an experimental setup, where only the time averaged signal of the 
electric field can be detected the intensity / is of interest, it can be formulated 


Cc: & 
Eou =, (Ei? + [Bal +2- E1 - E2) 3.2) 
Cc: £0 —?2 
l= | Futi (3.3) 


which is composed of the velocity of light c, the vacuum permittivity constant £0 
as well as the two electric field components of the interferometer arms Eı and E2 
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and is integrated over a given time. As the relative change between the two arms 
is of interest when measuring surface heights, the phase term holds the relevant 
information. Therefore, the intensity / (å) can be written as the spectral dependency 
of the phase with 


IQ) = Io- (1 + y (à) - cos go) (3.4) 
2m -d- Nair 


p = — (3.5) 


where Jp is the initial spectral intensity and y (A) is the spectral contrast of the inter- 
ference fringes. It is dependent of the OPD between the both arms denoted with ô, 
the wavelength X and the refractive index of air nair which was assumed to equal 
one for simplification. These equations are suitable under the consideration that the 
interferometer is dispersion-free. In this general case, the phase is changing propor- 
tionally with the wavelength. Through the introduction of a dispersive medium in 
one arm of the interferometer, Fig.3.1, the phase term g extends to [246, 249] 


[nPE A) — 1] tpe — ô 


=2 
p T 7 


(3.6) 


Assuming that pure material dispersion is the only effective mechanism?, the trans- 
formation of the interference signal is dependent on the wavelength-dependent 
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Figure 3.2 Simulation of a typical spectral interference signal / (A) from a single point having 
the OPD ô with a equalization wavelength A, =600 nm plotted in black and the corresponding 
phase signal (à) plotted in red 


? In the cases considered within this work, only material dispersion was used as a mechanism 
to encode surface information as the setup was implemented in free-space. By using a fiber- 
based setup, other mechanisms of dispersion such as waveguide dispersion might be utilized. 
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refractive index n? Ë (X) and the material’s thickness tpg. The periodicity of fringes 
tends to a minimum at the so-called equalization wavelength Aeq, which is depen- 
dent on ô, Fig.3.2. Using the interferometer as a profilometer, every height change 
in the sample’s surface changes the OPD and thus leads to a different equalization 
wavelength. The equalization wavelength can be calculated analytically from the 
intensity signal as the minimum of the derivative of the phase signal with respect to 
the wavelength 


1— nPE (heq) | toe +6 
G | P Oe] (3.7) 
A), a? 
with nPE (a) = nPE A) A. 2 (3.8) 
8 dà’ i 


where np (A) is the group refractive index of the dispersive element. It can be 
concluded from this equation that the dispersive element and the path difference 
ô have the most significant influence on the signal. In this setting, the DE defines 
the phase slope of the signal and furthermore, the axial measurement range where 
the equalization wavelength can be observed, Fig. 3.3 a). With increasing thickness 
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Figure 3.3 Simulation of the basic dispersion-related dependencies of the interferometric 
signal with a) dispersion-dependent behavior where (I) shows the signal with a DE of N- 
BK7 having a tpg = mm, (ID tpg =2 mm and (II) tpg =4mm as well as b) the dependency 
regarding the OPD for a N-BK7 DE with tpg =4 mm where (I) eq =0.5 um, (ID Acg = 0.6 um, 
(ID) Aeg =0.7 ym and (IV) representing the slope of 5 which is a function of the refractive 
index of the dispersive element n? 
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of the DE of the same n?*, the phase gradient around the equalization wavelength 
increases which results in a higher frequency intensity signal around A¢,. In practical 
terms, this signal modification is responsible for the effective axial measurement 
range. This measurement range Az(A) can be estimated within a spectral range of 
interest 
_ [DE 
Azü) = [n j= 1] ee 3.9) 


This characteristic of the DE defines the response of Aeq due to a change of the 
optical path difference, Fig.3.3 b) which therefore can be used as a measure for 
surface profile changes. During the design of a DE-LCI setup, the detector size and 
its resolution as well as its dynamic range in the desired spectral range set the initial 
boundaries for the axial measurement range and resolution. The characteristics of 
the DE enable fine tuning of the axial measurement range and resolution even after 
the initial design of the system. 

In a system like this, the axial measurement range and resolution are determined 
by the spectral probing range AA, the center wavelength Àc, the dispersive element 
and the detecting spectrometer configuration. In order to evaluate the possible axial 
measurement range as well as the corresponding resolution, the interconnection of 
the detector properties with the material properties of the DE have been studied. 
For this purpose, simulations for a system with a defined set of parameters have 
been performed, Tab. 3.1. As stated before, the center wavelength in relation to the 
refractive index of the DE and the spectral range have a characteristic influence on 
the axial height measurement range which follows an e.g. SELLMEIER-like slope. 
This behavior can be used to manipulate the measurement range in a pre-designed 
setup. According to Eq. (3.9) the measurement range Az(A) = Az can be calculated 
over a given spectral range Ad as a function of the group refractive index and the 
thickness of the DE, Fig. 3.4 a). The behavior, which is in most cases SELLMEIER-like, 
increases the measurement range for a DE, following this characteristic especially 


Table 3.1 Spectrometric system properties for initial simulations regarding measurement 
range and resolution of the designed DE-LCI 


component value 

number of spectrometer pixels n p 3648 

spectral range AA 500 nm 
center wavelength àc 610 nm 
spectral resolution Arspec 0.3 nm 
dispersive element (N-BK7) tpg 1,2 and 4mm 
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towards shorter wavelengths. With the aid of different dispersive characteristics 
from e.g. glasses, polymers or thin films, this behavior can be controlled in a wide 
range while the setup is kept constant, Fig.3.4b). 
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Figure 3.4 a) Simulation of measurable height ranges in dependence of the spectral range of 
the setup Az for different thicknesses tpg of N-BK7 and b) simulated measurement ranges 
for different materials with a tpg = 1 mm 


Light sources such as SLDs or swept-source lasers are commonly used in other 
LCI approaches. With regard to the system design of a DE-LCI, these light sources 
are considered suboptimal as the spectral bandwidth is usually about (80—120) nm 
with center wavelengths in the range of (800 - 1300) nm, [117, 118, 120]. The results 
of the simulation, Fig. 3.4 a), show that such a spectral range with a center wavelength 
at 1050nm would significantly limit the height measurement range. In fact, DE- 
LCI benefits from broadband light sources such as supercontinuum or laser-driven 
plasma light sources which provide a spectral power output from (300-2000) nm, 
[250]. 

With regard to the DE, a spectrally defined resolution parameter, rpg, can be 


calculated 
Az 


Ax 
Clearly, the resolution can change due to dispersion as a function of the center 
wavelength while the nominal spectral range is constant. Nonetheless, the center 
wavelength influences Az, see also Eq. (3.9). As the detection of spectral interference 
is usually performed in a spectrometer, its specific resolution properties have to be 
considered as well. Therefore, the resolution of the system, r,,;, can be calculated 
with 


FDE = (3.10) 
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Tsys =TDE Arspec- (3.11) 


Due to the construction of e.g. a grating spectrometer, its resolution r,pec is defined 
by the combination of the slit size, the grating constant, the detector size and its 
specific pixel size. In consequence, a set of typical values for the measurement range 
and axial resolution was simulated for three different dispersive elements, Tab. 3.2. 
From this simulation, it can be deduced that the axial height resolution is tied to the 
measurement range and scales linearly with the thickness of the dispersive element. 
Additionally, it becomes clear that the dispersive element holds a large potential 
to tune both measurement range and resolution in a setup where the remaining 
components are already selected and fixed as shown on Fig. 3.4. 


Table 3.2 Calculated system properties regarding measurement range and resolution of the 
designed dispersion-enhanced low-coherence interferometer for investigations in the spectral 
range of AA =360-860nm 


tpe [mm] Az [um] Tsys [nm] 
1 80.45 48 
2 160.9 97 
4 321.7 193 


In order to evaluate the simulations, a temporally controlled experiment was car- 
ried out with the setup seen in Fig.3.1. In contrast to classical temporal LCI, the 
experiment was designed to be controlled temporally but detected spectrally. For 
this purpose, a spectrometer according to the specifications of Tab.3.1 (AvaSpec- 
ULS3648 VB, Avantes BV, Apeldoorn, The Netherlands) was used to capture inter- 
ference data. The plain mirror in the reference arm acted as a sample which was 
translated to different OPDs in order to emulate height changes of a sample. The 
sample arm, which is also equipped with a plain mirror, was kept at a constant posi- 
tion. Aluminum coated mirrors with a flatness of 4/20 (EO Partno. 34360; Edmund 
Optics Inc., Barrington, NJ, USA), having a scratch-dig of 20-10, served as sample 
and reference reflectors. In the experiment, a dispersive element of tp g =2 mm (EO 
Partno. 49121; Edmund Optics Inc., Barrington, NJ, USA) was used. The height 
emulation was performed with a piezo-driven precision stage (SLC 2412, SmarAct 
GmbH, Germany). As a common reference, the setup was adjusted to an equaliza- 
tion wavelength of Aeg =500nm. Subsequently, the translation stage was used to 
move the reference mirror in steps of Aö=0.5, 1 and 2 um along a maximum mea- 
surement range of Az=60 um. The spectral interference at the spectrometer was 
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recorded and the equalization wavelength was determined for every emulated height 
step. Based on Aeq, the analysis of the measured height step, Zmeas, was performed 
relative to the previous position with Eq. (3.9), Fig.3.5. The results show that the 
emulated height steps could be measured with some statistical deviations. Specif- 
ically, heights of Zmeas =(0.497 + 0.098), (0.998 + 0.106) and (1.997 £0.108) um 
could be measured for the nominal emulated steps of Aö=0.5, 1 and 21m respec- 
tively. The detected deviations were only slightly larger than the calculated reso- 
lution of 0.097 um for the dispersive element of tpg =2 mm, Tab. 3.2. In reference 
to the work of Ruiz et al. [118], the dynamic range (DR) can be calculated as the 
quotient of the measurement range Az and the resolution r,y;. The simulations as 
well as the experiments show that the DR of about 1667 is constant for the different 
measurement ranges since it is limited by the hardware. Due to the limitations in 
available detector sizes and trade offs between measurement range and resolution, 
other analysis schemes have to be considered in order to enable a higher dynamic 
range. 
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Figure 3.5 Results for the measured emulated height steps zmeas for the targeted steps of 
Aö=0.5, 1 and 2 um respectively over a given travel range of a piezo stage in the reference 
arm 


It can be noted that the signal not only shows a dependency of ö towards a 
change in equalization wavelength but also towards the signal amplitude at Aeq. Ina 
more advanced analyzing scheme this behavior is the primary component of a two- 
part process in order to determine height profiles with high-dynamic range. When 
analyzing the intensity behavior of the signal at the equalization wavelength in a sig- 
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nificantly smaller ö-range such as 10: r,,; * 1 wm, an oscillatory behavior becomes 
visible, Fig.3.6a). This simulation reveals that the differences in the height depen- 
dent intensity behavior are negligibly small for the different dispersive elements. 
By deceasing the ö-range of interest even further towards the maximum hardware 
resolution of tpg =2mm to Aö=97nm, the intensity dependency becomes nearly 
linear and unambiguous, Fig. 3.6 b). This effect can be exploited in a more advanced 
analysis scheme where the first step consists of the determination of the equalization 
wavelength as a rough measure for height changes. During a second step, a fit of the 
intensity amplitude around this point can contribute to a fine-scaled analysis with 
resolutions beyond the hardware limit. Furthermore, the behavior can be considered 
independent of the thickness of the DE and hence of the measurement range. The 
simulation revealed that the intensity difference between the signals of tpg = 1mm 
and tpg =4 mm are about 0.012 % at maximum. Once the analysis becomes not only 
dependent on the wavelength accuracy but also on the intensity signal measured, the 
noise of this signal determines the height resolution. In the setup presented within 
this work, Fig.3.1, the actual measured spectrum, S;, is affected by four main noise 
components, of the light source, Pz s, the cameras chip, Peam, the cameras amplifier 
circuit, Pamp, as well as of the A/D converter, P4/p, which contribute to the mea- 
sured signal, So, Fig. 3.7. The impact of the light source as well as the camera sensor 
together with its amplification are highly influenced by the experimental conditions 
such as integration time and gain of the camera. On a more detailed level, the four 
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Figure 3.6 Simulated intensity change in relation to a height change A6 for different thick- 
nesses of the DE tpg of N-BK7 with a) oscillations visible over a Aö-range of 1 um and b) 
reduced region of interest for Aö where the slope is nearly linear and the maximum difference 
between different DEs is about 0.012 % 
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noise components represent more fundamental noise sources. The photon flux of the 
light source is the main source of statistical variation which leads to photon noise 
of this component. While receiving this fluctuating source of energy, the camera 
sensor converts incoming photons into electron-hole pairs with a given quantum 
efficiency, which in itself is a statistical process known as photo-electron noise. The 
camera is also source of photo-current noise, which arises when electron-hole pairs 
are converted to pulses of electric current. In addition to the before mentioned noise 
sources, the process is dependent on the area of the detector, its quantum efficiency 
and the integration time, [251]. If the incoming photons are Poisson distributed, this 
noise source is known as shot noise. Additionally, the amplification of the photon- 
induced electrical current will be statistically dependent as well as the quantization 
of the detector current in A/D conversion. Furthermore, thermal variations will cause 
deviations and random signal contributions in the electronic circuitry of the sensor, 
the amplification and the A/D conversion which are also known as receiver circuit 
noise. 
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Figure 3.7 Noise-affected components of the setup with the initial spectrum S; which 
is manipulated by the noise sources Pzs, Pcam; Pamp, Payp resulting in the measured 
spectrum So 


It should also be noted that thermal influences will affect the setup as a whole 
and lead to geometrical deformations of critical components. Due to the slow rate of 
change of these fluctuations and the expected short integration times of the detector, 
itcan be assumed that they have a neglectable influence on the signals. In determining 
a combined noise level of a typical measurement situation, these influences will be 
captured during a noise characterization measurement. 

In order to qualify the combined influence, an experiment was conducted to char- 
acterize the induced noise. For this purpose, the spectral intensity response of the 
system with a mirror in the reference arm (Thorlabs PF10-03-P01) and a typical 
sample surface (silicon height standard Simetrics VS, Simetrics GmbH, Germany) 
was recorded. The recording was performed at gain levels ranging from 0 to 12dB 
in steps of 3dB and with five equally spaced integration times per gain level. The 
integration time was used to control the relative magnitude of the intensity signal 
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in steps of 20% starting from 100 %°. The recorded spectral intensity was spline 
interpolated along the spectral range, which was recorded in order to gather the 
moving average of the data, Fig.3.8 a). Subsequently, the spline interpolated data 
was subtracted from the raw, gray-valued data. A normalized signal was computed 
to get information on the relative spectral noise, AJ, Fig.3.8b). For one particular 
spatial position, a constant fluctuation over the complete spectral range is visible. 
The analysis of the distribution of this data reveals that it can be modeled using 
a Gaussian function which fits the data set with a coefficient of determination of 
R? =0.979, Fig. 3.9 a). By analyzing the o of the Gaussian distribution, the averaged 
intensity noise was found to be AJ =+0.75 % (21.25 dB). To gain a better insight 
in the spatial dependency of the noise, the distribution of noise for a number of 
points along the spatial axis was evaluated. The coefficient of determination R? 
for a Gaussian function describing the noise at the individual point was calculated, 
Fig. 3.9 b). It can be seen that no spatial dependency is present. For this reason, the 
spatial domain was not investigated in further detail and all data was integrated over 
along this domain. 
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Figure 3.8 Analysis of the noise of the setup with a) raw measured spectral intensity in 
gray-values (abbrev. gw) and spline interpolation as well as b) relative spectral noise A/ with 
respect to the spectral range 


3 In this context, the 100% reference was defined as the integration time where the recorded 
intensity maximum was just not provoking overexposure. While this section only discusses 
the behavior of the camera for a gain of OdB, the appendix in the Electronic Supplementary 
Material (ESM) presents an analysis for the other gain levels. 
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Figure 3.9 Analysis of the statistical behavior of the captured intensity noise with a) distri- 
bution of the measured, relative intensity noise of all positions in the x-dimension for a gain 
of OdB and a relative signal magnitude of 100% with a Gaussian fit of the same having a 
mean R? =0.979 where AJ =0.75 % (21.25 dB) could be measured as the averaged intensity 
noise and b) spatially-dependent plot of coefficient of determination R? for Gaussian fits of 
the relative noise distribution 


In order to estimate the impact of this noise on the determination of the height of 
a surface, an analysis of the initial interferometric equation Eq. (3.4) was performed. 
The equation was solved for the path length difference ô, 


ge (nF) = 1) wee = -cos7! EF _ | (3.12) 


Where the derivative of the path length difference ô with respect to the intensity 
variation dI is given by 


dê dpf peg à d aif 2 
7-76 0) =1) toe] = 3:5 [oo (-1)| (3.13) 


and the path length uncertainty Aé is 


à 1 
Ae ——— AM. (3.14) 
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This relation describes the limits for the described profilometry approach, which is 
also dependent on the spectral range as well as on the relative intensity, Fig.3.10 
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a). As Eq. (3.14) describes, the resolution depends linearly on the wavelength. It 
could be determined that it has a slope of 0.0012nm per nm at Jy =0.5 arb. units. 
Furthermore, the resolution strongly depends on the probing intensity in a more 
complicated relation, Fig. 3.10b). In order to qualify this dependency, the intensity 
range in which the resolution increase is < 10 % of the minimal value, was calculated. 
It was found that this equals an intensity range of J =0.3—0.7 arb. units. For practical 
purposes, this means that the intensity range should be chosen together with an 
equalization wavelength as low as possible in order to achieve high resolutions. 
This can be ensured by adjusting the path length difference of the interferometer. 
In all experiments, the setup was adjusted to Aeq =562 nm. As the spectrometers 
detection range started at 447 nm, the chosen equalization wavelength enabled the 
acquisition of enough data for fitting in proximity of Aeq. This helped to minimize 
imaging distortions which are typically present close to the border of the detector. 
Using the values for Jp and A., as well as the intensity noise AJ in Eq. (3.14), a typi- 
cal resolution of Ad (Aeg ) =0.67 + 0.05 nm was calculated. The calculated detection 
limit according to Eq. (3.14) is valid for the analysis at one spectral position. The 
data analysis of the presented experiments utilized the fit of spectra in a region- 
of-interest (ROI) around Aeq which is explained in more detail in subsection 3.3.3. 
Specifically, the measured intensity data was fitted with Eq. (3.4) and (3.6), where 
the thickness of the dispersive element tpg and its refractive index nP= (A) were 
assumed to be known and ö(x, y) as well as y (à) were approximated. In order to 
account for utilization of the ROI in fitting, the single point detection limit Ad (Aeq) 
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Figure 3.10 Estimation of the resolution based on the measured intensity noise with a) three- 
dimensional dependency related to the measured spectral and intensity range according to Eq. 
(3.14) as well as b) detailed plot at a equalization wavelength of Aeq = 562 nm 
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was used to calculate the resolution of the fitted data, r fir, with the aid of an RMS 


approach as 
AS (Aeq)? 
ee uL (3.15) 
nf 


where n f is the number of spectral data points used for fitting. Within this work, 
n f =530 spectral data points were used so that the theoretical resolution was esti- 
mated as r fis =0.029 nm. With a measurement range of Az=79.91 um and the above 
calculated minimal resolution, the dynamic range is 2.75 x 10°. This value is sig- 
nificantly higher than the initial estimation which was solely based on the evaluation 
of Ae, which was DR = 1667. 

As this value is based on the model expressed through Eq. (3.14), the exper- 
imentally achievable DR might additionally be limited by other influences such 
as thermal fluctuations or the data processing routines which are not included in 
the model. Utilizing this analyzing scheme, the limitation on the axial resolution 
imposed by the thickness of the dispersive element was minimized. 


3.3 Signal Formation and Analysis 


As shown in the previous examinations, signal analysis in dispersion-encoded low- 
coherence profilometry needs to be based on a combined evaluation of the equal- 
ization wavelength Aeq and the amplitude at this wavelength in order to achieve 
sufficient resolution. This holds true especially for comparatively large measure- 
ment ranges. Different approaches to analyze the information in recorded spectra 
have been developed and assessed within this work. 


3.3.1 Fitting of Oscillating Data 


Conventional fitting approaches such as the Levenberg-Marquardt algorithm, con- 
verge typically fast for periodical signals with constant phases. In the case where 
the phase of the data is varying over one dimension as in DE-LCI, fits can converge 
fast as well but most likely on a local minimum rather than on the global minimum, 
Fig.3.11 a). In case of this simulation the actual OPD was 45; =535.233 um. A 
brute force calculation of spectra according to the model described with Eq. (3.4) 
and (3.6) and the respective squared sum of errors osse was calculated in a range 
of ö,ng =(534—537) rm. It can be seen that the osse value oscillates and shows 
two distinct local minima at 534.3 and 536.1 um apart from the global minimum. 
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An approximation with the Levenberg-Marquardt method using an initial guess of 
Öguess =536 ym converged in < 10 iterations to a value of ô fit =536.118 um which 
falls on one of the local minima. It is clearly visible that the resulting spectrum is 
significantly different from the one which was actually present, Fig.3.11b). One 
way to circumvent this problem is to compute spectra within a very large range of 
possible values for the fit parameters. As this approach can be time and memory 
consuming, one has to consider strategies like downsampling of the measured data 
or the introduction of more advanced fitting approaches, [252-254]. Here, Monte- 
Carlo-based methods can be used in order to guess different start parameters and 
perform fitting in the parameter range with the highest likelihood of convergence on 
the global minimum, [255]. Furthermore, another strategy is approximating spectral 
data in multiple stages with coarse variations of the fit parameters to determine the 
area of interest and finer variations to exactly converge on the final parameters. In 
either way, it is desirable to determine the start value ög4ess with high precision. 
It should be emphasized that under the consideration of a precise determination of 
good start values any established fitting approach such as Levenberg-Marquardt can 
be used instead of the brute-force approach. 
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Figure 3.11 Evaluation of classical Levenberg-Marquardt-based fitting approach versus 
brute force calculation for dispersion-encoded interferometric data with a) error progression 
of spectra using different path differences 6 showing several local and a global minima where 
fitting is likely to fall on a local minimum if initial fit parameters are not carefully chosen and 
b) resulting spectrum of the fit routine which converged on a value of ô fis =536.118 um for 
an initial guess of öguess =536 um in relation to the simulated spectrum 
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3.3.2 Frequency Analysis 


In order to evaluate ôguess from measured data, the equalization wavelength Aeq can 
be used in conjunction with Eq. (3.9) which follows from ae = 0, Eq. (3.7) 


Öguess = [nF Gea) = 1] ‘IDE; (3.16) 


where the value for the group refractive index of the dispersive element is called 
np E (eg) as well as its thickness tpg are used to calculate the initial estimate for the 
path length difference 59,255. However, this approach relies on the correct estimation 
of Aeg from measured data. Although this method describes a theoretical way to 
gather the necessary data, it is problematic in its implementation for measured data. 
In order to analyze the phase signal, a cos~!-operation can be performed which 
results in wrapped data ranging in [—z, 7x]. To avoid phase wrapping, the data 
processing established within this work relied on a frequency-based analysis. The 
analysis of periodic signals is typically performed by frequency analysis approaches 
such as Fourier transforms like FFT. This is the standard method in e.g. FD-OCT, 
where interference at different axial positions can be separated as different frequency 
components from the spectral domain signal, [105]. As a number of publications 
demonstrated, dispersion in the experimental setup or in samples, which might be 
composed of different materials, leads to measurement uncertainties. For this reason, 
a variety of dispersion compensation approaches have been researched, [256-258]. 

In distinction to approaches performing dispersion compensation of the gathered 
spectral data, DE-LCI aims to make use of the dispersion within the system. In order 
to gain information on the phase and its wavelength-dependent slope, a Fourier- 
based analysis was performed on small slices Sj, of the signal. Known as windowed 
or short-time Fourier-transform, this method assumes that the phase is constant in 
a small slice of the signal, [259]. It composes a resulting spectrogram by sliding 
an observation window in steps over the signal having a fixed length and overlap, 
Fig.3.12 a). By changing the window shape, window width Aw as well as the 
overlap between subsequent windows Ao, the resolution with regard to the phase 
minimum can be controlled. In consequence, a spectrogram of the stacked and 
Fourier-transformed slices §(S/,) is used to analyze the phase minimum and the 
resulting equalization wavelength Aeq, Fig. 3.12 b). The approach can also be used 
in analyzing signals composed of multiple frequencies from different axial positions. 
In the context of DE-LCI, the approximation of the equalization wavelength Aeq was 
used to calculate an initial value of öguess in the interferometer in order to enable 
the fit to converge fast on a global minimum. 
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Figure 3.12 Principle depiction of the STFT where a window of the width Aw is continu- 
ously slid over the signal in the direction ds to form spectral slices $/, with an overlap of Ao 
and b) resulting spectrogram of the stacked, Fourier-transformed slices §(S;, ) 


3.3.3 Two-Stage Fitting 


As described in subsection 3.3.1, the usage of conventional fitting routines on oscilla- 
tory data with varying phase can lead to problems. For this reason, the here developed 
fitting routine was constructed as a two-step process. Using a range of Ad) =+ 1 ym 
with a step size of 2nm centered around the previously calculated öguess, a set of 
simulated spectra based on Eq. (3.4) and (3.6) was calculated in a brute-force fash- 
ion. The determination of the error sum of squares (SSE) of these calculated spectra 
with respect to the measured spectrum enabled the estimation of a more precise 
value for the path length difference ô; at the minimum of the SSE curve, Fig. 3.13. 
The calculated 5; was used in a second iteration of the routine to calculate another 
set of spectra with a finer spacing in Ad) = + 140 nm with steps of 0.02 nm*. Compa- 
rably, the SSE of the calculated spectra was evaluated with respect to the measured 
spectrum. The minimum SSE indicates the path length difference 62 which can be 
used to compute the height at a point of the sample, see also Fig. 3.14b). 

The described method was chosen instead of other established fitting algorithms 
to ensure the convergence on the global minimum rather than a local minimum 
which can be the case due to the oscillating nature of the data. The iterative fitting 
approach bears further potential for optimizations regarding the processing time. 


4 These values were typically used throughout this work. They are dependent on the used DE 
as well as on the sample to be measured. 
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Figure 3.13 Visualization of the simplified two-step fitting process based on SSE deter- 
mination for simulated data sets in two ranges (Ad; and Aô2) for the path length difference 5 


Opposing to doing a fit of the whole spectrum, fitting in a ROI was carried out 
in order to reduce processing time, as only about 25 % of the gathered data had 
to be processed. The processing time of a whole profile with this spectral ROI 
was about 2 seconds’. The ROI was selected as a fixed set of 530 data points 
distributed symmetrically around the equalization wavelength. The size of the ROI 
was determined in preliminary experiments in order to include at least one spectral 
modulation to each side of Aeg which is dependent on the used DE. It is not expected 
that the size of the ROI influences the resolution of the setup. 

As described above, the fit of the measured data in close proximity of the equal- 
ization wavelength significantly enhances the resolution, see Eq. (3.14). In order to 
prove the resolution limit experimentally, a sample data set was evaluated regarding 
its regression error towards simulated data sets within a range of path length dif- 
ferences Ad, Fig.3.14a). Due to the ability to perform a search for the regression 
minimum based on a SSE approach, the best fitting data set can be used to deter- 
mine the path length difference as a basis for the height calculation. This method 
assumes that the calculated minimum value of regression SS Emin is the center of a 
confidence interval which has a variance o. The variance can be computed using 
the number of fitted parameters m f =3 and the number of points in the fit interval 
n f =530 with, [260] p. 287 


2 SSEmin 
o E 


= = 4.89 x 107. (3.17) 
(nf —mf) 


5 This value is strongly dependent on the used hardware and the optimization performed with 
regard towards parallelization. Within this work neither specialized hardware nor software 
optimizations have been used. This is up to future work. 
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Figure 3.14 Example of the two-stage fitting routine with a) plot of a measured intensity 
signal at \., with a selection of fit curves where the OPD is separated by Aö=5nm for 
each iteration and b) corresponding error sum of squares (SSE) for the different Aô and a 
interpolating curve plotted in black where the arrows indicate a magnified plot of the area in 
close proximity of the minimum SSE 


By using the slope calculated in Fig.3.14 b), the variance can be used to determine 
the interval of path length difference with Aé = + 0.27 nm. Using typically N=10 
consecutive measurements, a mean height profile deviation can be calculated with 
Aö and Eq. (3.15) to ar =0.085nm. The intensity of the data at the equaliza- 
tion wavelength influences the resolution as shown in Fig.3.10. Taking this into 
account, a further analysis can be done where the data set of Fig.3.14a) showed a 
relative intensity of I=0.97 arb. units. From the theoretical calculation of Eq. (3.14), 
Fig. 3.10 it can be deduced that this intensity corresponds to a theoretical resolution 
of r fit = 0.087 nm. This is well aligned with the expected experimental value. Fur- 
ther experimental evaluation in the spatial domain was performed in subsection 3.4. 


3.3.4 Error Estimation of the Data Processing 


Besides the previously discussed error influences due to the optical setup, the imple- 
mentation of the analysis algorithm introduces further noise. In order to quantify 
its influence, a simulation and subsequent analysis were conducted. For this pur- 
pose, a simulated height profile with a nominal height of Znom =0 was constructed. 
It was represented by 500 spectra per profile and 10 repetitions. Based on a typi- 
cal equalization wavelength of Aeq =562 nm, a set of spectra in close proximity to 
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Aeg was calculated having relative intensities of I9(A = Ae,)=O-1 arb. units with 
AI =0.05 arb. units, Fig. 3.15. The dispersive element was assumed to be of N-BK7 
with tpz =2 mm within this simulation. Every calculated spectrum was subsequently 
obstructed by white, Gaussian noise having a SNR SN Rsimy =20-45 dB in steps 
of ASNR=5dB. In total, this resulted in 630,000 spectra processed by the analysis 
algorithm. 

The analysis was concentrated on two important features, the averaged height 
Zmean aS Well as its corresponding standard deviation over the length of the 500 
spectra per data set ^Zmin as an indicator for the resolution. By analyzing the mean 
height and its standard deviation for the spectra at J) =0.5 arb. units with respect to 
the added noise levels, some initial insight can be drawn, Fig.3.16 a). As expected, 
the mean height is close to zero while the standard deviation is dependent on the 
noise level. A mean value Zmean =—2-9 x 101? nm of across all noise levels could 
be measured. On a closer look, the standard deviation of the height measurement, 
which can be regarded as a measure for the resolution, can be fitted by an inverse 
logarithmic model, which corresponds well with the logarithmic scale of the noise 
levels. Using this fitted equation, the influence of the used algorithms could be 
estimated. 
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Figure 3.15 Example of the simulated data set for the characterization of the influence due 
to noise from the data processing routines where the curves are separated by AJ =0.05 in a 
range of /o=0-1 


It is known from the analysis of the spectral intensity signal of the light source, 
Fig. 3.8, that the typical noise is 21.25dB. By evaluating the fitted curve of the 
standard deviation of zmean, a value of AZmin(21.25 dB) =0.16nm can be found. 
As this value is significantly apart from the calculated minimal resolution for this 
noise level of r ¢;;=0.029 nm, a further investigation was performed. The power P 
of the averaged height over the range of S=500 spectra was used to calculate the 
signal-to-noise ration SN Rmeas of the height profile z using 
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Figure 3.16 Depiction of the analyzed simulated data in relation to the different noise levels 

at a relative intensity of Jọ =0.5 arb. units with a) mean height data of averaged profiles Zmean 

and related standard deviation Azmjn where the nominal value was Znom = 0 as well as the 

slope of Azmin with a fitted inverse logarithmic relationship with respect to the added noise; 

b) Plot of the measured SNR from the power of Zmean over the length of 500 analyzed spectra 

per data set with respect to the simulated noise levels displayed on top of a plot of the effective, 
alg 


noise-dependent influence of the algorithm on the resolution Az,,;, 


N 
1 
P=<: i 3.18 
5 2 izl (3.18) 
SN Rmeas = 10 * log P; (3.19) 


where it was compared to the initially simulated noise SN Rsimy, Fig.3.16b). It 
can be seen that the relation follows a linear slope as expected, but a general off- 
set of about 5dB is present. Under consideration of this offset, the resolution for 
the typical noise of the system in an experimental situation was re-calculated as 
AZmin (21.25 dB)=0.088nm. As value takes the complete processing of data into 
account, it is supposed to be the expected minimal resolution of the system at this 
noise level. A subtraction of the noise-related resolution according to Eq.(3.14) 
from the resolution calculated in this simulation— see bottom plot of Fig. 3.16 a)— 
results in the effective influence of the algorithm on the resolution Agi bottom 
plot of Fig. 3.16b). Obviously, it is dependent on the effective noise of the system. 
According to Eq. (3.14), the resolution AZmin of the DE-LCI setup is also dependent 
on the relative intensity at the equalization wavelength AJ. Consequently, the anal- 
ysis of the height measurement, the standard deviation and the SSE in relation to the 
relative intensity. In correspondence to the finding of Fig. 3.16b), this investigation 
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was performed at a noise level of 25 dB, Fig. 3.17 a). From the height measurement 
it can be derived that the measured values are normally distributed around the nom- 
inal value Znom =0. More interestingly, the values for the resolution with respect to 
the relative intensity can be fitted using the already derived relationship from Eq. 
(3.14). This fit supports the finding that the algorithm has a significant influence on 
the resolution which prevents it from achieving a minimal theoretical resolution of 
rfir= 0.029 nm. 
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Figure 3.17 a) Plot of the measured resolution Az, in with respect to the relative intensity at 
the equalization wavelength AZ at a simulated noise level of 25 dB and b) three-dimensional 
representation of the influence of simulated noise and change in relative intensity on the height 
resolution 


As a result of the simulations, the noise behavior of the data processing routines 
can be plotted against all three influences, the relative intensity, the original SNR 
level as well as the standard deviation of the measured average height, Fig. 3.17 
b). From this plot, the described trends can be analyzed by direct comparison. The 
mean standard deviation for the height estimation follows an inverse logarithmic 
relationship with an increase in SNR. Simultaneously, the expected relationship to 
the relative intensity AJ according to Eq. (3.14) is visible. 

In typical measurement scenarios, the characterization of height profiles is 
desired in one lateral dimension or even in an areal fashion. As other approaches 
have shown, gathering line profile data either by scanning the beam or by mov- 
ing the sample relative to the beam is common. In either case, the accuracy of the 
moving parts has an influence on the accuracy of measurement. In order to evaluate 
this influence, an experiment with a plane mirror on a translation stage (Z812B, 
Thorlabs Inc., USA) was conducted, Fig.3.18 a). The translation stage was moved 
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continuously over a range of 2.5 mm during the experiment along the x-axis. Simul- 
taneously, data was collected with the interferometer. As the same mirrors with 
high flatness were used on the translation stage as well as reference mirror, the 
result should ideally show no height differences between both arms, Fig.3.18 b). 
It is visible in the data that the translation stage introduces a significant amount 
of height differences to the measurements. Over a scanning range of 2.5mm the 
stage oscillates in z-direction up to + 1.5 um. In order to achieve a resolution in the 
nanometer range, such deviations have to be diminished or calibrated. 
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Figure 3.18 a) Principle of a setup for the characterization of axial influences from lateral 
movement where a TLS—translation stage moves a SM— sample mirror in the x-dimension. 
The sample arm beam is recombined with the beam from the RM—reference mirror using 
a BS—beamsplitter in order to characterize the errors introduced as height changes in the 
z-dimension interferometrically and b) Result of the translation stage evaluation showing the 
introduced, oscillating height error over a scanning range in the x-dimension 


3.4 Two-Dimensional Approach and Characterization 


The estimations of section 3.2 demonstrated the capability to perform high-resolution 
profilometry with the dispersion-encoded low-coherence approach. But they also 
demonstrated the scanning mechanism’s influence on deviations of the result. In 
order to avoid these and other unwanted deviations from e.g. thermal or vibrational 
influences, an imaging approach was developed to gather two-dimensional data 
without scanning, Fig.3.19 a). In this configuration the reference arm is composed 
of an element with known dispersion (here Schott N-BK7, tpg =2000 um) and a 
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Figure 3.19 a) Experimental setup with WLS—white light source, BS—beamsplitter, DE— 
dispersive element (having the thickness tpg and the refractive index nPE(,)), REF— 
reference mirror, SMP—sample profile including the points zı (x1, y1) and z2(x2, y2) which 
are imaged with a given magnification M (typically M=1.3 or 4) by the LI—imaging con- 
figuration, relayed by a FM—folding mirror onto the slit of the IMSPEC—imaging spec- 
trometer as magnified points zi (x1, y1) and za (x2, y2) and a detailed view of the same in b) 
with SPT—measurement spot, SLT—slit, L2—collimating lens, GRT—grating, L3—imaging 
optics which is used to realize an internal magification (refer also to the appendix in the Elec- 
tronic Supplementary Material (ESM)) and CAM—camera where the spectral information 
for every point on the line in x-dimension is recorded 


plain mirror. The sample arm holds a sample with a varying height profile along 
the x-y plane noted with z(x, y). The recombined collimated light from sample and 
reference arm is imaged on the slit of an imaging spectrometer. Within this spec- 
trometer the light is spectrally decomposed and imaged onto the two-dimensional 
complementary metal-oxide-semiconductor (CMOS) array of a camera, Fig. 3.19 
b). In contrast to a single-line detector of a standard spectrometer, this configura- 
tion enables the detection of spectra at every point on a line in the x-dimension of 
the measurement spot. The information is only selected from one position in the 
y-direction which means that the acquisition of height profiles along a single line 
at once becomes possible. Following this approach, the recorded signal enables the 
detection of spectral interferograms along a line of the x-dimension which can be 
described analogous to Eq. (3.4) and (3.6) with 


I(x, A) = DA) [1 + cos g(x, %)] (3.20) 
with p(x, 4) = 27 Faz@)-Ilne-3@), (3.21) 


where /o(A) is the initial spectral intensity before the beamsplitter and (å, x) is the 
absolute phase of the signal at every point in the x-dimension which is dependent 
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Table 3.3 Calculated parameters and components for the designed imaging spectrometer. 


Name Value Component 

Wavelength range AA 447-780 nm - 

Slit 10 um Thorlabs SIORD 

Collimation lens f3 25 mm Thorlabs AC127-025-A-ML 

Grating 300 lines / mm, 500 nm blaze | Thorlabs GR25-0305 

Off-axis parabolic mirror fa | 101.6 mm Thorlabs MPD149-GO1 

CMOS camera 2048 x 2048 px / 11.8 x 11.8 | Basler acA2040-90um-NIR 
mm 


on the OPD between both arms denoted with ö(x) and the wavelength A. The signal 
detected by the camera of the imaging spectrometer is composed of stacked spectral 
interferograms where the resolution of the x-dimension is dependent on the magni- 
fication M, of the interferometer and the spectrometer construction, see Fig. 3.19. 
Typically, magnifications M=1.3 and 4 where used within this work. For all inves- 
tigations within this work, an imaging spectrometer with the following parameters 
was designed and built, Tab.3.3. A detailed calculation including ray-tracing and 
optimization of optical components for the designed imaging spectrometer as well 
as of the calibration methods can be found in the appendix in the Electronic Sup- 
plementary Material (ESM). 


3.4.1 Height Standard Evaluation 


Additionally, experiments have been conducted to evaluate the accuracy of the 
developed system for the determination of small height steps. For this purpose, a Si- 
based step standard (VS 0.10, Simetrics GmbH, Germany) was examined, Fig. 3.20. 
The results of this examination revealed a good ability to resolve nm-sized height 
steps with a measured height of (101.8+0.1)nm which is in good agreement to 
the nominal value of (100+ 7) nm quoted by the manufacturer. The corresponding 
RMS error with regard to the nominal value was 1.1 nm. The roughness, Ra, of the Si 
surface could be measured with 0.8 nm which scales with a factor of about 8.7 to the 
roughness measure Rt=7.0nm, [261]. This is within the range of 6- 10nm quoted 
by the manufacturer, [262]. The recorded and measured profiles show bat-wing 
effects at the sharp edges, [263]. The measurement error increases in the regions 
of these effects due to diffraction and deflections. It was visible that deviations of 
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Figure 3.20 Results of the measurement of a (100 +7) nm nominal height Si-standard with a) 
simulated spectral interference signal over a spectral range of 333 nm and lateral dimension 
of 450 um with the equalization wavelength Aeq marked and b) corresponding measured 
spectral interferences data with visible intensity modulations due to diffraction (marked with 
red ellipses) and c) calculated mean height profile from the raw data with diffraction-induced 
bat-wing effects at the sharp edges (marked with red ellipses) having a measured height of 
(101.8 £0.1)nm 


up to 20nm occur, marked with red ellipses in Fig. 3.20 b) and c). These deviations 
were attributed to diffraction effects visible as additional intensity modulation in the 
spectral interference raw data, Fig. 3.20 b). For calibration purposes, the oscillations 
of the diffraction can be modeled as Fourier filtering by the aperture of the capturing 
optical system, [263]. In relation to the simulated raw data, Fig. 3.20 a), it was visible 
that not only diffraction occurs, but other distortions as well. In case of a flat, properly 
aligned sample, the spatial distribution of the maxima and minima is parallel to the 
x-axis of the plot, see Fig. 3.20 a). It can be seen in the actual measured data, that this 
was not the case, Fig. 3.20 b). This was the result of a slight tilt of the sample (about 
0.11 nm/um) in relation to the sample arm. It was corrected during post-processing 
of the final measured profiles assuming a linear tilt. 

On the same standard, a single edge as well as a series of steps having a 
pitch of 250um were studied, Fig.3.21. For both sample positions heights of 
(104.35 + 0.11) nm and (99.88 +0.11)nm were measured respectively. It is clearly 
visible that both of these features show the same diffraction effect, which confirms 
that it is independent of the sample position but dependent on the feature size and 
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Figure 3.21 Plot of measured structures of the Si-based height standard showing diffraction 
effects with a) single edge having a mean height of (104.35 + 0.11)nm and b) series of steps 
with a mean height of (99.88 +0.11)nm and diffraction effects of up to 50 um lateral size 
from each edge which influence the step width significantly, see inset 


slope. Fig. 3.21 a) leads to the note that the effect has a length of influence of about 
li=50um into the profile. When measuring structures with lateral feature sizes 
smaller than 2 - l; information of these structures can be obscured. It can be seen in 
the measured profiles of multiple successive steps with a width of only 125 um that 
an evaluation of e.g. roughness is influenced by this effect, see inset Fig. 3.21 b). 


3.4.2 Repeatability and Resolution Characterization 


The error of the system can be analyzed by utilizing two measures where one is 
the repeatability, defined by the standard deviation o,(x) of multiple profiles z; (x) 
gathered in a short time frame. The second measure is the resolution, calculated as 
the standard deviation ^Zmin of a feature such as height h;, Fig.3.22. a). In order 
to analyze the repeatability, the structure presented in Fig. 3.20 c) was measured 
N =10 times in a row without any other delay than the acquisition and data transfer 
time. The analysis of the standard deviation of the profiles with respect to their 
mean, (x), allows one to conclude on the repeatability, 


N 


oz(x) = — E haw- 0). (3.22) 


i=1 
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Figure 3.22 a) Depiction of the calculation of repeatability as the standard deviation o; (x) of 
multiple profiles z; (x) according to Eq. (3.22) as well as the resolution Az,,;„ as the standard 
deviation of the feature height h; according to Eq. (3.23) and b) plot of the spatially resolved 
repeatability of the nm-sized height standard of Fig. 3.20 c) where the impact of diffraction is 
visible at x = 150-210 um and x = 300-350 um with an inset to visualize the magnitude 
of o; (x) between x = 50-150 um 


It can be noticed by analyzing the standard deviation in relation to the lateral dimen- 
sion 0;(x) that the error significantly increased due to the sharp edges and the 
diffraction effects, Fig. 3.22 b). In order to estimate the repeatability of the setup in 
this configuration, the standard deviation was evaluated without the data points that 
were affected by defraction (x =150-210um and x = 300-350 um). A mean 
value of o; (x)=0.13 nm was calculated. The value of 5;=0.13nm, measured on a 
low scattering nm-height standard, is expected to be the lower limit of the setup, 
as sections of the sample with disturbing influences were excluded from the calcu- 
lation. In order to characterize the repeatability of the setup further, an experiment 
was designed where the path length difference was altered with a translatable sam- 
ple, Fig. 3.23 a). The designed sample was a flat piece of a Si-wafer (15x 15 mm?) 
diced and mounted on a piezo stack (A = 1.8 um per 100 V, PA3JEW, Thorlabs Inc., 
USA) which was attached to a bulk glass substrate. During the experiment, the 
voltage of the piezo was discretely controlled (power supply QL355TP, Aim and 
Thurlby Thandar Instruments, United Kingdom) and monitored (digital multimeter 
DM3068, RIGOL Technology Co., Ltd, China) to adjust the path length differ- 
ence in defined steps of Aö=0.2nm. For every adjusted step, N= 10 consecutive 
measurements were taken, Fig.3.23 b). The data for the three pictured positions 
demonstrates that the setup is capable of resolving steps of 0.2 nm as the standard 
deviations hardly overlap. The mean standard deviation of every position was cal- 
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Figure 3.23 a) Depiction of the modified setup to characterize the repeatability where the 
incoming electric field E;n is split by aBS— cube beamsplitter, while 50 % of the light passes 
the DE—dispersive element (N-BK7, tpg =2 mm) and gets reflected on the REF—-reference 
mirror. In the second arm, light gets reflected of the designed sample which was a piece of a 
Si—silicon wafer attached to a PZO—piezo stack (A=1.8 um per 100 V, PA3JEW, Thorlabs 
Inc., USA) which was mounted on a bulk glass substrate. The voltage of the PZO was remotely 
controlled in order to change the path length difference in steps of Ad =0.2 nm where b) plot of 
the recorded averaged surface profiles along the x-dimension and the corresponding standard 
deviation shaded around the slopes for three different positions of the piezo stack 


culated as o (Aö = 0)=0.11 nm, o (Ad = 0.2) =0.13nm, o (Ad = 0.4) =0.10nm. 
These values correspond well with the measurements on the Si height standard 
shown previously and demonstrate the limit of the current setup with regard to 
stability. 

The gathered data was further analyzed using a two-sample Student’s f-test in 
order to statistically evaluate if the measured averaged slopes are significantly dif- 
ferent from each other, [260]. A mean value of f}2 = 3.282 and f3 =3.619 was calcu- 
lated between the profiles with Aö=0 / Aö=0.2nm and Aö=0.2nm / Ad=0.4nm 
respectively. In contrast to the value of t, =3.250 (with n=9 at 99% probability) it 
was found that ? > t; which rejects the null hypothesis. Consequently, the averaged 
profiles are significantly different in this experiment. Furthermore, this test can be 
used to estimate the minimal repeatability for the case where ?=t,. It was found to be 
Admin =0.12 nm for a 95 % probability and Aômin =0.18 nm for a 99 % probability. 

The creation and measurement of samples which are close to the proposed reso- 
lution of the setup is complex. For the calibration of AFM instruments, height step 
samples exist which have heights in the size of one atomic layer of silicon, [264, 
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265]. However, the steps on these samples usually have widths below 1 um which 
are not resolvable with the presented approach. While the repeatability is a measure 
for temporal fluctuations that occur from one measurement to the other, the ability 
to resolve structures along the spatial domain (here denoted as the x-coordinate) 
is independent from these fluctuations. A measure for the resolution can be found 
in the standard deviation Az in of feature sizes such as the height of structures h; 
relative to the mean height of multiple measured features, h. It can be assumed that 
in-between short time frames of the acquisition time of single data sets the sample 
does not change, 


N 
1 2 
Azmin = NI = 2 (hi = h) n (3.23) 


In case of the nm-sized, Si-height standard, the height was measured as the difference 
between the two base levels, x; = 100—150 um and x2 = 350—400 um, and the 
top plateau of the step at x3 = 225-275 um. The quadratic mean of AZmin for 20 
measured heights, and therefore the resolution, was found to be 0.1 nm. The standard 
deviation of the feature size represents a cumulative measure for the resolution which 
includes influences of the optical setup, the electronics, the calibration routines and 
the data processing alike. During the data analysis of the presented results it became 
obvious that the difference between the calculated minimal resolution of 0.088 nm, 
see bottom plot of Fig. 3.16 a), and the minimal measured resolution of 0.1 nm are 
partly due to data processing routines. As the recorded profiles usually were tilted 
by a minor degree, an appropriate tilt correction was performed based on the linear 
fit of every captured surface profile. Although the tilt correction was optimized, a 
minor influence on the standard deviation cannot be excluded. 

As a supplement, a step with a nominal height of Anom =3 nm was created using 
standard semiconductor procedures, Fig.3.24a). A silicon substrate was prepared 
with polyimide tape and sputtered with ITO of about 3 nm thickness. Afterwards, 
the tape was removed and a second sputtering step using titanium was performed. 
Due to the initial application of tape, the sample has a defined height difference 
which is maintained after subsequent generation of the titanium layer. This layer 
was applied in order to generate a uniformly, high reflecting surface and prevent 
any thin-film interferences which a sole layer of ITO would have caused. 

The created sample was analyzed using the setup similar to the measurements 
on the height standard where N=10 consecutive measurements were taken and 
averaged, Fig. 3.24b). Using the same scheme as before, the height of the step was 
determined as hstp =3.45 +0.19 nm. The result fortifies the claim that the proposed 
DE-LCI method is capable of measuring surface profiles with sub-nm resolution. 
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Figure 3.24 a) Depiction of the processing chain to fabricate a sample with a step of 3 nm 
nominal height where (I) polyimide tape is placed as mask on a Si-substrate and a layer of 
ITO is sputtered on the substrate having the desired thickness. (II) The removal of the tape 
finishes the step formation. (III) In order to generate a sample with a uniformly reflecting 
surface, a layer of 40 nm Ti is sputtered onto the sample which maintains the step. b) Plot of 
the averaged height profile of the step with hy,» =3.45 +0.19 nm using the DE-LCI setup 


In particular, this is supported by the repeatability of o =0.12 nm and the resolution 
of the height measurement AZmin =0.19 nm. 

On the basis of the calculated resolution, the DR of the setup was calculated 
with the measurement range of Az=79.91 um and the resolution of AZmin =0.1 nm 
as DR=7.99 x 105. Compared to the latest findings of other areal profilometer 
approaches such as of Reichold et al., [119], the achieved dynamic range is about 
5.8 times higher. 


3.4.3 Edge Effects 


The occurrence of edge effects on sharp edges is caused by different sources. While 
diffraction and scattering play an important role, also filtering effects of the aperture, 
shadowing from the sample’s steep slopes as well as interferometric mixing of 
components of the different height levels due to the lateral resolution are relevant. 
A rough estimation of the influence of the aperture revealed that its contribution 
is neglectable, such that it is assumed that diffraction due to the spatial coherence 
properties of the light source are dominant. As these influences are mixed with 
an unknown ratio, no single model can be used to filter the signal appropriately. 
For this reason, a deconvolution according to the WIENER approach was found to 
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Figure 3.25 Depiction of the WIENER-based deconvolution routine with a) components of 
the process where U represents the measured signal (see Fig.3.21a), Z is the ideal profile 
and H is the estimated impulse response function as well as b) comparison of the measured 
profile U with the deconvoluted profile Z’ 


be suitable, [266], Fig.3.25a). This approach typically assumes that besides the 
measured profile U, a response function of the system H is known. Based on these 
functions, a deconvolution for the ideal filtered profile Z’ can be performed. Using 
this notation, a measured profile is the convolution of the ideal profile with the 
response function® 

U=H.-Z. (3.24) 


Following WIENERS idea, a wild-card function G acting on the measured signal is 
utilized to minimize the error between the ideal, but unknown profile Z and the 
deconvoluted signal Z’ 


Z'=G-U (3.25) 
min(Z’ — Z). (3.26) 


As the response function in the typical use case of the WIENER deconvolution is 
known, the wild-card function can be constructed of H 


H* 


Ge (3.27) 
IHI? + (zur) 


6 These functions are assumed to be handled in Fourier-space. 
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In the case investigated here, a modification has to be performed as H is unknown 
but Z is known. Hence, in order to perform the deconvolution, a wild-card function 
Q has to be introduced based on Z 


Zz 
Q = — (3.28) 
ZI? + (sie) 
which can be used to compute the response function 
H=Q.-U. (3.29) 


With the aid of this response function, finally the deconvolution of the original 
measured signal U can be done using Eq. (3.27) and (3.25) resulting in the filtered 
signal Z’. This computation was performed on the measured data of a single silicon 
edge, originally presented in Fig.3.21a), Fig.3.25b). It is obvious that the edge 
effects can be filtered well. The outlined procedure can be made part of a calibration 
routine when measuring similar structures repeatedly. 

Additionally, the dependency of edge effects with regards to the used light source 
was evaluated. For this purpose, a supercontinuum white-light source (SC) (with 
Aà =380-1100nm) and a laser-driven plasma light source (LDP) (with AA =200- 
1100nm) were used in comparison. Due to the limited coherence of the broadband 
light sources the interference contrast of the wavelength-integrated signal forms a 
considerable envelope function which can be detected by using a spectrometer. This 
data is the basis for the calculation of the coherence length. The coherence length 
le can be determined exactly as the integral of the area under normalized degree of 
coherence. Usually, an approximation can be found by estimating the width of the 
interferogram at an intensity of !/e, [105]. The respective coherence lengths were 
determined as le =1.81 um for the LDP light source and le =1.58 um for the SC 
light source. Consequently, the same step standard of 100nm nominal height was 
measured with both light sources. While the measured height was comparable with 
both light sources, the behavior on the edge of the standard showed differences, 
Fig. 3.26. The analysis reveals that using a LDP source, a much steeper slope can be 
measured. A linear approximation results in a slope of 0.94 um/mm for the SC source 
and a slope of 1.82 um/mm for the LDP source. Furthermore, it can be seen that the 
formation of edge effects is stronger in the measurements of the SC light source. As 
the coherence length is very similar for both light sources, a possible explanation is 
that the spatial coherence of SC light source is higher. Due to its operation principle, 
highly spatially coherent light of a microchip laser is used to generate a broadened 
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Figure 3.26 Result of the measured slopes using a SC versus a LDP light source 


spectrum . In contrast, the LDP light sources generates its broad spectrum from a 
random process both temporally and spatially. 

Further knowledge of the precise spatial coherence properties of the used mea- 
surement light source as well as of possible sample geometries can help to develop 
appropriate filter models for the correction methods demonstrated above. The spa- 
tial coherence properties of the light source and resulting effects such as diffraction 
and scattering are the main limiting factors for measurements requiring high lateral 
resolutions. 


3.4.4 Roughness Evaluation 


As surface quality and roughness in particular can be essential for the function 
of technical products and components, their in-line assessment is of high interest, 
[21, 82, 267-269]. Various norms and guidelines exist for a number of established 
measurement technologies such as tactile profilometers, confocal microscopes and 
AFM, [17-19]. Most commonly, quantification methods based on the distribution 
of heights like the averaged roughness Ra and the root-mean-square roughness Rq 
are determined. 

As the demonstrated DE-LCI approach is capable of capturing precise height 
profiles over a large lateral measurement range of several hundred micrometers, an 
application for roughness evaluation of samples is possible. For the initial qualifi- 
cation, a PTB-traceable surface roughness standard (KNT 20170/3 superfine, Halle 
GmbH, Germany) was analyzed, Fig. 3.27. According to the specifications of PTB, 
surface profiles were captured on the slices A, B and C in segments of (I), (ID and 
(IID) in each slice. Finally, the results of the roughness measurement of all segments 
were averaged. According to ISO 4288-1996, surface quality is assessed by the appli- 
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Figure 3.27 Schema of the measured PTB-traceable surface roughness standard where mea- 
surements are taken in the slices A, B, and C with segments of (I), (II) and (IID all having the 
same evaluation length of L.=5-A. 


cation of appropriate processing steps and filters, [17]. Following these processing 
steps, the captured surface profiles were form corrected (A,-filtering according to 
ISO 3274-1996 [270]) and referenced to their respective mean values in order to 
generate the so-called primary profile, Fig.3.28a). Based on this profile, a Gaus- 
sian filter was applied to calculate the waviness and the roughness of the primary 
profile, Fig. 3.28 b) and c) respectively. Here, the application of the low-pass filter 
results in the waviness profile, while the application of the high-pass filter results in 
the roughness profile. The cut-off wavelength A, is determined in correspondence 
with the evaluation length of the profile which should be approximately 5-A.. In 
case of the roughness standard, the captured evaluation length was Z.=1.25mm 
as the filter length was A; =250 um. In order to exclude possible deviations on the 
roughness calculation that are due to filtering effects, the actual measured length 
was expanded equally by 1/2-A.. at the beginning and the end of the profile. By 
applying this methodology to the data captured from the roughness standard, values 
of Ra=(21.15+0.8)nm and Rq=(26.58+1.0)nm were calculated. The value of 
Ra is within the specifications given by the PTB calibration which measured a mean 
value of Ra=(22.4+0.5) nm with an uncertainty of +5 %. As a second measure of 
comparison, a roughness evaluation using a confocal microscope was performed. 
Using the same probing scheme highlighted in Fig. 3.27, data was captured in stitch- 
ing mode as the 50x magnification objective was needed to achieve sufficient axial 
resolution and is restricted to a lateral measurement range of 220 um. The captured 
data was filtered for noise using a A, =8 um Gaussian low-pass filter and analyzed 
in the same way as the interferometric data. The roughness parameters were cal- 
culated with Ra= (21.42 +0.6)nm and Rq=(26.81 + 0.7) nm. These values support 
the interferometric values within the respective errors. 
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Figure 3.28 Result of the measurements on a PTB-traceable surface roughness standard with 
a) primary, form corrected profile, b) waviness profile (Gaussian low-pass filtration of a)) and 
c) roughness profile (Gaussian high-pass filtration of a)) both using a cut-off wavelength of 
A.=250um 


It is known that roughness parameters which are based on amplitude values of the 
height distribution, such as Ra and Rq, are not directly comparable in different mea- 
surement techniques, [267]. This is particularly due to the different spatial bandwidth 
limitations of the individual techniques, [271]. For the purpose of enhancing compa- 
rability, alternative methods to determine the RMS roughness were established. The 
most commonly used ones are based on the determination of the integral of either 
the auto-correlation function (ACF) or the power-spectral-density function (PSDf) 
of a profile, [272-274]. In order to perform a comparison of optically measured sur- 
faces with tactile measured ones, data from the same standard sample was captured 
using both methods, Fig. 3.29 a). It can be seen from the roughness profiles that both 
data sets show similar amplitudes. However, the data from the tactile measurement 
shows significant noise which could be measured with + 2.6nm. This noise can be 
attributed to the typical noise of a tactile measurement system consisting of mechan- 
ical and electronic components, [275]. With the expected roughness levels of about 
20 nm, this noise has an influence on the roughness measurements. For the com- 
parison of both methods, the respective auto-correlation functions were determined 
using a Fourier-based approach, Fig.3.29 b). The acquired data was subsequently 
fitted using a Gaussian approximation in order to perform further analysis, [276]. A 
fundamental analysis is the determination of the correlation length tx. Most com- 
monly, it is measured as the distance x where the approximated ACF reaches a value 
of g(x) = !/e, [277]. While the analysis of the tactile measurement yielded in a value 
of Ty = 11.37 um, the nine analyzed measurements of the interferometric evaluation 
resulted in a value of t, =(11.69+1.0) um. Obviously, the measurements of the 
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Figure 3.29 Results of the comparative roughness evaluation with a) roughness profiles of 
the PTB-calibrated height standard from an interferometric and tactile measurement as well as 
b) calculated auto-correlation functions and appropriate fits using a Gaussian approximation 


different technical approaches show a high similarity within the standard deviation 
of the measurement. Furthermore, the amplitude of the fitted ACF can be inter- 
preted as the squared RMS roughness of the measured data. The comparison of 
these properties reveals the influence of noise on the tactile measurement. A value 
of Orms =27.6nm was calculated for the tactile measurement, while a mean value of 
Orms =(21.46+ 1.8)nm was calculated for the interferometric measurements. Tak- 
ing the noise of the tactile measurements into account, both measurements are very 
similar within the respective error bars. 

In addition to measurements on a standard, the characterization of an industry- 
relevant configuration, in particular an aluminum mirror coating on a float glass 
substrate (Layertec GmbH, Mellingen, Germany), was performed. The coating was 
applied using magnetron sputtering on one half of the circular substrate for evalua- 
tion purposes, Fig. 3.30. Using the averaged data of three measurements on different 
positions along the coating edge with ten measurements at every position, a mean 
height of zaz =(99.47 + 0.12) nm was measured. In comparison with the data ana- 
lyzed for the height standard, the edge of the mirror was less steep which resulted in 
significantly lower edge effects. Furthermore, the evaluation of sub-nm roughness 
differences as a part of production accompanying characterization is of interest. For 
this purpose, data from the float glass substrate as well as from the aluminum coated 
part of the mirror was analyzed with the described roughness methodology, Fig. 3.31. 
The separation of roughness (III) and waviness (II) from the primary form-corrected 
profile (I) reveals a distinct difference. While the waviness of both profiles is in the 
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Figure 3.30 Plot of a measured aluminum mirror edge on a float glass substrate with height 
profile as mean value of ten measurements on one position along the edge having a height of 
har =(99.47+0.12) nm 


same order of magnitude, the roughness of the aluminum coated part of the sample 
is larger. The mathematical analysis resulted in values of Ra = (0.27 + 0.01) nm and 
Rq=(0.35 + 0.01) nm for the substrate area. The roughness of the aluminum coated 
part was Ra=(0.38 + 0.02) nm and Rq= (0.47 + 0.02) nm. Compared to Ra=0.31 nm 
for the substrate and Ra=0.40nm for the coating, which are quoted by the manu- 
facturer, the measured values correspond well. 

In contrast to other technologies such as AFM, scanning white-light interfer- 
ometric microscopy and confocal microscopy, DE-LCI is able to capture data for 
roughness evaluation on a large lateral measurement range of a few mm in one single 
data acquisition. The lack of the necessity to scan a sample eliminates problems of 
stitching, vibration and speed. 


3.4.5 High-Dynamic Range Measurements 


In order to measure the performance of the setup with a high-dynamic range where 
the measurement range is > 10 um while the achievable height resolution should still 
be in the nm-range, a precision-turned height standard (EN14-3, PTB, Germany) 
was examined. The standard provides grooves of defined heights with steps of 
1, 5 and 20um which were subject to a series of measurements, Fig.3.32. The 
recorded data includes measured steps of (971.26 +0.31), (4951.40 +0.28) and 
(19924.00 + 0.36)nm. This results in an overall averaged RMS error of 26.9 nm 
with regard to a measurement on a tactile profilometer. According to the calibration 
certificate of the standard, these values are within the quoted uncertainty for the 
nominal height steps of + 33 nm. Furthermore, the high axial resolution leads to 
the ability to capture roughness data in the nm-range on all height steps, see inset 
in Fig.3.32a). A RMS value of Rq=26.7nm was calculated. Some edge effects 
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Figure 3.31 Result of roughness measurements on a) an aluminum mirror surface and b) a 
float glass substrate with (I) primary form corrected profile, (II) waviness profile (Gaussian 
low-pass filtration of (I)) and (III) roughness profile (Gaussian high-pass filtration of (I)) using 
a cut-off wavelength of A, =25 um (mean of ten measurements) 


and noise occur in slopes of the steps. In the current optical design, that uses a 
NA=0.06 imaging system, a large lateral measurement range could be covered but 
data on the slopes with an 70° angle could not be gathered reliably. Depending on 
the application, the setup can be optimized to increase the sensitivity on these parts 
of the sample. Reference measurements with a tactile profilometer confirmed these 
heights but emphasized the fact that the transitions between the different levels are 
formed by segments of 70°, Fig.3.32 b). In contrast, the tactile profilometer is able 
to gather a much higher number of data points in these areas while having a lower 
overall resolution. 

Additionally, the repeatability according to Eq. (3.22) was analyzed by investi- 
gating N= 10 consecutive measurements of the profile. The sample showed a slightly 
increased averaged standard deviation of oz =0.52 nm with respect to the measure- 
ments on a low-scattering silicon sample, Fig.3.22 b). This can be attributed to 
influences of noise due to diffraction, scattering and other effects affecting the mea- 
sured repeatability. In a comparative manner to the analysis of the silicon height 
standard, see subsection 3.4.2, the heights on the zm-sized standard, Fig. 3.32 a), 
were evaluated for 10 measurements as the difference of the two closest base levels 
(z = 0) to the particular step. From the data, quadratic means of AZmin for the three 
steps with nominal heights of 1, 5 and 20 um were calculated as Azyinı =0.31 nm, 
AZmins =9.28 nm and AZmin20 =0.36 nm respectively. The result shows that the res- 
olution is not dependent on the size of the measured step. 
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Figure 3.32 Plot of an averaged line profile where recorded depths of (971.26 0.31), 
(4951.40 +0.28) and (19924.00 +0.36)nm could be measured with a mean RMS error of 
26.9nm with respect to a measurement on a tactile profilometer while having the ability to 
capture roughness information which is shown in the inset where a value of Rq=26.7nm 
was calculated as well as b) overlay of a dataset from the same sample taken with a tactile 
profilometer which shows significantly better capabilities to capture data on steep edges 


In contrast to the silicon standard, the um-sized height standard has a significantly 
higher roughness which leads to scattering. In consequence, the measurements on 
this sample were affected by noise which led to a number of outlier data points. 
The implemented post-processing routines took these outliers into account and cor- 
rected them. The outlier correction was performed by the detection of rising edges 
using data of the first derivative of the profile with respect to the x-coordinate in 
combination with the correction of the difference between the outlier and the mean 
value of five previous data points. In the analysis of the step heights and its stan- 
dard deviation it could be detected that the outlier correction scheme influences 
the profile on a sub-nm level. As outliers occur on different spatial positions for 
consecutive measurements, a higher standard deviation was measured opposing to 
measurements where less outliers occurred, as on the silicon height standard. 


3.4.6 Dual-Channel Approach 


As surface profile evaluation is crucial at different industrial processing steps, tomo- 
graphic evaluation of structures becomes also interesting. For this purpose, the DE- 
LCI approach was adapted to the NIR to perform tomographic examinations of 
Si-based structures, Fig.3.33 a). The setup was designed to work in a dual-channel 
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configuration where a broadband light source illuminates the sample and a dichroic 
mirror at 800nm separates the recombined light for the two analysis channels. In 
contrast to conventional DE-LCI, the light of the NIR spectral range holds informa- 
tion from inside the sample while light of the VIS spectral range holds only surface 
information. This is of course valid for samples like e.g. silicon which are transmis- 
sive in the NIR range but not in the VIS. Consequently, an appropriate imaging spec- 
trometer was calculated and designed for a spectral range of AA =(1133-1251)nm 
based on an indium gallium arsenide (InGaAs) camera (Bobcat 640, Xenics Ltd., 
Belgium). A detailed description of assumptions, parameters and components is 
given in the appendix in the Electronic Supplementary Material (ESM). 
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Figure 3.33 a) Schema of the extended optical setup for dual-channel interferometry with 
WLS—white-light source, BS—beam splitter, SMP—sample having a thickness tsmp and a 
refractive index n°”"?(A) where zı (x1, yı) and z2(x2, y2) are two points, one on the surface, 
one on the back side, REF—fixed reference mirror, DE—dispersive element with the thickness 
tpg and nP” (X), L1—lens to image the sample with a given magnification M (typically M 
= 1.3 or 4), HP—high-pass filter @800nm to reflect the VIS part of the spectrum where 
VISSPEC—VIS imaging spectrometer detects the surface information of the magnified point 
z (x1, y1) and FM- -folding mirror relays the NIR part of the spectrum to NIRSPEC—NIR 
imaging spectrometer which detects the depth information of the magnified point 25 (x2, y2) 
as well as b) simulation of possible measurement ranges of materials suitable as dispersive 
elements for NIR investigations 


In order to generate spectral power densities in the spectral range of (1000- 
1400) nm, an amplified supercontinuum light source (ASC) was used, [250]. The 
source utilizes an Yb?* doped photonic crystal fiber as medium for non-linear 
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spectral broadening and as gain medium within a fiber amplifier configuration, 
Fig. 3.34 a). The configuration utilizes a pump diode laser (A pmp =976 nm) to core 
pump the fiber while a passively Q-switched microchip laser (Aseeq =1064 nm) with 
pulse durations of 1.3 ns and a variable repetition rate of up to 20 kHz is used to seed 
the system. The amplification enables several W/nm pulse peak power in the desired 
spectral range and beyond, Fig.3.34 b). The high pulse peak power as well as the 
short pulse duration make the light source interesting for dynamic measurements 
to e.g. observe MEMS movements. Stroboscopic illumination can be envisioned to 
achieve high penetration depths and high temporal resolution’. 
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Figure 3.34 Schema of the optical setup for the generation of amplified supercontinuum 
(ASC) in the NIR range with MCL—microchip Laser, PLD—pump laser diode (A pmp = 976 
nm), Ol—optical isolator, M—mirror, NM—notch mirror, PCFYb—Yb:doped PCF— fiber 
as well as b) the optical spectral pulse peak power density of the light source (both adapted 
from [250]) 


Measurement range and resolution in NIR evaluation 

With its dispersion characteristics, N-BK7 is a very suitable material to be used for 
DE-LCI in the VIS where typical measurement ranges of 79.91 um are achieved 
(tpe =2mm). Due to lower spatial resolutions of cameras used for imaging spec- 
trometers in the NIR the covered spectral range is usually low. Furthermore, most 
materials have a rather flat n(A) slope in this spectral region. Both factors limit the 
possible axial measurement range. Using N-BK7 with tpg =2 mm in the range of 
Ad =(1133-1251)nm, yields in an axial measurement range of Az=1.12 um. A 


7 These aspects were not studied within this work and are subject to further research. 
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possibility to increase the range is the utilization of the sample as dispersive material 
if its thickness is known. Under the assumption that a sample is a silicon wafer with 
a thickness of 100 um, the resulting range would be 7.03 um. Depending on the 
application, this range can be suitable to detect buried marks in a wafer. Similar 
measurement ranges are achievable by substituting N-BK7 for a higher refractive 
glass. The usage of FK51A would enable a range of Az = 1.8 um while SF11 would 
lead to Az=5.8 um assuming a thickness of 2mm. Significantly higher measure- 
ment ranges, which are comparable to the VIS approach, can only be achieved with 
non-glass materials, Fig. 3.33 b). As already discussed, Si can be used to extend the 
range as it can reach Az= 146.04 um for tpg =2 mm. Even higher ranges can be 
observed with materials like gallium arsenide (GaAs) where Az=238.6 um. How- 
ever, it has to be noted that both Si and GaAs are non-transmissive in the visible 
spectral range. Therefore, these materials are not suitable in a dual-channel approach 
where it is desired to simultaneously gather data of both VIS and NIR channels. A 
possible measurement mode would incorporate a highly dispersive material which 
is non-transmissive in the VIS as a reference mirror. In this way, a glass DE can 
be used for measurements in the VIS while the reference mirror acts as a DE for 
the NIR investigations. Alternatively, a transmission mode operation can be imple- 
mented with a material which is transmissive in the VIS and NIR range. A material 
which meets this requirement is zinc selenide (ZnSe). It has high optical transmis- 
sion starting at 550 nm. In the NIR, a measurement range of Az =47.88 um can be 
achieved by the application of a DE with tpg =2mm*. Consequently, this leads to a 
maximal achievable measurement range in the VIS of about Azy zs =605.84 um, if 
the DE is used in transmission mode. All presented measurements in this work have 
been performed using a ZnSe element (tpg =2 mm) in transmission mode. Further- 
more, it has to be noted that the measurement range only describes the ability to 
detect surface height changes on one particular surface. In a tomographic measure- 
ment, the separation of the equalization wavelengths for multiple surfaces is also 
important. A simulation showed that a sample has to have a minimum thickness of 
60 um of silicon in order to capture the equalization wavelengths in both spectral 
channels separately for the particular setup described here. 

While the equations derived in section 3.2 for the estimation of the profile height are 
still applicable for the surface information, the back reflected data from structures 
within the sample follow a different relation. In this configuration, the phase is 
transformed in slight variations regarding Eq. (3.20) with an additional component 
for the samples refractive index n°”? (A) =n°}""? and its thickness tsmp, where n°”? 


8 In a real measurement setup, the measurement range will be extended by the dispersion on 
the sample. The amount of increase is dependent on the sample depth to be analyzed. 
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Consequently, all data captured from the inside of a sample is scaled by the depth- 
dependent refractive index which has to be corrected, if height information should 
be obtained. In order to perform this, blind dispersion compensation techniques 
known from OCT can be utilized, [256-258]. 

Analogous to the measurements and calculations performed previously, see sec- 
tion 3.2, the resolution limit of the NIR approach was characterized. A measurement 
of the noise yielded in an average value of AJ = 16.0dB which was found to be nor- 
mally distributed along the spatial and spectral dimension of the NIR imaging spec- 
trometer. The measured noise was used to calculate the single point resolution limit 
Aô =4.73 nm. This calculation was based on Eq. (3.14) where the wavelength range 
was AA=(1133-1251)nm while the equalization wavelength was Aeg =1189nm 
and the relative normalized intensity at this point was Jo =0.5 arb. units. Under the 
assumption that n=300 points were used for fitting, a resolution of the NIR system 
of r fir =0.27 nm was calculated, referring also to Eq. (3.15). By extrapolating the 
estimated influence of the algorithm for the measured noise AZ, see subsection 3.3.4, 
an influence of 0.22nm can be computed. This leads to an expected resolution of 
0.49 nm for this experiment in the NIR spectral range. In relation to the measure- 
ment range of Az =95.76 um, which can be achieved by using a dispersive element 
of ZnSe with tpg =2 mm, a dynamic-range of DR=1.95 x 10° was calculated. 


Results of tomographic profilometry 

Using the above described setup and configuration of the DE, an experiment for 
the tomographic imaging of a thinned wafer was conducted to capture the surface 
profile as well as the backside profile of the sample in a simultaneous measure- 
ment, Fig. 3.35. The sample was a Si-wafer with a partially coated area which was 
thinned previously. It was mounted on a thin glass substrate. The measurement 
shows the stepped surface profile where the coating height was determined with 
18.70 + 1.42 um. By analyzing the interferometric signal in the NIR channel and 
correcting the measured OPD with the refractive index of silicon, the wafer thickness 
was determined with 84.71 +0.38 wm. The edge of the coating can be identified in 
the VIS as well as in the NIR signal with some significant noise (x= 110—140 um). 
The refractive index of the coating introduced a dispersion-induced deviation in 
the NIR measurement. This deviation was not corrected in this measurement as the 
material composition of the coating was unknown. It is also visible that the VIS 
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Figure 3.35 Result of the profilometry measurement using a dual channel approach where 
the VIS channel shows the surface profile on top of a thinned Si—wafer with a coated area 
and a tomographic profile from the backside of the wafer using the NIR channel of the setup 


signal shows significant overall noise. This is due to the relatively small spectral 
power density of the used ASC source in combination of the low optical transmis- 
sion of ZnSe in this spectral range. This led to a reduced SNR of 16.5 dB for the VIS 
measurement compared to the experiments presented before. Future developments 
will account for this and develop methods to increase the SNR. Furthermore, the 
development of automatic dispersion correction for tomographic data according to 
known approaches will be worked on. 


3.5 Areal Measurement Approaches 


In order to gather areal information, different approaches were developed. While 
two methods were developed only theoretically, one approach was implemented, 
characterized and tested. 


3.5.1 Translation-Based Areal Information 


In this implemented approach, information was obtained by constantly translating 
the lens LZ in order to image different parts of the sample on the slit, as already 
depicted in Fig. 3.1. In consequence, a stack of two-dimensional line profiles were 
gathered and analyzed in order to receive three-dimensional data. As the imaging 
lens, placed after the light of both interferometric arms, was recombined and the 
sample was not moved between measurements, negative influences on the measure- 
ments were kept to a minimum. 
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Figure 3.36 a) Simulated stack of two-dimensional spectra gathered by the translation of 
the imaging lens L1 in the y-direction to capture 3D information and b) plot of the three- 
dimensional surface of a precision-turned groove standard (Gaussian filter applied to reduce 
edge effects for display purposes) with measured depths of (971.26 +0.31), (4951.40 + 0.28) 
and (19924.00 + 0.36)nm 


Three-dimensional information of the precision-turned height standard used for 
the high-dynamic range evaluation, Fig.3.32a), was gathered in steps of 25 um 
along the y-direction and a rather small magnification to enable a lateral measure- 
ment range in the x-direction of 1.5 mm, Fig. 3.36 b). As noted before, see subsec- 
tion 3.4.5, the high axial resolution leads to the ability to capture nanometer-sized 
roughness data on all height steps while maintaining a large axial measurement 
range of 79.91 um. Furthermore, steps of (971.26 +0.31), (4951.40 +0.28) and 
(19924.00 +0.36)nm were measured over an area of 1500x250 um? without the 
need for stitching, which distinguishes the approach clearly from other techniques 
such as confocal microscopy. In the current optical design, which utilizes an imag- 
ing system with a NA of 0.06, a large lateral measurement range could be covered 
while it was not possible to gather data on the slopes with an 70° angle reliably. 
This is due to the comparatively low lateral resolution of 5 um. Depending on the 
application, the setup can be optimized to increase the sensitivity on these parts of 
the sample. The results of this sample also highlight the capability of the approach 
to decouple the axial resolution from the lateral measurement range as nm features 
can be detected while measuring over a range of 1.5 mm. 

For comparison, the height standard was also analyzed using a confocal microscope 
(Smartproof 5, Carl Zeiss Microscopy GmbH, Göttingen, Germany), Fig. 3.37 a). It 
has to be noted that the capturing of a profile having the same length as the DE- 
LCI measurement relied on stitching of multiple images, as a magnification of 20x 
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Figure 3.37 Plot of the measurements on a PTB-traceable height standard using a confocal 
microscope with a) three-dimensional representation of the sample, b) profile plot in the 
middle of the data set with inset highlighting a typical stitching error with the segments (I) 
and (II) having different distributions of noise n and Gaussian fits 


with an lateral field of view of 450x450 um? was necessary to achieve a compara- 
ble resolution. This approach is not only time consuming (approx. 10 minutes per 
areal image) but also tends to be prone to errors as stitching inconsistencies exist, 
see inset Fig. 3.37 b). It can be seen that these errors on the nm-scale influence the 
representation of the surface topography in the stitched regions which in turn can 
have an effect on quantitative analysis. The roughness distribution in these areas is 
no longer a Gaussian one, which makes it unusable e.g. for roughness evaluation, 
Fig. 3.37 c) and d). Furthermore, the RMS value of the profile data was significantly 
larger compared to the other methods. This fact is usually addressed during rough- 
ness evaluation by filtering the signal with an appropriate low-pass filter, known as 
micro-roughness filtering, [270]. 

Further three-dimensional evaluation was performed by analyzing a commercially 
available echelle grating in the Littrow configuration, Fig.3.38a). A total of five 
facets of the grating could be imaged in a lateral range of 100x200 um?, having 
a mean height of (9.66 +0.40) um. The data on the edges of the steps is notably 
noisy. A SEM scan of the grating was performed to examine individual steps as 
a reference, Fig.3.38b). From this image it can be seen that each edge has a very 
coarse structure in the size of about 6 um. These lead to very low SNRs during the 
DE-LCI measurements which are the reason that a larger area of about 6.75 um is 
obstructed on each side while only about 5.5 um of the plateaus are visible. Apart 
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Figure 3.38 Result of the measured echelle grating with a) DE-LCI measurement in Littrow 
configuration with five consecutive steps with a mean height of (9.66 + 0.4) um and an inset 
which shows the nm-fine structure of the first step at the position y=100 um and b) SEM 
image that captures the coarse structure of the individual step edges which are the reason for 
the noise in DE-LCI 


from this, the center of the plateaus could be resolved clearly with sub-nm surface 
structures, see inset Fig. 3.38 a). 


3.5.2 Alternative Spectral Encoding for Areal Measurements 


In order to gather full areal surface profile data, hence three-dimensional informa- 
tion, without any need for scanning two alternative approaches have been developed. 


Multi-slit approach 

In this approach, the measurement spot is spatially expanded and the single slit is 
substituted with a set of parallel slits in order to make use of a large area of the grating 
in the imaging spectrometer, Fig. 3.39. This arrangement allows the decomposition 
of the measured spot both spatially and spectrally. In consequence, several spatial 
parts of the measurement spot can be analyzed in the same fashion as described in 
Section 3.4 while the individual spectral slices are stacked in the y-dimension of the 
spectrometer, Fig. 3.39 b). 

In a practical realization, the multi-slit approach can be implemented by using 
imaging fibers in a linear arrangement in order to simplify the setup and to avoid 
diffraction effects from tight fitted slits. The approach holds the potential for mea- 
surements with large axial resolution while the lateral resolution of both dimensions 
is dependent on the physical size of the camera used. 
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Figure 3.39 a) Schematic representation of a modified imaging spectrometer for three- 
dimensional encoding with SPT—measurement spot which is spatially segmented by an 
LA1—lens array onto the SLI—multi-slit arrangement. A second LA2—lens array images 
the slits onto the grating which spectrally decomposes the light of every facette of LA2 while 
individual LE, —imaging elements image the components on the CAM—camera as well as 
b) simulation of a signal on the CAM where the axial information on the profile height z is 
spectrally encoded in the cameras x-dimension, individual spectral slices S, (à) are encoding 
information of one lateral dimension x, while the combination of these slices hold information 
of the second lateral dimension y 


Spatial combiner approach 

In a further approach, the encoding of a second lateral dimension is performed in the 
spectral domain. For this purpose, the measurement spot is composed of different 
spatial components which inhibit individual, discrete spectral ranges. These spectral 
slices, which are formed with low coherent light sources, are used to illuminate 
discrete regions of the sample, Fig.3.40. The data evaluation in this approach is 
similar to the conventional, two-dimensional analysis. 

The x-dimension of the camera acquires spectral information while the y- 
dimension stores information of one lateral dimension. However, the data of certain 
spectral slices with ranges of AX; to AA; corresponds to the height of the respective 
lateral dimension x,„, so that each spectral slice has to be analyzed separately. The 
resolution in the axial as well as in the lateral domain are controlled by the size 
of the spectral slices, the detector size as well as the imaging magnification of the 
measurement spot. 
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Figure 3.40 Schema of a modified imaging spectrometer for three-dimensional encoding 
with LS,—low-coherent light sources which are coupled to the interferometer by a spatial 
combiner where a BS—beamsplitter delivers appropriate beams to the SMP—sample and 
REF—reference mirror where it is also manipulated by the DE—dispersive element. Finally, 
the SPT—measurement spot which consists of spatially separated spectral ranges A, which 
are imaged onto and analyzed by the IMSPEC—imaging spectrometer 
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As discussed in section 2.2, the characterization of cross-linking and especially 
its spatial distribution is crucial for the fabrication of micro optics, MEMS, and 
semiconductors. Some studies have found that the measurement of the refractive 
index over the cross-linking process can be used as an indicator for the degree 
of cross-linking of a sample. According to Kudo et al., [130], the cross-linking 
of a polymeric material leads to a densification which can be directly related to an 
increase in refractive index using the Lorentz-Lorenz equation. The characterization 
of the refractive index is therefore a suitable measure for the degree of cross-linking, 
[139]. 

Based on the spectral interferometric approach utilized for surface profilometry 
in section 3.2 of this work, a characterization method was developed, tested and eval- 
uated. It is based on the fact that the wavelength-dependent refractive index nsmp (A) 
of a sample in correspondence with its thickness t,,, determines the spectral output 
of a low-coherence interferometer, [278, 279]. In contrast to dispersion-encoded 
profilometry, in cross-linking characterization, the sample itself is the dispersive 
element. Hence, the dispersion characteristic of n(A) is the unknown quantity. 

In the most basic configuration a two-beam interferometer with spectral detection 
is used to analyze a sample, Fig. 4.1 a). In this case, the signal at the spectrometer 
IQ.) can be described as an adaptation of Eq. (3.20), 
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IQ) = (à) - [1+cosp(A)] (4.1) 
[nP (A) Ze 1] tsmp — ô 


with g = 27 1 


; (4.2) 


where J (A) is the spectral profile of the light source and ¢ the phase. In the assumed 
simple case, the thickness of the sample tsmp is a constant whereas the path difference 
6 can be altered with a translation stage. 


4.1 Temporal Approach 


A temporal approach of the aforementioned two-beam interferometer was realized in 
a MICHELSON configuration where the sample is a transmissive part of one interfero- 
metric arm while the reference arm mirror is translatable to achieve temporal control, 
Fig.4.1 a). According to Eq. (4.1), the signal depends on the wavelength-dependent 
refractive index n°”? (A), the sample thickness fs, and the path difference between 
the arms ô. As n°"? (à) and tsmp are material constants for a supposed bulk material, 
the path difference is the only variable which can be used to evaluate the refractive 
index. A variation of the path difference leads to a deformation of the phase and 
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Figure 4.1 a) Simple setup for the characterization of transmissive samples where a WLS— 
white-light source is splitted by a BS—beam splitter so that in one arm the light trans- 
mits through the SMP—sample with the thickness f;,,, before and after it is reflected from 
SMP-M—sample mirror while in the second arm the light is reflected from REF-M—teference 
mirror which can be adjusted in the x-dimension; both signals are analyzed after recombina- 
tion using a SPEC—spectrometer. b) Simulated, spectrally-resolved phase data from such a 
setup with a sample of N-BK7 with a thickness of 5 mm for a number of different delays of 
the reference arm noted in ö 
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most notably to a shift of the phase minimum according to Eq. (4.1), Fig. 4.1 b). 
The minimum is described by its wavelength, the equalization wavelength Aeg and 
can be tracked as a function of the path difference 5(A) or the temporal delay t (À). 
This information can be used to calculate the group refractive index of the material 
using 


arg OL Ne (4.3) 


tsmp tsmp 
The measurement can be performed either in a relative or absolute way. For relative 
measurements, the delay introduced relative to a starting position (noted with 29) 
is used to calculate the relative group refractive index Ang ” (à). Absolute values 
of ng” (à) can be obtained, if the delay is referenced to the stationary phase point 
of the interferometer in a dispersion-free status. 
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Figure 4.2 a) Measured temporal delays for N-BK7 samples of 1, 3 and 5mm nominal 
thickness respectively in the spectral range from (400-1000)nm and b) calculated group 
refractive indices for the three nominal thicknesses in relation to the values known from 
literature according to [280] 


In an initial experiment, a reference mirror was placed onto a precision stage 
which was then used to introduce defined delays to the signal in form of path differ- 
ences ôn. By tuning the delay to reach a certain equalization wavelength gg, repeated 
measurements to calculate ng” (à) were possible. In order to quantify the method, 
measurements on a set of samples of N-BK7 glass with nominal thicknesses trom of 
1,3 and 5 mm were performed, Fig. 4.2 a). In an experiment with ten measurements 
per nominal thickness it was found that the mean standard deviation over all equal- 


ization wavelengths for trnom =1 mm was o1(T)=1.86 x 10-3 ps, for thom =3 mm 
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was 03(T)=1.99 x 10-3 ps and for trom =5 mm was o5(T) = 3.75 x 10-3 ps. The 
mean values of all three delay slopes were used to calculate the group refractive 
index according to Eq. (4.3), Fig. 4.2 b). The errors of the measured data have been 
calculated for every equalization wavelength relative to the respective literature val- 
ues for N-BK7, [280]. The root-mean-square error for the calculated mean group 
refractive index of the nominal thickness t,o, =1 mm was Ang; = 2.08 x 1074, 
Ang3 = 6.54 x 1075 for trom =3 mm and Angs = 1.45 x 1074 for trom =5 mm. 
The mean standard deviation over all equalization wavelengths was o1(ng) = 
5.72 x 1074, 03 (ng) = 1.93 x 1074 and 05(ng) =2.13 x 1074 for the three nominal 
sample thicknesses. As the delay was acquired as primary information, it was used 
in conjunction with the thickness of the sample in order to calculate the refractive 
index for each equalization wavelength. This calculation is not reliant on any knowl- 
edge about the underlying material model. In case of a known material composition, 
a model can be chosen in order to fit the measured data. For the evaluation data of 
N-BK7, a fit using the basic Sellmeier equation was performed using the parameters 
Bi, Bo, B3, Gis Ca and C3, 


2) 14 Bir? ‘i Br? B3A2 
n = ý 
22 — Ci A2—Cy 2 — C3 


(4.4) 


It was found that the root-mean-square error of the fitted data in relation to 
literature data for the refractive index was calculated An . 4 =2.71x 1074, An n = 
2.06 x 1074 and Anf% = 2.30 x 1074 for the three nominal sample thicknesses. 
These results prove the accuracy of the method to determine the refractive index 
of transmissive samples which in turn can be utilized to characterize cross-linking 
of polymers as typical cross-linking differences in the range of An =0.001 — 0.02 
are expected, [138]. An error propagation for the temporal approach was performed 
and documented in section 4.3.1. 

For further evaluation, two epoxy-based samples in different states of cross-linking 
have been investigated with the described method. Both samples were distinc- 
tively different in their appearance as one sample was made of bulk material 
(Araldite epoxy, eo. = 3.04 mm) while the second one was a thin-layered SU- 
8 EB =0.23mm). The properties of both samples have been studied using the 
setup described in Fig. 4.1 a) over a spectral range from (0.4-1) pm. The dispersion 
related temporal delay was recorded and referenced to the sample thickness for 
comparison, Fig. 4.3 a). The materials show significant differences in their refrac- 
tion behavior over the spectral range. While Araldite shows a delay ranging from 
1000 ps/m (@900 nm) to over 1200 ps/m (@450 nm), the SU-8 sample just causes 


delays of 850 ps/m to 1050 ps/m at the same spectral points. Both data sets could be 
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Figure 4.3 Results of the temporal approach measurements for Araldite and SU-8 samples 
with a) recorded delays due to dispersion and corresponding CAUCHY fitted curves as well as 
b) calculated dispersion D from fitted curves 


well fitted with a CAUCHY model to approximate their behavior using the parameters 
Al, A2 and A3, [281] 

Az 
24 
The difference in the relative temporal delay is an indicator of the different degrees 
of cross-linking. The first derivative of the temporal delay t in relation to the wave- 
length A was used to calculate the corresponding dispersion D = D(A), Fig. 4.3 
b), 


A 
tA) = Ait + (4.5) 


DQ = 


dT(A) 1 (E =) l (4.6) 


dÀ- tsmp tsmp 3 25 


This representation of the measured optical properties focuses on the slope of the 
dispersion induced delay. The differences between both materials are still evident 
but much smaller than the relative delay. Especially for shorter wavelengths, the 
differences are significant which indicates a different absorption behavior, hence a 
different molecular composition. This results indicates that the two materials are 
similar in their composition, but show differences in cross-linking. In the spectral 
region below 600 nm a higher dispersion of SU-8 is visible. A reason for that might 
be the different mechanism of the cross-linking initiation. While the polymerization 
of SU-8 is initiated by UV-light radiation, the polymerization of Araldite is initiated 
chemically. In the wavelength range from (600 — 750)nm the dispersion of both 
materials is approximately the same. In the near infrared region, Araldite shows 
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Table 4.1 Conclusive representation of key properties measured for both polymeric samples 
with the CAUCHY coefficients, the micro hardness and the sample thickness 


Araldite SU-8 
CAUCHY coefficients 
Ail? ] 956.9 + 3.1 793.9 + 2.7 
Al BH] 9.9+1.0 15.3 +0.7 
A3 [2E ] 0.70.1 0.10.1 
hardness [H V] 1.40 + 0.01 1.57 + 0.07 
thickness 1,7,» [mm] 3.04 + 0.01 0.23 + 0.02 


a slightly higher dispersion. Both curves are significantly different in terms of the 
measured standard deviation. 

In order to support the observations made in temporal low-coherence interfer- 
ometry, an additional material characterization were undertaken by VICKERS micro 
hardness tests. The results were compared to the obtained optical properties, Tab. 
4.1. Hardness measurements have been performed on a micro indenter device oper- 
ating with a force of 0.245 N in 5 indentations per sample. The micro hardness tests 
support the measurements of the optical properties as only minor differences are to 
be noted. This corresponds with the fact that both materials have a similar chemical 
composition. The slight differences therefore have to be caused in a difference of 
the cross-linking. SU-8 on the one hand is harder than Araldite but on the other 
hand constantly shows lower refraction over the recorded spectral range. An effect 
similar to that has been shown in literature, where the refractive index of SU-8 was 
determined at different baking steps, [282]. Here, the samples which were subjected 
to longer baking have shown higher cross-linkage as well as lower overall refraction. 

However, it should be noted that the experimental approaches for the mechanical 
and optical properties are fundamentally different. Since micro hardness is measured 
at distinct points on the samples’ surfaces, it can be heavily influenced by residual 
stresses and local inhomogenities. In contrast to that, the dispersion measurements 
integrate the properties spatially over the sample’s thickness as well as over the 
cross-section (dependent on the spot size of the light setup). Therefore, the hard- 
ness measurements are only used as a tool for the classification of the determined 
optical properties. With the aid of these measurements, it was possible to determine 
the existence of differences in cross-linkage; these should be observable in optical 
measurements as well. 
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While the presented results were gathered from samples of slightly different 
materials, further experiments were carried out that only examined one material 
with defined degrees of cross-linking. A sample set of ethylene-vinyl acetate (EVA) 
was used for this purpose. The samples where prepared from sheet material. The 
sheets were exposed to a temperature of 150°C in an industrial laminator. The 
exposure time determined the state of cross-linking. Beforehand, the samples have 
been analyzed with state-of-the-art methods by Hirschl et al., [124]. 

Especially in the field of photo voltaics (PV), EVA, which are used as encap- 
sulants there, have to maintain their properties over an operation time of 20-30 
years, [283]. Mainly, these encapsulants serve as a protection to prevent damage 
from mechanical, electrical and humid sources. They have to provide high strain 
and temperature stability to compensate for the different thermal expansion coef- 
ficients. Besides, they have to compensate for stresses and prevent cracks of the 
substrate materials. Another important function is the optical coupling of the light 
in the desired wavelength region. That demands a transmission of >90% with tol- 
erated losses of maximal 5% in 20 years, [284]. EVA is a random co-polymer of 
ethylene and vinyl acetate with a percentage of vinyl acetate typically in the range 
from 28 to 33 weight-% for PV module applications. The native EVA would not 
fulfill the thermo-mechanical requirements due to its melting range between 60 and 
70°C. By chemical cross-linking utilizing hydroperoxides during PV module lam- 
ination, the moldable EVA sheet is transformed into a highly transparent elastomer 
with the required thermo-mechanical stability up to 100°C, [124, 284]. It shows 
good adhesion, high transmission in the interesting wavelength region and it is suf- 
ficiently long-term stable regarding its properties. In order to establish the desired 
properties of polymers like EVA, it is necessary to develop and control appropriate 
curing processes. 

The prepared EVA samples of different cross-linking states were cut from sheets 
in proper pieces with a thickness of 400 um. The temporal delay t of each sample 
in relation to the white-light point was recorded in repeated measurements with 
10 repetitions each. The data was then normalized to the corresponding sample 
thickness fs, and plotted as an averaged curve in relation to the wavelength, Fig. 
4.4 a). The data points were also fitted using the CAUCHY equation Eq. (4.5). From 
the plot it becomes obvious that the differences in the delay due to dispersion are 
<10Ps/m. It is also obvious that the differences between an un-laminated sample 
(NLT) and a laminated sample (LT 1 min) are rather high. That must be considered 
especially in relation to a longer lamination period from LT | min to LT 8 min. This 
fact leads to the assumption that the curing reaction starts fast. 

In order to gather additional information, the dispersion parameter D was cal- 
culated as first derivative of the fitted data according to Eq. (4.6), Fig. 4.4 b). The 
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Figure 4.4 a) Results for the measured temporal delays t due to the different lamination 
times normalized to the material thickness fs, for three representative lamination times 
with NLT—no lamination, LT 1 min—1 minute and LT 8 min—8 minutes lamination and 
the corresponding fits using a CAUCHY model according to Eq.(4.5) and b) derivation of the 
temporal delay relative to the wavelength according to Eq. (4.6) 


results reveal that there is no particular difference between the varying degrees of 
cross-linking. Although some differences in the wavelength range of (0.4-0.6) um 
can be observed, the corresponding errorbar proves that the dispersion slope in rela- 
tion to the wavelength is constant for different degrees of cross-linking. This leads 
to the assumption that the magnitude of the temporal delay and therefore also of 
the group refractive index can be used as a measure for cross-linking differences. 
In contrast, the slope of the wavelength-dependent curves is not a suitable measure 
to make out cross-linking differences. 

The most important information, the degree of cross-linking, can be extracted by 
plotting the differences in temporal delay At/tsmp versus the lamination time, Fig. 
4.5. The data acquisition over a broad spectral range enables the analysis at different 
wavelengths. For comparison, the relative delay at three probing wavelengths (0.5, 
0.7 and 0.9) um was analyzed. It is visible that the data is equal within the standard 
deviation of the measurements especially for the spectral probing points of (0.7 and 
0.9) um. The data at 0.5 um shows a slight deviation. In order to compare the results 
to other methods, the measurements have been fitted using a pseudo-first order 


reaction kinetics model! of the form y = a + b - e~* where k is the characteristic 


! In situations where the amount of cross-linker is small compared to the amount of polymer, 
the reaction can be described as a first order reaction. In this context, the term pseudo-first 
order is used, [124]. 
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Figure 4.5 Plot of the differences in At/tsmp) over the lamination time selected probing 
wavelengths 


constant of the reaction process and a, b describe the absolute position of the slope. 
This model is commonly used to describe cross-linking behavior determined by 
characterization approaches, [124]. As good fits (R? = 0.976, 0.985 and 0.989) could 
be obtained at all probed wavelengths, an averaging over a range of wavelengths can 
be used to determine the degree of cross-linking with higher statistical confidence. 
In comparison to other technologies, the presented results show reasonable errors in 
the range of 6.35-8.39%. Also, the trend of the data is in good analogy to reference 
technologies like soxhlet-extraction, Raman spectroscopy, DSC or DMA, [124]. 
EVA samples show fast cross-linking at the beginning of the lamination process 
which significantly slows down after 3-4 minutes. The calibration of the At/tsmp 
data to a degree of cross-linking on a percentage scale can be done by choosing an 
appropriate reference technology. 

In consequence, the temporal approach showed the ability to gather model-free 
refractive index data over a large spectral range, (0.4-1) um, with deviations to lit- 
erature values in the range of 2.06 x 1074 - 2.71 x 1074. It could be shown that 
cross-linking of polymers for industrial applications can be evaluated on the basis 
of the refractive index measurements. The results proofed to be comparable, and 
in terms of their standard deviation, more reliable than the results of established 
technologies such as soxhlet-extraction or Raman spectroscopy. The additional sta- 
tistical confidence through the measurement over a large spectral range as well as 
the ability to gather spectrally-resolved refractive index data are advantages of this 
approach over established methods. 
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4.2 Scan-free Approach 
4.2.1 Wrapped-phase Derivative Evaluation (WPDE) 


One significant drawback of the temporal approach to estimate the degree of cross- 
linking with the refractive index is the need for mechanical scanning of one interfer- 
ometric arm. It possibly introduces additional errors and increases the measurement 
time. 

From Eq. (4.2) it is known that the wavelength-dependent refractive index 
n“”"P (À) is contained in the phase of the interferometer output in combination with 
the sample thickness tsmp and the path difference ô. Under the assumption that fsmp 
as well as 6 are known, the relevant cross-linking information can be found in the 
refractive index. 

By rewriting Eq. (4.2) the measured phase-term Ømeas, containing the refractive 
index, can be extracted 


7 + poff- (4.7) 


Pmeas = cos! (= 2 E ı) = 20 aA) = Dismp Z8) 
Inherent to this approach is the ambiguity of the resulting values as ¢ is not limited to 
the range of 0 - x. Other works have shown methods to perform the correct quadrant 
selection in order to resolve this ambiguity, [285]. In contrast, an alternative method 
to avoid quadrant selection was developed by performing a local signal analysis in 
the spectral range close to Aeq, Fig. 4.6 a). In the first stage, this approach determines 
the phase minimum and defines a ROI around the minimum. For this purpose, the 
raw measured data is analyzed using a STFT where a FFT is performed in one small 
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Figure 4.6 a) Simulated data of a wrapped phase due to the cos”! operation according to 
Eq. (4.12) with marked equalization wavelength Aeg and ROI for the extraction of loc 
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window of the complete data set which is then slid over the signal successively along 
the wavelength dimension, see subsection 3.3.2. This approach accounts for the non- 
uniform frequency of the signal. As a result, the minimum of the extracted frequency 
slope can be determined from the power spectrum. It represents the position of the 
phase minimum which also occurs at Aeq. The ROI is defined as a local wavelength 
range Ajoc in the proximity of the detected Aeq where only unambiguous phase data 
is included. This so-called local phase, groc, is subject to a phase offset, poff, with 
regard to the absolute phase due to the cos~!-operation, Eq. (4.12). 

A second analytical step implements a newly developed approach called WPDE, 


where Ploc is differentiated with respect to the wavelength, noted with a. 


IPloc u ð 27 [n(Atoc) — 1]tsmp — d(x) 
OA ~ oA Aloc 


+ vats) : (4.8) 


This eliminates the phase offset Yoff and enables the evaluation of the cross-linking 
characteristics in terms of the group refractive index ng” (x, A) as well as the relative 
derived optical thickness (RDOT) to PT 


ig OA) =1-— (4.9) 
smp 
1 -A2 
withk = Ploc (4.10) 
27 
Ploc A? 
topr =ng P (x, A) - tsmp = tsmp = ö (4.11) 


where gj, is calculated from the measured data using the difference quotient with 
AA as interval. This case holds true when experiments, as sketched out in Fig. 4.1 a), 
are performed where one simple sample is part of the interferometer as well as the 
primary source of dispersion. In situations where samples with low dispersion are to 
be measured or the simultaneous measurement of the samples surface profile should 
be realized, a modified setup with additional dispersion is favorable, Fig. 4.7. In this 
case, the phase term of Eq. (4.7) has to be expanded by an appropriate term for the 
dispersive element, 


sate) (4.12) 
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According to this equation, the derivative in order to access the group refractive 
index can be noted as 


Ploc z Ir À + Poff ? 


(4.13) 
which leads to a new description of the group refractive index and the RDOT 
to pr for the approach with additional dispersion 


ð (z [n (x, A) = 1) timp | + [(n?P F(a) = 1) toe | —6 


WE 
z @,A) = 1- > 4.14 
Ng (x, A) a smp ( ) 
; 27 
with = pl. — 55 I(1-nP#) toe +5] (4.15) 
a2 
topr = ng Gs A) . tsmp = tsmp = = (4.16) 


A detailed derivation of Eq. (4.9) and (4.11) for the sample-only approach as well 
as for the approach with additional dispersion resulting in Eq. (4.14) and (4.16) can 
be found in the appendix in the Electronic Supplementary Material (ESM). 
In order to evaluate the algorithm, the group refractive index of a N-BK7 sample 
with a nominal thickness of 5 mm was determined. The averaged standard deviation 
of 10 consecutive measurements of the sample was found to be 9.97 x 1075. The 
averaged group refractive index data of these 10 measurements was fitted using a 
SELLMEIER equation, Fig. 4.7 b). In resemblance to the literature values, [280], a 
root-mean-square error of 1.65 x 1074 and of 3.36 x 1075 was achieved for the 
averaged measured and for the fitted data respectively. Compared to the measure- 
ments using the temporal approach, subsection 4.1, this demonstrates an improve- 
ment as the RMS error was An i =2.30 x 1074. An additional advantage over the 
temporal approach is the ability to gather the wavelength-dependent group refractive 
index without the need for mechanical scanning. The result shows that the WPDE 
approach achieves a comparable accuracy to state-of-the-art refractive index mea- 
surement technologies. Furthermore, the refractive index resolution is sufficient to 
characterize cross-linking in waveguide polymers, where differences in the range 
of An = 0.001 - 0.02 are expected, taking the respective sample thickness into 
account, [138]. 

This result is calculated only within the ROI and is dependent on the amount of 
dispersion, represented by ne” (à). Therefore, it is valid only within a small spectral 
range. Different approaches have been considered to gather information over the 
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Figure 4.7 Setup for the imaging approach to cross-linking characterization with WPDE 
incorporating a WLS- white light source which is split into a reference arm with a REF-M 
reference mirror and a sample arm which consists of a SMP—sample of fsmp on a 
SMP-M—sample mirror and an optional DE—dispersive element of tpg. The beams of both 
arms are recombined by the BS—beamsplitter and imaged by a LE—lens onto the IMSPEC— 
imaging spectrometer. LE can be translatable in the y-dimension to gather areal cross-linking 
information and b) Plot of the averaged measured group refractive index of N-BK7 (t = 5mm, 
10 measurements) which was calculated using the WPDE approach and its corresponding 
SELLMEIER fit in comparison to the literature values according to [280] 


complete spectral range of the data set. On the one hand, the WPDE analysis algo- 
rithm can be applied to other ROIs within the data. The advantage is that the group 
refractive index can be calculated without an a priori knowledge of the underlying 
material model. On the other hand, one can calculate the group refractive index over 
the complete spectral range, if the material model of the sample is known. 


4.2.2 Spatially-resolved Approaches 


All approaches described so far have been based on point-wise measurements 
of spectra and relied on scanning either one arm of the interferometer to gather 
wavelength-dependent information (temporal approach) or on scanning the sample 
in order to gather information from different locations of the sample (WPDE). 


Scan-free temporal approach 

The evaluation of the refractive index of a sample and therefore of the degree of cross- 
linking relied on mechanical scanning of one interferometer arm in order to scan the 
spectral domain for certain equalization wavelengths eq. The underlying principle 
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Figure 4.8 a) Proposed setup for the scan-free temporal evaluation of samples with a WLS— 
white-light source split by a BS—beam splitter into one arm that focuses light into a SMP— 
sample which then gets back reflected by a mirror and a second arm, known as REF—reference 
arm where light is focused in the same way before it gets back reflected on a mirror which 
introduces controlled temporal delays in the x-domain 6, and a IMSPEC—imaging spectrom- 
eter in order to analyze the signal in a spatially-resolved manner and b) a simulated signal of 
the imaging spectrometer with different temporal delays ô and corresponding equalization 
wavelengths Aggn 


is the introduction of different temporal delays to the setup. A possible method 
to introduce the delays all at once, relies on the spatially encoding of them and 
the appropriate detection, Fig. 4.8. For this purpose, the standard one-dimensional 
spectrometer is replaced by an imaging spectrometer with appropriate optics. The 
imaging setup enables the detection of spectral information in one spatial domain. 
In the proposed setup, the probing beam is focused in the sample volume and re- 
collimated onto a mirror. Consequently, the beam of the reference arm is also focused 
and re-collimated without a sample being present. The mirror of the reference arm 
is designed to introduce temporal delays to the beam with a spatial distribution. 
After recombination of both beams and their spectral detection, this distribution of 
delays can be recorded as spectra with different adequate equalization wavelengths. 
Analogous to section 4.1, these wavelengths can be used to calculate the refractive 
index of the sample at these discrete points. Depending on the number and size of 
the delays as well as the construction of the spectrometer, a large spectral range 
can be covered. Model-based fits can be calculated accordingly. The introduction of 
delays can be done by means of a stepped mirror, a deformable micro-mirror array 
or a transmissive element with a refractive index gradient. 
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Imaging WPDE 

In a second approach, the aforementioned implementation of an imaging spectrome- 
ter was also used to perform WPDE with spatial resolution using the setup described 
in Fig. 4.7 a). Within the setup, a sample with the thickness tsmp is placed on a reflect- 
ing substrate and is used as a mirror for one interferometer arm. Correspondingly, 
the reference arm only compromises a mirror and no additional dispersive element. 
Depending on the thickness and amount of dispersion of the sample, an element 
with additional dispersion with the thickness tpg might be necessary in the sample 
arm in order to enhance the measurability. A detailed explanation on the usage of an 
additional dispersive element is given in subsection 4.3. The setup is also equipped 
with a translation stage for the imaging lens LE. This lens enables the recording 
of areal cross-linking information. As the imaging spectrometer allows capturing 
refractive index data along a line, the translation of the imaging lens in the y- 
dimension enables the stacking of these line profiles in order to receive information 
on the whole two-dimensional plane, refer also to Fig.3.36a). Using this method 
neither the sample nor the reference arm have to be moved during measurements 
which prevents obstructions of the interferometric measurement due to movement. 
The data analysis was performed analogously to the WPDE approach described in 
section 4.2.1. As already pointed out, in case of the usage of an additional dispersive 
element, a modified set of equations, Eq. (4.14), has to be applied. 

This approach was used to characterize lithographically generated structures in 
a photo-resist, Fig. 4.9. The resist was spin-coated on a Si-wafer with a thickness 
Of tsmp =750 um. Afterwards, it was exposed to visible light, (400-420) nm, for 
primary cross-linking and to UV-radiation, (300-360) nm, in a secondary cross- 
linking process. The secondary cross-linking was performed through a mask to 
generate structures of rectangular refractive index patterns which also lead to the 
shrinkage of the cross-linked areas. The goal of the investigations was to determine 
the surface height profile that is altered due to shrinkage as well as the refractive 
index profile which is due to different degrees of cross-linking. 

The surface height profile of the sample was characterized with the described 
setup using the front-surface reflex and the profilometry approach described in 
chapter 3, Fig. 4.9 b). It is obvious that apart from a slight overall waviness, the 
sample shows a regular height pattern with the expected pitch of 50 um. The depth of 
the shrunken areas is about 120 nm, which lies in the expected range. In consequence, 
these calculated height profiles enable the separation of shrinkage from the refractive 
index information for every sample individually and simultaneously. 

With the knowledge of the surface height profile of the sample due to shrinkage, 
the correct thickness along the spatial domain, tsmp = tsmp(x) ~ z(x), can be 
calculated. Therefore, the surface height profile was measured in relation to the 
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Figure 4.9 a) Detail of the photo-resist sample under test with a lithographically gener- 
ated, spatially-dependent refractive index pattern and a resulting surface height profile due to 
shrinkage. After transmission through a DE—dispersive element with tpg a FSR—front-side 
reflex from the sample can be used to interfere with light from the REF-arm—reference arm 
in order to calculate the surface height profile z(x) ~ tsmp(x) while a BSR—back-side reflex 
can be used in conjunction with fsmp (x) to calculate the refractive index profile n(x) as a mea- 
sure for the degree of cross-linking across the sample b) plot of the measured surface profile 
z(x) from a polymer sample under investigation utilizing a wavelength-calibrated imaging 
spectrometer 


substrate. By the application of either Eq. (4.14) or (4.16), the group refractive 
index or the relative derived optical thickness can be calculated corresponding to 
its position on the sample, Fig. 4.10 a). For the results pictured above, the RDOT 
profile of the sample was calculated for a single wavelength of 557 nm. The spatial 
profile allows a resolution of cross-linking differences of 4 um in the lateral domain. 
Although the results are affected by noise and batwing-effects, [263], a dynamic 
range of +1.5 um in the RDOT for the given sample was revealed over a lateral 
range of nearly 550 um, while a section of 250 um is displayed here. Furthermore, 
it also has to be noted that the plateaus do not show completely flat RDOT profiles. 
This behavior was attributed to a mixture of effects ranging from diffraction during 
exposure of the structures to deformation during shrinkage and diffraction during 
measurements. As the profile was taken at a specific wavelength, it represents only a 
fraction of the captured information, which was originally analyzed over a spectral 
range of 20nm. 

In order to estimate the effect of cross-linking, the RDOT differences have been 
measured over the complete spectral range as a mean value of two different exposed 
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Figure 4.10 Results of the measured RDOT a) spatially resolved along one sample dimen- 
sion of a lithographically structured photo-resist layer with structures having a nominal pitch 
of 50 um on a Si-substrate at a wavelength of 557 nm and b) mean values and fitted data for 
two marked areas with different degrees of cross-linking over a spectral range 


areas, Fig. 4.10 b). An RDOT difference of about 3 um between the differently cross- 
linked areas could be resolved while the RDOT slope for every area was determined 
over 10nm. The results are affected by noise in the original data which is amplified 
by the process of taking the derivative. Some smoothing with a Gaussian filter was 
applied to the data. 

One of the main advantages of the described approach is the lack of neces- 
sity for a model in order to calculate the spectrally-resolved refractive index. As 
some compromise towards the size of the spectral measurement range was made by 
the choice of the dispersive element, subsection 4.2.1, the application of a refrac- 
tive index model might become interesting in post-processing. In the context of 
(photo-)polymers, a variant of CAUCHY’s equation was selected, [281, 286]. Using 
this model, the group refractive index ng” (à) can be calculated according to Del- 
barre et al. [246] with 


sm dn(A) 3A 5A 
ne”? (A) = n(a) rm Ait en 2. 


(4.17) 


By appropriate fitting, the CAUCHY coefficients A1, A2 and A3 were determined 
which enable the calculation of the refractive index and the group refractive index 
over any given spectral range where the CAUCHY model is valid. 
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4.3 Influences and Limitations 


The measurement range as well as the accuracy of the cross-linking determination is 
dependent on the accuracy of the determination of the equalization wavelength Aeq, 
the sample thickness f,,,, and the path difference of the interferometer ô amongst 
other parameters, Eq. (4.1) and Eq. (4.7). While the position of the equalization 
wavelength depends on the path difference between both interferometer arms, the 
width of the fringe around Aeg is determined by the amount of optical dispersion. 
Hence, thicker materials show tighter fringe spacing than thinner samples of the 
same refractive index, Fig. 4.11. The different measurement approaches discussed 
in this chapter demand different signal types for analysis. While the detection of 
just the equalization wavelength with the temporal approach works with very lit- 
tle dispersion, hence a very wide fringe spacing, the WPDE method requires a 
tighter fringe spacing in order to resolve the phase minimum in the power spectrum, 
Fig.4.12. The estimation of the phase minimum is performed as tracking for the 
minimal frequency of the component carrying the most power in the spectrogram. 
An analysis of the frequency content for three distinct wavelengths, Aeq = 0.6 ym 
as well as 0.5 and 0.8 um, visualizes the problem of low dispersion, Fig.4.12b). 
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Figure 4.11 Simulated signal in order to demonstrate the influence of dispersion on the 
relative signal intensity and phase slope at an equalization wavelength of Aeq =0.6 um and 
sample of N-BK7 with a) tsmp =1 mm, b)tsmp =3 mm and c) tsmp =5 mm 
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Figure 4.12 a) STFT analysis of the simulated signals in order to demonstrate the influence 
of dispersion at an equalization wavelength of Aeq =0.6um and a sample of N-BK7 with 
(1) tsmp = 1mm, (ID) tsmp =3 mm and (III) tsmp =5 mm with colored markings to indicate the 
data used for b) visualization of the limitation to determine Aeq for setups with dispersive 
elements <1 mm where the frequency resolution of the STFT is not sufficient to resolve the 
peak for the phase minimum (here at 0.6 um) from the rest of the phase signal 


The lower the dispersion in the setups, the closer the not infinitely sharp peaks of the 
power spectrum move towards the lower frequencies. For samples of N-BK7 with 
thicknesses below 1 mm, a distinct separation is not possible anymore. Potential 
solutions to this problem can be the application of an adaptive windows size during 
the STFT which would decrease the peak width for frequency peaks as well as the 
introduction of additional dispersive elements. The exact thickness which necessi- 
tates additional dispersion is dependent on the sample’s refractive index. Also, it 
has to be noted that the increase in dispersion reduces the spectral range as only 
the data within one phase jump around A., is used in the algorithm. In order to 
compensate for the loss in resolution due to the additional dispersion, the before 
mentioned precision fit of the spectral data around the equalization wavelength can 
be utilized. 

A detailed description of significant influences and error sources is given in the 
following section. 


Influence of shrinkage 

The curing and cross-linking of polymeric materials results in a change of the 
refractive index as well as in dimensional changes like shrinkage. As pointed out in 
Eq. (4.1), both properties influence the measurable signal in low-coherence interfer- 
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Table 4.2 List of key properties of an exemplary positive photo-resist before and after 
processing in a UV-harden cross-linking step which were used during the simulation, [H. 
Aßmann\ Th. Albrecht, personal communication, 12.07. & 30.08.2018] 


property value 

resist type positive photo-resist 
tsmp before processing 1878.32nm 

tsmp after processing 1755.49 nm 

Atsmp /tsmp 6.5% 

n*P before processing (@675 nm) 1.6168 

n*P after processing (@675 nm) 1.6702 

An"? In"? 3.2% 


ometry. In order to provide a meaningful metrology tool, the influence of both prop- 
erties on the measured data was to be studied. For this purpose, simulations using 
an industrial photo-resist material system were performed. The material properties 
after different processing, e.g. softbake/hardbake and UV-hardening, were deter- 
mined with the reference techniques tactile profilometry and spectral ellipsome- 
try, [H. Aßmann\ Th. Albrecht, personal communication, 12.07. & 30.08.2018], 
Tab. 4.2. It can be seen, that significant changes in thickness as well as in refractive 
index occur in a counteracting fashion. In order to calculate the influence of these 
two properties on the measurement data, two simple simulation cases have been 
studied. 


(A) transmissive sample: A transmissive sample with a given thickness tsmp and a 
reflective index n*”? is placed in one arm of the interferometer. By changing one 
of the parameters during simulation while leaving the second fixed, the individual 
influences can be estimated, Fig. 4.13 a). 

(B) reflective sample: A defined layer of photo resist is spin-coated on a reflective 
surface (e.g. Si or glass). In this configuration, the expected measurement signal 
can be composed of the different components such as the front side reflex of 
a shrinked and non-shrinked sample as well as of the back side reflex from a 
sample with different states of cross-linking, Fig. 4.13 b). 


When analyzing case (A), a thickness change due to shrinkage of 6.5% was used 
to calculate the influence on the phase signal and equalization wavelength by also 
using a fixed refractive index of n*’”?(A) for the sample., Fig.4.14a). The result- 
ing phase signal showed a relative change of about 0.13% in its amplitude for the 
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Figure 4.13 a) Setup for the simulation of the influence of a thickness or refractive index 
change on the signal in transmission mode of case (A) with LS—light source path, SMP— 
sample with refractive index n“”"? (A) and thickness tsmp, SMP-M— sample mirror and REF— 
reference beam path; b) setup for the simulation of shrinkage and cross-linking of a polymer 
on a reflective sample of case (B) with DE—dispersive element having a refractive index 
nPE (A) and a thickness tpg, non-cross-linked sample with thickness tnc; and cross-linked 
sample with thickness tnc; as well as the respective FSR—front side reflexes which influence 
the path difference ö,„.ı and 6,; and the BSR—back side reflex which is influenced by np (A) 


smp ncl 
ong (A) 


equalization wavelength, whereas the equalization wavelength showed a change of 
about -0.038%. It can be noted, that the changes are small relative to the thickness 
change. Furthermore, the equalization wavelength change is an order of magnitude 
smaller than the change of the phase signal. In order to evaluate the influence of 
a change in refractive index An, the thickness of the simulated material was kept 
constant, while the refractive index was changed by 3.2%, Fig. 4.14 a). In this case, 
the phase change was found to be -0.18%, while the change in equalization wave- 
length was 0.074%. These results make clear that the influence of the refractive 
index change alone is stronger than a thickness-induced change. Furthermore, as 
both effects are counteracting, a separation of the changes during a measurement 
might be obscured. Although shrinkage-induced changes of thickness and refractive 
index are expected to be in the single percentage regime, [138], this effect has to be 
taken into account. This result also reveals the necessity of a more profound data 
analysis as shown in chapter 3. Especially the fit of the measured data in the region 
of the equalization wavelength, enables high precision in the measurement of both 
thickness and refractive index of a (polymeric) sample. 
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Figure 4.14 Simulated phase signals of the influence of changes in sample thickness and 
refractive index for a) study of case (A) where a sample is transmitted and both changes are 
analyzed separately and b) study of case (B) where front side (FSR) and back side reflections 
(BSR) from a non-hardened and hardened sample have been investigated 


A more realistic assessment of the signals has been performed by the study of case 
(B), Fig. 4.13 b). This case describes four relevant signal components in the analysis 
of a single layer of polymer, where in one part the sample is assumed to be non-cross- 
linked (initial) and in the other fully cross-linked (hardened). As described before, 
all signals passed through a known dispersive element. In case (B).1 the front side 
reflex (FSR) of the non-cross-linked material interface interferes with the reference 
signal while in case (B).2 the front side reflex of the cross-linked material is analyzed. 
In this case, as pure shrinkage was studied, a change in the equalization wavelength 
of 0.19% and a change in the phase signal of -0.52% could be observed. In the 
cases (B).3 and (B).4 the effect of the refractive index is studied for both samples 
while the shrinkage effect is included in the signal of the cross-linked material. By 
analyzing the back reflected signal from both, the cross-linked and the non-cross- 
linked sample, it is possible to investigate not only the effect of shrinkage but also 
the effect of refractive index variation due to cross-linking. The results show clearly 
the counteracting nature of shrinkage and refractive index alternations during cross- 
linking with a change in the equalization wavelength being 0.037% and a change 
in the phase signal being -0.036%. This highlights that a refractive index induced 
change can be nearly obscured by the influence of the shrinkage of a material. For 
this reason, all experiments within this work performed an evaluation of the FSR 
as well as of the BSR signal. This approach enables separating signal changes due 
thickness variations from those that are caused by refractive index variations and 
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can deliver information on the surface profile of the sample. During the evaluation 
of thin materials (ft < 100 um), the signals of front and back side reflection mix. In 
order to perform the correct interpretation of the measured data from cross-linked 
and non-cross-linked samples, both signal parts have been separated. This can be 
done by filtration or temporal separation. 


4.3.1 Error Parameters of the Temporal Approach 


The determination of a sample’s group refractive index, based on the temporal 
approach, relies on the translation of the reference mirror in order to capture the 
path length difference for a number of equalization wavelengths, see subsection 4.1. 
The final calculation relies on the path length difference 6, and the thickness of the 


sample tsmp 
en 
nÀ) = ( ) (4.18) 


tsmp 


In order to estimate the error for this measurement, an error propagation was per- 
formed with respect to both relevant features 


A Ome ABCA) i PITON i (4.19) 
ng = - - i i 
8 952) I P 
The calculation of the derivative of the respective terms leads therefore to 
2 2 
1 5(A) 
Ang = -AdA)) +|- Atsmp | - (4.20) 
Ismp temp 


Based on the parameters used for the experiments, the error limits were calculated 
for the dispersive elements of tsmp =(1,3and5)mm while the measurement accu- 
racy on these thicknesses was Atsmp = 20 nm when measured with a tactile pro- 
filometer, [275]. The covered path length differences were ö(%)=(49.8,185.5 and 
345) um respectively while the resolution of the translation stage was Ad(A) = 1 nm, 
[287]. This resulted in values of An} =1.41 x 1076, An3=5.30 x 1077 , An} = 


3.41 x 1077 for the different measured samples in section 4.1. 
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4.3.2 Error Propagation in WPDE 


Independently of the kind of WPDE, sample-only or with an additional disper- 
sive element, the accuracy of the refractive index calculation is dependent on the 
deviation of the measured input parameters such as the sample thickness tsmp, the 
wavelength X or the thickness of the dispersive element tpg. Therefore, a propa- 
gation of deviations of these parameters was performed to estimate their relative 
influences following the scheme 


2 2 2 
| ech Eu EL UNSER ER on (4.21) 
g EIER smp Jint ad ; ; 


For all investigations on this topic, samples are supposed to have thicknesses varying 
from (0.1-5) mm which were measured with a deviation ranging from (0.16-20) nm, 
[275]. The dispersive element was made of N-BK7 glass with the given thicknesses. 
Wavelength dependencies were investigated in a spectral range of (400-1000) nm 
which was determined with a deviation of 0.1 nm, [288]. 


Sample-only WPDE 

When using the WPDE method for simple transmissive measurements of bulk mate- 
rials, see section 4.2.1, the refractive index depends only on the deviation of the 
sample thickness At, , and the wavelength AA 


sngl 2 sngl 2 
I on on 
Ang = ( T £ Ar) + ( F k s) . (4.22) 
smp 


According to this notation the partial derivative of the group refractive index relative 
to the sample thickness 


i A2 
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reveals a simple quadratic dependency which is countered by the effective thickness 
and the relative path difference 6 that contribute to the measured phase signal 
Pioc: In the given range of sample thicknesses the deviation of Ang(tsmp) could 
be estimated at the sodium D1 line (A=589.592nm), Fig. 4.15 a). The quadratic 
dependence becomes visible. The data here was plotted on a semi-logarithmic scale; 
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it has to be noted that the error is smaller than 1 x 1075 for samples thicker than 
500 um. In a measurement situation this might be considered as a limiting factor 
for the characterization of thin samples or layers. 
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Figure 4.15 Error contribution of the sample thickness to the group refractive index 
Ang(tsmp) for the sample-only WPDE approach with a) over a given thickness range at 
the sodium D1 line (A=589.592 nm) and b) over a spectral range with tsmp =5 mm and the 
measurement error for Atsmp=20nm 


Furthermore, the contribution to the error for a sample of N-BK7 having the 
thickness fsmp=5mm was analyzed over a spectral range from (0.4-1) um, Fig. 
4.15 b). The calculated error for this sample at the sodium D1 line is 2.16 x 10°. It 
changes less than 0.3 x 1076 over the given wavelength range. In consequence, for 
high precision measurements the analysis wavelength should be as low as possible 
in order to achieve results with minor error although the gain in accuracy is small. 
Accordingly, the error contribution resulting from the wavelength measurement 
uncertainty 


anins! š ic _ 2a [a Snes 4. 5] 
ar vie l t oz Met = 
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was estimated for a N-BK7 sample having the thickness tsmp =5 mm over the spectral 
range of 0.4 to 1 ym, Fig.4.16a). Interestingly, the error contribution is quiet high 
with +5 x 1075 but becomes zero at a wavelength of 494.9 nm. This results from 
the ratio of the group refractive index to the wavelength squared which in case 
of N-BK7 results in a zero error within the wavelength range of interest. In this 
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Figure 4.16 Spectral dependency in a range from 0.4-1 um with a) error contribution of 


the wavelength to the group refractive index Ang ® ; (A) for the sample-only WPDE approach 


with tsmp =5 mm and the measurement error for AA =0.1 nm and b) total error on the calculated 
group refractive index of the sample with all discussed error contributions according to Eq. 
(4.22) 


case, experiments should be designed to measure the group refractive index with 
the WPDE method in a wavelength range close to this minimum. The error varies 
from 5.17 x 1075 at 0.4 um to 4.93 x 1075 at 1 um. In case of different materials, 
simulations have to be performed to find an optimized probing wavelength. 

Consequently, the overall error according to Eq. (4.22) shows a minimum at this 
wavelength, Fig.4.16 b). It can be seen that for a sample thickness of tsmp =5 mm 
the contribution of the wavelength is quite significant and introduces the main error. 
Furthermore, it is known from the simulations that the sample material as well as 
its thickness can have an even more significant influence on the error of the calcu- 
lated group refractive index. This behavior should be considered when designing 
experiments for arbitrary materials. 


WPDE with additional DE 
As discussed within this chapter, WPDE can be applied in some cases with an 
additional dispersive element in the setup, Fig.4.7. In this case, the propagation 


of deviations is based on a different equation for the determination of the group 
refractive index nE? E Eq. (4.12) 
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Furthermore, it additionally depends on the measurement uncertainty of the thick- 
ness of the dispersive element Atpg 
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Similar to the case in subsection 4.3.2, the deviation relative to the error of the 
sample thickness is quadratically dependent on the sample thickness itself 


antDE 3 2E 
8 
- Atsmp = —— | ———— | At, 4.28 
ðtsmp MPO Stomp \ 27 i mee Geet) 
1 a? ` Pioc [( DE 
= l-n )toe +8] - Atsmp 
oe ( 27 2 
with 2 
I 
Ploc = 77 [a-n +a -np E +] (4.29) 
anger l ff ff 
sm e e 
£ A Atsmp mer we [a —Ng P) . tsmp + ad = nPE) (Bf = toe) | i Atsmp- 
Otsmp temp 
(4.30) 
It has to be noted that the calculation of g),. using the effective thicknesses 


i, and e is performed only in context of the error propagation in order to 
simulate the signal which is usually measured. Using the same range of sample 
thicknesses, with a dispersive element of tp g = 5mm N-BK7, a thickness-dependent 
error can be calculated for the sodium D1 line (à = 589.592 nm), Fig. 4.17 a). It is 
clear from the equation and the plot that the contribution of the dispersive element 
to this particular error is minimal. The measured thickness tpg as well as the effec- 
tive thickness i , which contributes to the composition of the phase signal ¢/,_.., 
only affect the error with their difference to each other. Therefore, analogous to 
the approach with a single sample material, the error contribution with respect to 
the sample thickness dominates. It is smaller than 1 x 1075 for sample thicknesses 
above 420 um for measurements at the sodium D1 line. With regard to the spectral 
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range of the measurement, the eg behavior is smaller than 3 x 1075 and thus 
neglectable in comparison to the influence of the sample thickness, Fig. 4.17 b). 
An error variation of 1 x 1076 can be observed over a spectral range of 0.6 um. In 
consequence, the measurement wavelength range should be as low as possible in 
order to optimize the error contribution of the sample thickness. 
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Figure 4.17 Error contribution of the sample to the group refractive index Ang(tsmp) for 
the WPDE approach with additional dispersion a) over a given sample thickness range at the 
sodium D1 line (A =589.592 nm) and b) over a spectral range with tpg =5 mm, tsmp =750 um 
and the measurement error for Atsmp =4nm 


When partially deriving Eq. (4.25) with respect to the thickness of the dispersive 
element, 


ön+PE (1-nPE) 
4 AtpE = — 8 7. Atpg (4.31) 
OtpE Ismp 


itbecomes obvious that the error contribution of the dispersive element is determined 
by its refractive index and not by the thickness itself. However, the thickness of the 
sample shows a 1/tsmp influence on this error contribution, Fig. 4.18 a). The error 
contribution becomes significant for samples thinner than 200 um for the sodium 
D1 line where it is larger than 5.41 x 1075. The differences in the error contribution 
for different wavelengths are in the range of +.09 x 1075 over the spectral range 
of 0.4-1 um while its maximum value is 1.56 x 1075 for a wavelength of 0.4 um, 
Fig. 4.18 b). This leads to the same consequence as for the error contribution of the 
sample thickness alone. In order to minimize the error in the calculation of the group 
refractive index, the spectral range for measurements should be as low as possible. 
This can be achieved by the choice of the light source, the spectrometer and by 
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tuning the path length difference of both arms as to shift the ROI appropriately, see 
Fig. 4.7. 
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Figure 4.18 Error contribution of the dispersive element to the group refractive index 
Ang (tpg) for the WPDE approach with additional dispersion a) over a given sample thickness 
range at the sodium D1 line (A=589.592 nm) and b) over a spectral range with tpg =5 mm, 
tsmp = 750 um and the measurement error for Atpg =20nm 


Furthermore, the group refractive index was also partially derived with respect 
to the wavelength X in order to estimate its influence on the error. For this operation, 
the equation was rewritten and broken down into its main wavelength dependent 


components X and Y 


ant PE a ra? [Ln PE (aye ô 
& A= Ploc 8 me -AA (432) 
oA aA 27 - tsmp tsmp tsmp 
Eo A2 
where “0” x (4.33) 
T tsmp 
[1 — nPE Ot 
ma E PE y (43A 


tsmp 


In this notation the derivative of X with respect to A can be written as ax 
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while the derivative of Y with respect to A can be written as oy 


= 4.36 
oA tsmp OA eee) 


tsmp 


ay a (oe) tpg anPE 
aa u 


In the case examined here, the dispersive element is made of N-BK7 glass which 
leads to the use of the Sellmeier equation for the refractive index, [280], 


1. (4.37) 


„PE _ Ar? 4 A22 A" A3A2 
— Bı 22 — Bo A2 — B3 


This equation is 
group refractive index is computed using, [246], 


dn 
=n——-d. 4.38 
Ng =n rr ( ) 


By caleulaine the derivative of the refractive index for glass with respect to the 


wavelength “ Sz» the group refractive index for N-BK7 can be formulated as 
G 
DE DE 
= — A 4.39 
i Sn 2. DE (4:39) 
—2A; Biù 
with G = 3 0I- Bye By?" (4.40) 


In consequence the partial derivative of the group refractive index can be formulated 
as 


DE 2 
ang G (igatG).n nPE- DE A 4.41 
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Eq. (4.41) is used to eventually determine the error contribution of the wavelength 


dependency to the calculation of the nn group refractive index with Eq. (4.32), 
DE 


the solution of 2X X from Eq. (4.35) and X X from Eq. (4.36) in conjunction with —— rg 
from Eq. (4. 41) 
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It has to be noted that the calculation of g/,. is only performed in context of the 


error propagation with the effective thicknesses E and eT . In an experiment 


this would be the signal to measure. The wavelength contribution to the error was 
studied in a spectral range of (400-1000) nm for a dispersive element of tp g =5 mm 
and a sample of N-BK7 with tsmp=750um, Fig. 4.19 a). It is obvious that the 
error contribution of the wavelength to the group refractive index An E (A) is one 
order of magnitude larger than the influences of the sample and the DE thickness. 
Furthermore, it shows a significant minimum at a wavelength of à =683.75 nm. This 
behavior is similar to the sample-only approach of WPDE and is due to the spectral 


dependence of the refractive index of N-BK7 as sample and DE material. 
3.5e4 I T I I I 3.5e4 I I I I I 
3.0e4 tsmp = 750 um / | 3.0e4 tsmp = 750 pm / 
> tpe = 5 mm / — . tpg =5 mm / 
T 256% j AA =0.1 nm H a N-BK7 
< i 7 
<j 2.0e4 H I H H H H ea | a 2.0e4 H H F 
& 1.5e* j TP 1.56% 
& 1.06% h 1.004 
5.06" F 5,06" 
0 0 
0.4 0.5 0.6 07 08 09 1 0.4 0.5 0.6 07 0.8 09 1 
à [um] à [um] 
a) b) 


Figure 4.19 Spectral dependency in a range from 0.4-1 um with a) error contribution of the 
wavelength to the group refractive index Ang? E (X) for the WPDE approach with additional 
dispersion with tpg =5 mm, ftsmp =750 um and the measurement error for Aà =0.1 nm and 
b) total error on the calculated group refractive index of the sample with all discussed error 
contributions according to Eq. (4.27) 
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Finally, the total error of the group refractive index calculation for a 750 um thick 
sample of N-BK7 and a 5mm thick dispersive element of the same material has 
been determined as an absolute combination of all error contributions according 
to Eq.(4.27), Fig. 4.19. Due to the large influence of the error contribution of the 
refractive index measurement, the total error is dominated by it. Therefore, it shows 
the same characteristic minimum at A=683.75 nm with a value of 1.46 x 1075. It 
is important to note that the error increases exponentially for probing wavelengths 
smaller than the minimum with values larger than 3 x 1074, while the error slowly 
tends towards a boundary value for probing wavelengths above the minimum. The 
boundary value is significantly smaller than 5 x 1075. As this behavior is dominated 
by the refractive index of the materials present in the setup, an error estimation should 
be performed for every experiment. 

The newly developed approach for the characterization of cross-linking in poly- 
mers differs from state-of-the-art technologies in the accuracy of measurement as 
well as in the variety and quality of possible measurements. In comparison to clas- 
sical approaches such as Soxhlet extraction, DSC and DMA, the interferometric 
approach determines cross-linking in an indirect manner, but works non-destructive. 
This results in the ability to measure samples in in-line situations as well as during 
the lifetime of a product. It was shown that the typical measurement error was about 
0.1% in the scan-free approach. These errors are significantly lower than an error 
of 2-4% that classical approaches posses. Compared to other recently researched 
approaches, such as Raman or luminescence spectroscopy, the developed approach 
shows significant advantages relating its measurement speed where single profiles 
can be gathered in about 50 ms compared to multiple minutes. In summary, the 
developed approach combines high accuracy in the determination of cross-linking 
with high acquisition speeds and the ability to work in a scan-free fashion over lateral 
measurement ranges of multiple millimeters. Additionally, it allows non-destructive 
evaluations and combines material characterization with surface profilometry. In 
comparison with the identified state-of-the-art technologies for the determination 
of the degree of cross-linking, the DE-LCI approach shows significant advantages 
regarding accuracy, measurement time, inline capabilities and others, Tab. 4.3. 
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Table 4.3 Comparison of state-of-the-art technologies with the developed DE-LCI approach 
for the degree of cross-linking characterization in polymers 


approach meas. time | error spatial inline non- comments 
resolution destruc- 
tive 

Soxhlet 42h 2-4% no no no gold 

extraction standard 

(122, 124] 

DSC [145] 90 min. 10% no no no usable 
with many 
polymer 
types 

DMA [143, 2-6 h 10% no no no direct 

145, 149] determina- 
tion of 
degree of 
cross- 
linking 

spectroscopic |50-100s | 4-6% yes yes yes single 

[155, 156] point mea- 
surement, 
time 
consuming 

DE-LCI 50ms 0.01% yes yes yes indirect 
but very 
fast and 
precise 
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The third intended application for the proposed dispersion-encoded low-coherence 
interferometry is the evaluation of thin-film characteristics on substrate materials. 
Due to the usage of thin-film technologies in high-volume production in e.g. the 
photovoltaics and semiconductor industry, process monitoring becomes relevant in 
order to ensure functional parameters such as solar cell efficiency, [289]. In this 
context, film thickness as well as film homogeneity over large areas are important 
criteria for quality assurance. This section describes the modifications and devel- 
opments of the DE-LCI approach in order to measure film thickness of nm-sized 
films on comparatively thick substrates. Furthermore, results of different sample 
measurements are presented. 


5.1 Setup Considerations 


In order to achieve the goal of scan-free, spatially resolved film-thickness measure- 
ments of transmissive samples, a Mach-Zehnder interferometer was combined with 
an imaging spectrometer analogous to the approaches used in surface profilometry 
(see chapter 3) and polymer characterization (see chapter 4), Fig. 5.1. The collimated 
beam of a white-light source (AA =(400-1000)nm) was divided by a broadband 
plate beamsplitter into sample and reference arm. The reference arm contained a 
dispersive element (N-BK7, tpg =6.23 mm) while the sample arm was equipped 
with the particular sample which was a thin-film on a transparent substrate. If the 
sample contained a film thickness gradient, it was mounted in such a way that the 
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thickness gradient was aligned with the x-axis. The thickness and refractive index 
properties of the sample were dependent on the composition of the substrate mate- 
rial (ts, ns) and on the number of film layers (tfn, n fn) so that tsmp = f (tsub, tfn) 
and n°"? = f (Nsup, N fn). The transmitted light was superimposed with light from 
the reference arm after passing the secondary beamsplitter BS2 and imaged onto 
the slit of an imaging spectrometer. The used imaging magnification was typically 
M=0.6. The spectrally decomposed signal was detected with a two-dimensional 
CMOS-array of a camera. 

The signal formation and analysis were largely based on the theory of dispersion- 
controlled low-coherence interferometry established for precision profilometry and 
cross-linking characterization in polymers, chapter 3. However, due to the compo- 
sition of the sample as a multilayered system, a mathematically correct attribution 
has to be performed for all signal components to the sample composition. Typically, 
a transfer-matrix formalism is used to describe the transmission and reflection of 
light at every material interface, [290], Fig. 5.2. The basic idea of this formalism 
is that one singular matrix is used to describe the propagation of the electric field 
through a whole system of multiple layers. As of the design of the experiment 
within this work, only light with normal incidence regarding the front surface of 
the material system was taken into account. The formalism calculates the electric 


IMSPEC 


Figure 5.1 Experimental setup based on a Mach-Zehnder configuration with a WLS—white 
light source which is splitted 50:50 by BS 1—first beamsplitter into a reference beam which 
is directed by M1—treference arm mirror and manipulated by a DE—dispersive element 
(N-BK7, tpg = 6.23 mm) and the sample beam which is directed by M2—sample arm mir- 
ror through the SMP—sample with thin-film (thickness slope tf(x)) on a substrate. Both 
are recombined by BS2— second beamsplitter and imaged using the IL—imaging lens on a 
IMSPEC—imaging spectrometer where the analysis is performed 


M2 SMP BS2 IL 
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Figure 5.2 Simplified schema of the transfer-matrix formalism used to calculate the film 
thickness of a layer 1,5 on a substrate using the information of the present boundaries B;; and 
materials the electric field propagates through P;, [290] 


field components for a forward traveling wave (positive x-direction) E;„ and for a 
backwards traveling wave Eour in order to account for transmission and reflection 
at every material interface due to the material properties by using the wavenumber 
k}, and the propagation vector x 


Ein = Ge" 4 Met, Bay = Te ™ + Ret, (5.1) 


where the propagation coefficients 5, IT, Y, 2 hold information on the materials 
and interfaces. Considering the fact that the tangential components of the electric 
field must be continuous at a material interface as a boundary condition and that 
the matrices of every single layer have to be multiplied, a notation of the electric 
field component of every material and subsequent material boundary can be found 
as a matrix describing the electric field when entering, M;, and exiting the material 
system, Mo, 


eikx etikx = B eikx etik yY 5 2) 
-ike'kX +iketik® J (m) ~ sikei +iket® } | Q 


ln a ln 
M; Mo 


which leads to the final transfer matrix notation M7, 
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(5.3) 


i 


—1 -1 -1 
Mr = Mo; : Moo: Moo ` Mii * -M 5-10 ° 


M; 
The description of the electric field in this notation represents the resulting wave in 
the sample arm. In combination with a notation for the reference arm, the complete 
propagation of light through the Mach-Zehnder interferometer can be described. In 
combination with the estimated equalization wavelength, measured data sets could 
be fitted according to the methods mentioned before, see section 3.3. A detailed 
derivation of a simple, one-layer sample material which was used for analysis of the 
majority of measurements presented in this chapter, can be found in the appendix 
in the Electronic Supplementary Material (ESM). 

One important requirement in the practical implementation of this approach is 
the consideration of the substrate material as it shows incoherent behavior. Other- 
wise, the calculated signal inhibits higher order interferences which could reduce 
the signal-to-noise ratio significantly. Methods to minimize the influence of this 
problem would either be averaging the propagation in the substrate material by 
applying random phases or utilize a net-radiation method to get rid of the disturbing 
interference overlays, [291, 292]. 


5.2 Characterization of Thin-films on Bulk Substrates 


In order to evaluate the characteristics of the setup regarding its capabilities in thin- 
film characterization, samples of single-layer ITO coatings on polished float glass 
substrates (CEC020S, PGO GmbH, Germany) were prepared. The samples were 
half-sided chemical etched in a ferric chloride bath at 230K for 3 hrs. 

The determination of height profiles on substrates with coated thin-films can 
be performed by classical methods like spectral photometry or tactile profilometry 
either in a point wise or scanning fashion, Fig.5.3a). The acquisition of surface 
profiles using a tactile profilometer over a range of 4mm can be performed within 
about 30s, depending on the desired lateral resolution. The results show that the 
thickness gradient was captured with a mean height of 65.9 +12.3 nm. Addition- 
ally, some significant noise with a value of +3.6nm was observed in the data. In 
contrast to the scanning acquisition, the thickness gradient could be determined 
in a single acquisition along the lateral domain using the DE-LCI approach, Fig. 
5.3 b). In direct comparison, the captured profile, representing the average of ten 
successive data acquisitions (10x30 ms), shows significantly lower noise of about 
1.8nm. Accordingly, the thickness gradient over 4mm lateral measurement range 
was determined with (63.0 + 1.6) nm. The lateral resolution was about 4.2 um, deter- 
mined by the designed magnification of M=0.6. The significantly lower noise as 
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Figure 5.3 a) Plot of the thickness gradient, measured with a stylus profilometer as well as 
b) plot of the thickness gradient by DE-LCI 


well as the higher resolution and measurement speeds are clear advantages of the 
DE-LCI approach in thin-film characterization. 


5.3 Characterization of Flexible Substrate Materials 


As described above, the transfer-matrix formalism is used as a mathematical basis 
for the fitting of experimentally acquired spectra. The influence of the substrate 
material parameters is usually low when analyzing thin-films on bulk substrates. 
Some applications of thin-films are fabricated on thin, flexible substrates such as 
polymer sheets, [289]. Thus, a more accurate knowledge of the material properties, 
especially of the thickness of these substrate materials, is necessary. In order to 
account for this, a dual-channel variation of the originally developed Mach-Zehnder 
interferometer was developed, Fig. 5.4. 

In this setup, light is guided as before, see Fig. 5.1. Additionally, light from 
all interfaces but specifically that from the substrate was collected in a secondary 
detection path at the first beamsplitter. This path was equipped with a high-resolution 
grating spectrometer (Ospec = 0.3 nm, Avaspec ULS3648 VB, Avantes BV, Apel- 
doorn, The Netherlands). In this context, the surfaces of the substrate material can 
be considered as a resonator where interference occurs dependent on the distance 
of the surfaces, hence the thickness. By utilizing the spectrometer to record the 
spectral modulations of this interference signal, the actual thickness was calculated 
according to the principles of Fourier-domain optical coherence tomography, [105]. 
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Sheets of polyethylene terephthalate (PET) foil substrate (ts = 135 um) with 
a coated ITO layer (nominal thickness og = 150 nm), typically applied in the 
photovoltaics industry, were used as samples. The ITO coating was partly removed 
from the sample by means of chemical etching to generate a film-thickness gradient 
for spatial investigations. From the measured back-reflected interference signal, a 
Fourier analysis with appropriate x-axis scaling could be performed, Fig. 5.5 a). The 
data shows a significant amount of DC-signal components and noise, but also very 
clear peaks. These peaks could be attributed to the first and second order reflec- 
tions of the substrate surfaces. The peak information was extracted by applying a 
Blackman-Harris window function and fitted with a Gaussian function. From the 
fitted data, the optical path distances of the reflections could be estimated with 
(450.4 and 443.2) um respectively. With the knowledge of the refractive index of 
the substrate material, the thickness was calculated for both reflections with (137.3 
and 135.1) um. The thickness of the substrate material could be confirmed by mea- 
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Figure 5.4 Dual-channel setup based on a Mach-Zehnder configuration with a WLS— white 
light source which is splitted 50:50 by BS1—first beamsplitter into a reference beam which 
is directed by M1—reference arm mirror and manipulated by a DE—dispersive element 
(tpg = 6.23 mm) and the sample beam which is directed by M2—sample arm mirror through 
the SMP—sample with thin-film (thickness slope t ¢(x)) on a substrate (ts = 135 um) and both 
are recombined by BS2—second beam splitter. The analysis is performed using a CH1— 
IMSPEC—primary, imaging spectrometer on which the recombined beam is imaged to with 
a LE—lens, here marked in blue as well as with the CH2—SPEC—secondary spectrometer 
which records the interference of the back-reflected signal from the substrate, here marked in 
red 
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surements on a tactile profilometer. With this result, an appropriate start value for 
the calculation of the film thickness with the imaging spectrometer is given. The 
method only allows the thickness calculation of the substrate material as the thin- 
film generates high-frequent interference which is not resolvable with a standard 
spectrometer. 
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Figure 5.5 Results of film thickness measurements using a dual-channel interferometer with 
a) Fourier-analyzed data of the secondary channel with the x-axis scaled to the optical path 
difference (OPD) showing peaks from interference of two reflections which were used to 
calculate the substrate thickness as 137.3 and 135.1 um respectively and b) slope of the ITO 
coating on this substrate measured in the primary channel having a ITO thickness of 151.6 nm 
in comparison to the slope measured on a tactile profilometer where the height was measured 
with 152.4 nm 


By making use of the measured substrate thickness, data from the transmission 
measurement could be analyzed with the transfer-matrix formalism and appropriate 
fitting according to the model described with Eq. (3.4) and (3.6), Fig. 5.5 b). The 
results show that the slope of an ITO film on a PET substrate can be sampled with 
the appropriate resolution. From this slope the film thickness was measured with 
151.6 nm. This corresponds well with a measurement of the sample with a tactile 
profilometer (Talysurf i-Series, Taylor Hobson Ltd, UK) showing a measured thick- 
ness of 152.4 nm. In the data, it is noticeable that the tactile profilometer shows 
significant, periodic noise which can be attributed to the surface roughness of the 
substrate material. In direct comparison the surface roughness is less pronounced 
in the interferometric data. Furthermore, the interferometric system was capable of 
resolving the slope as well as the film thickness with high precision. This simulta- 
neous measurement helped to improve the accuracy of the fit model for the thin-film 
thickness determination without compromising the actual measurement. 
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While ellipsometric approaches can be error-prone on e.g. flexible substrates, 
the demonstrated approach is capable to measure film thickness in this setting 
with high accuracy. In comparison to spectral reflectometry, another state-of-the-art 
technology, the developed approach enables comparable resolutions on film thick- 
ness measurements while also maintaining a far larger measurement range of about 
80 um. This enables the simultaneous capture of the film thickness as well as the 
thickness of substrate materials. In summary, the DE-LCI approach shows signifi- 
cant advantages regarding accuracy, measurement time and the variety of possible 
samples, Tab. 5.1. 


Table 5.1 Comparison of state-of-the-art technologies with the developed DE-LCI approach 
for thin-film characterization 


approach meas. time | resolution | spatial inline flexible comments 
[nm] resolution substrates 


reflectometry | 60s 10nm no yes yes single 
[179] point mea- 
surements, 
large 
integration 
times 
necessary 


ellipsometry | 8s/A 0.01 nm yes yes no gold 
[211] standard, 
high 
accuracy 


DE-LCI 50ms 0.1 nm yes yes yes combines 
advantages 
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Conclusion 6 


The characterization of features such as geometrical dimensions or functional 
parameters like roughness or degree of cross-linking plays an important role dur- 
ing production in industries such as semiconductors, organic-electronics and the 
photovoltaics industry. The ability to gather measurements in a process-integrated 
fashion is of high interest. 

The scope of this work, in this context, was the development of an optical metrol- 
ogy tool which is capable to provide information on surface features and material 
properties in a fast and versatile way. In contrast to existing technologies, a novel 
approach was developed which delivered information on surface line profiles with- 
out the need for mechanical scanning and with higher resolution. Furthermore, 
in some aspects the approach was also able to gather information which was not 
accessible in a spatially-resolved fashion before. The developed setup was based 
on a modified, dispersion-enhanced low-coherence interferometer (DE-LCI). The 
properties of interest were characterized in the three major applications surface 
profilometry, polymer and thin-film characterization. 

In the main part of this work, the DE-LCI approach was designed and imple- 
mented for surface profilometry. Here, it was shown how controlled dispersion can 
be used to encode path length differences in an interferometer and therefore sur- 
face height information in the spectral domain. The dispersion was controlled by a 
dispersive element, i.e. a glass window, which was used as one possibility to adjust 
the axial measurement range of the setup. Additionally, it was demonstrated how 
an appropriately designed imaging spectrometer can be utilized to gather the sur- 
face height information along a line profile in a single data acquisition without the 
need for mechanical scanning. The introduction of an additional, movable lens after 
the interferometric recombination extended the setups capabilities and allowed the 
acquisition of three-dimensional surface height information. The development of a 
custom data analysis and fitting routine led to an estimation of the axial measure- 
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ment range of Az=79.91 um while a theoretical axial resolution of 0.088 nm was 
calculated for a selected experimental configuration. With these values, a dynamic 
range in the axial dimension of DR=9.0 x 10° was estimated. 

The characterization of a setup with the calculated properties was performed 
by analyzing the surface profiles of different measured step height standards. Most 
notably, the determination of the profile and the height of (101.8+0.1)nm on a 
silicon height standard demonstrated the capabilities of the setup. It could be shown 
that the setup has a repeatability of o7=0.13nm while the axial resolution was 
found to be Az jn =0.1nm. In relation to the available measurement range of 
Az=79.91 um, the experimentally determined dynamic range was DR=7.99 x 10°. 
This value is about 6 times higher than comparable current approaches known from 
literature, [119]. 

Furthermore, it could be demonstrated that functional parameters such as sur- 
face roughness are measurable with the same axial resolution like surface profiles 
although this data was gathered over a lateral measurement range of up to 1.5mm. 
This separates the novel approach distinctively from established technologies such 
as tactile profilometry or confocal laser scanning microscopy which rely on time 
consuming and error-prone methods of scanning or stitching to enable lateral mea- 
surement ranges of the same order. 

The capabilities of the setup regarding the acquisition of three-dimensional sur- 
face height information were evaluated with a measurement of a um-sized, PTB- 
calibrated height standard. By imaging an area of 1.5 x 0.25 mm? with sub-nm resolu- 
tion, the measurement of steps with heights of (971.26 + 0.31), (4951.40 + 0.28) and 
(19924.00 + 0.36) nm was performed while additional features such as the rough- 
ness of each step could be acquired simultaneously. 

Anextension of the setup with an imaging spectrometer for the NIR spectral range 
(AA =(1133-1251)nm) enabled tomographic measurements of a silicon sample. 
Specifically, the front and back surface of a thinned wafer could be investigated. 

In summary, the developed dispersion-encoded approach to interferometry for 
surface profilometry proved to be of high resolution in the axial dimension while cap- 
turing large measurement ranges in the lateral dimension. Due to the high-dynamic 
range and fast data processing, an application can be envisioned in process-integrated 
metrology for industrial production. 

The determination of the degree of cross-linking is an important criterion in 
quality assurance of polymer processing. The degree of cross-linking determines 
important mechanical properties of the fabricated products as well as their long-term 
durability. Additionally, the cross-linking process bears potential for optimization 
during production in regard to speed and properties. Therefore, production accom- 
panying monitoring is desirable. 
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Based on these conditions, the DE-LCI approach was adapted accordingly and 
tested for its ability to characterize polymers. The characterization was based on the 
measurement of the wavelength-dependent refractive index n(A) as a measure for the 
degree of cross-linking. The DE-LCI approach was utilized in a temporal scanning 
as well as in a scan-free configuration. Using the temporal scanning configuration, 
n(A) could be analyzed over large spectral ranges of AA =(400- 1000)nm while 
a resolution in terms of the group refractive index of o5(ng)=2.13 x 1074 was 
achieved. The measurement took several seconds and no spatial resolution could be 
accomplished. In contrast, the scan-free configuration was capable to measure the 
refractive index with a resolution of 3.36 x 1075 on a profile of 250 um length in 
50 ms. The spectral range was about 20 nm. In context of the scan-free configuration, 
a novel mathematical method for the analysis of phase data based on wrapped-phase 
derivative evaluation (WPDE) was developed and qualified. Both configurations 
were tested using typical samples from the photovoltaics and semiconductor industry 
respectively. 

In contrast to existing technologies to determine the degree of cross-linking, the 
novel approach is fast, non-destructive and capable to perform spatially-resolved 
measurements with a lateral resolution of about 5 um. The combination of charac- 
teristics is unique to the novel approach. 

The precise control of single-layer thickness in the production of thin-film sys- 
tems is crucial for the performance of these systems e.g. in the organic electronics 
or photovoltaics industry. 

The DE-LCI setup was adapted as a Mach-Zehnder type interferometer in order 
to characterize the thickness of single-layer thin-films. By measuring a layer of 
ITO on a bulk glass substrate, it was demonstrated that the layer thickness can be 
measured with a resolution of 1.6nm. The modification of the setup by using a 
secondary spectrometric detection channel allowed to capture back-reflected light 
of the sample. This way, it was possible to in-situ evaluate the substrate thickness 
of a flexible substrate (tsup = 135um) while measuring the film thickness of a 
layer of ITO (t779 =151.6nm) simultaneously. This improved the robustness of the 
underlying transfer-matrix model which was utilized to calculate the film thickness 
from recorded spectral data. 

The ability to capture film thickness data with spatial resolution within a single 
frame acquisition separates the DE-LCI approach from existing technologies and 
determines its usability as a process-integratable tool in future applications. 

As demonstrated, the scope of this work was to show a range of possible applica- 
tions of the developed DE-LCI approach while examining each application in depth. 
While answering the major questions by experiments and analysis, some additional 
questions arose which could not be addressed within the scope of this work. Most 
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notably, alternative approaches might be addressed in future works in order to ana- 
lyze the captured data. Preliminary experiments have shown that methods of image 
processing and statistical analysis are promising with regard to resolution, speed and 
robustness. Furthermore, it can be noted that the analysis of the captured data by 
means of FD-OCT is possible, [98]. Here, the different dispersion in the interferom- 
eter arms has to be compensated by multiplying the spectral data with an appropriate 
phase function, [258]. After re-sampling and Fourier transformation the data con- 
tains three components: the desired peak in depth space, second the broadened 
mirror term and third the broadened DC term. Due to the signal construction in an 
DE-LCI approach, these terms overlap and need more sophisticated compensation 
in order to perform analysis known from full-range OCT approaches, [111]. Never- 
theless, this approach bears potential for future work as it will increase the possible 
axial measurement range of the approach. Additionally, the development of fitting 
routines which can deal with data where the equalization wavelength lies not within 
the spectral range of the detector are interesting to enhance with respect to the axial 
measurement range. Beside these approaches to data analysis, advanced design and 
construction efforts with regard to mechanical stability and thermal management are 
interesting to further improvements of e.g. the repeatability of the setup. In terms of 
hardware aspects, the utilization of the triggerable high-power NIR supercontinuum 
light source makes possible the characterization of dynamic surface profilometry of 
oscillating samples in a stroboscopic acquisition mode. The experimental validation 
of theoretical designed methods to gather three-dimensional surface information in 
a scan-free fashion is of high interest as well. 

Furthermore, the characterization of cross-linking mechanisms during actual 
processes with high temporal resolution as well as the monitoring of internal stresses 
during cross-linking by polarization enhanced DE-LCI is an interesting topic for 
future research. 

The extension of the data analysis in thin-film characterizations towards multi- 
layered systems can be envisioned. 

The combination of advantageous features such as the versatility of applications, 
the high dynamic range in the axial as well as the lateral domain, the capabilities 
to tune the measurement ranges easily and the high resolution characterize the 
developed approach. In conclusion, it can be said that the DE-LCI approach was 
developed and qualified in its main features and bears the potential for interesting 
applications in research as well as in industry. 
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